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Abstract 

Unwanted  noise  is  a  by-product  of  many  industrial  processes  and  systems.  In 
active  noise  cancellation  (ANC),  one  introduces  a  secondary  noise  source  to  generate 
an  acoustic  field  that  interferes  destructively  with  the  unwanted  noise,  and  thereby 
attenuates  it.  Noise  reduction  is  important  to  protect  listeners  in  high  noise  envi¬ 
ronments  from  hearing  damage,  to  enhance  speech  communication,  and  to  reduce 
noise-induced  fatigue.  These  adverse  effects  of  noise  can  cause  accidents  and  reduce 
the  productivity  of  workers. 

By  formulating  the  ANC  problem  as  an  optimal  feedback  control  problem,  we 
developed  a  single  approach  for  designing  stochastically  optimal  feedback  controllers 
for  pointwise  and  distributed  active  noise  cancellation.  In  the  pointwise  case,  we 
develop  an  optimal  feedback  controller  for  attenuating  the  acoustic  pressure  at  the 
locations  of  a  finite  number  of  microphones.  In  the  distributed  case,  we  develop  an 
optimal  feedback  controller  for  attenuating  the  total  acoustic  energy  in  an  enclosure. 
In  either  case,  the  control  signal  is  pointwise  and  drives  an  ordinary  loudspeaker,  and 
the  input  of  the  controller  is  pointwise  and  is  obtained  as  the  output  of  an  ordinary 
pressure  microphone. 

Thesis  Supervisor:  Alan  V.  Oppenheim 

Title:  Distinguished  Professor  of  Electrical  Engineering 
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Chapter  1 


Introduction 


Unwanted  noise  is  a  by-product  of  many  industrial  processes  and  systems.  In  active  noise 
cancellation  (ANC),  one  introduces  a  secondary  noise  source  to  generate  an  acoustic  field 
that  interferes  destructively  with  the  unwanted  noise,  and  thereby  attenuates  it. 

Noise  reduction  is  important  to  protect  listeners  in  high  noise  environments  from  hearing 
damage,  to  enhance  speech  communication,  and  to  reduce  noise-induced  fatigue.  These 
adverse  effects  of  noise  can  cause  accidents  and  reduce  the  productivity  of  workers. 

Passive  silencers  are  low-pass  acoustical  filters  commonly  used  on  engines,  blowers,  com¬ 
pressor,  fans  and  other  industrial  equipment.  These  silencers  are  very  effective  in  attenuat¬ 
ing  high  frequency  noise;  however,  they  are  usually  bulky  and  ineffective  at  low  frequencies. 
This  is  because  at  low  frequencies  the  acoustic  wavelengths  become  large  compared  to  the 
thickness  of  a  typical  passive  silencer.  A  sound  wave  of  60Hz,  for  example,  has  a  wavelength 
of  5.66  meters  in  air  under  normal  conditions.  It  is  also  difficult  to  stop  low  frequency  sound 
being  transmitted  from  one  space  to  another  unless  the  intervening  barrier  is  very  heavy. 
Passive  mufflers  that  are  used  to  silence  engine  noise  produce  back-pressure  by  obstruct¬ 
ing  the  turbulent  flow  out  of  the  engine,  and  this  back-pressure  significantly  reduces  the 
efficiency  of  the  engine.  To  overcome  these  problems,  researchers  have  given  active  noise 
cancellation  considerable  attention  recently. 

The  fundamental  problem  in  ANC  is  to  generate  an  acoustic  field  that  interferes  de¬ 
structively  with  the  unwanted  noise  field  at  the  points  of  interest.  A  typical  ANC  system 
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utilizes  several  microphones  to  monitor  the  attenuated  field  and  several  canceling  sources  to 
generate  the  canceling  field.  The  outputs  of  these  microphones  are  used  as  inputs  to  some 
sort  of  electronic  controller.  This  controller  provides  the  inputs  to  the  canceling  sources  in 
such  a  way  that  the  acoustic  field  generated  by  these  sources  interferes  destructively  with 
the  unwanted  noise  field  at  the  points  of  interest. 

Almost  all  existing  ANC  systems  are  designed  explicitly  to  minimize  the  sum  of  the  noise 
power  at  a  finite  number  of  spatial  points.  We  shall  refer  to  this  type  of  performance  criterion 
as  a  “pointwise  performance  criterion”.  Similarly,  we  refer  to  the  resulting  ANC  system  as 
a  “pointwise  active  noise  cancellation  system”.  In  pointwise  ANC,  it  is  commonly  assumed 
that  attenuating  the  noise  at  a  finite  number  of  spatial  points  will  results  in  attenuation  of 
the  noise  at  aU  other  points  of  interest.  However,  it  has  been  shown  in  numerical  studies 
that  attenuation  of  noise  at  a  finite  number  of  spatial  points  can  result  in  amplification  of 
the  noise  at  other  points  [15]. 

Existing  controllers  for  pointwise  ANC  can  be  characterized  as  feedforward  or  feedback 
controllers.  Adaptive  feedforward  controllers  have  been  developed  by  a  number  of  authors, 
and  these  controllers  are  designed  to  take  advantage  of  the  statistics  of  the  signals  involved. 
However,  most  existing  feedback  ANC  systems  are  designed  deterministically  and  do  not 
take  advantage  of  the  statistical  characteristics  of  the  noise. 

A  distributed  active  noise  cancellation  system  is  defined  as  an  ANC  system  that  mini¬ 
mizes  a  distributed  performance  criterion.  For  example,  a  distributed  ANC  system  might 
minimize  the  total  acoustic  energy  in  an  enclosure.  To  the  best  of  our  knowledge,  a  mathe¬ 
matically  rigorous  procedure  for  designing  distributed  ANC  systems  has  not  been  developed 
yet. 

By  formulating  the  ANC  problem  as  an  optimal  feedback  control  problem,  we  develop 
a  single  approach  for  designing  both  pointwise  and  distributed  ANC  systems.  The  key 
strategy  is  to  model  the  residual  signal/field  as  the  sum  of  the  outputs  of  two  linear  systems. 
The  unwanted  noise  signal/field  is  modeled  as  the  output  of  a  linear  system  driven  by  a 
white  process.  Similarly,  the  canceling  signal/field  is  modeled  as  the  output  of  another 
linear  system  driven  by  the  control  signal.  Finally,  the  residual  signal/field  is  modeled  as 
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the  sum  of  the  outputs  of  these  two  linear  systems.  We  show  that  the  control  signal  that 
minimizes  a  certain  class  of  performance  criteria  is  a  linear  feedback  of  the  estimated  states 
of  these  two  linear  systems.  These  state  estimates  are  computed  using  a  Kalman  filter,  and 
the  feedback  gain  matrix/operator  is  obtained  by  iterating  a  Riccati  equation.  Note  that 
in  the  pointwise  case  and  in  the  distributed  case,  the  control  signal  is  not  distributed  and  is 
used  as  input  to  an  ordinary  loudspeaker.  Moreover  in  both  these  case,  the  control  signal 
is  generated  based  on  the  outputs  of  ordinary  microphones  that  monitor  the  residual  field. 

While  we  focus  specifically  on  the  acoustic  noise  cancellation  problem,  the  results  devel¬ 
oped  in  this  thesis  can  be  applied  to  other  active  cancellation  problems.  Vibration  control 
is  an  example  of  a  non-acoustic  problem  to  which  our  results  can  be  applied. 

1.1  Outline  of  the  Thesis 

Chapter  2  motivates  our  study  of  active  noise  cancellation  and  presents  a  summary  of 
the  previous  work  done  in  this  area.  The  physics  of  the  problem  are  discussed  in  this 
chapter  with  emphasis  on  three  dimensional  acoustic  fields.  The  different  sound  fields 
generated  by  ANC  systems  designed  to  achieve  different  acoustic  objectives  are  studied. 
Several  commonly  used  control  strategies  for  ANC  are  presented,  and  the  performance  of 
these  control  strategies  are  analyzed  and  compared.  The  relationship  between  active  noise 
cancellation  and  signal  estimation  is  also  discussed. 

The  optimal  feedback  controller  for  attenuating  the  acoustic  pressure  at  a  discrete  num¬ 
ber  of  microphone  locations  is  developed  in  Chapter  3.  Our  approach  is  to  model  the 
residual  signals  at  the  microphones  as  the  sum  of  the  outputs  of  two  finite-dimensional 
linear  systems.  The  unwanted  noise  signals  at  the  microphones  are  modeled  as  the  outputs 
of  a  multiple-input /multiple-output  linear  system  driven  by  a  white  process.  Similarly,  the 
canceling  signals  at  the  microphones  are  modeled  as  the  outputs  of  a  single-input/multiple- 
output  linear  system  driven  by  the  control  signal.  Finally,  the  residual  signals  at  the  micro¬ 
phones  are  modeled  as  the  sum  of  the  outputs  of  these  two  linear  systems.  We  show  that 
the  control  signal  that  minimizes  a  certain  class  of  performance  criteria  is  a  linear  feedback 
of  the  estimated  states  of  these  two  linear  systems.  We  also  show  that  our  formulation  can 
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be  used  to  minimize  the  frequency  weighted  power  of  the  residual  signal.  Such  frequency 
dependent  weighting  might  be  important  because  of  the  difference  in  the  sensitivity  of  the 
ear  to  sounds  at  different  frequencies.  These  state  estimates  can  be  computed  using  a 
Kalman  filter,  and  the  feedback  gain  matrix  can  be  computed  by  iterating  a  matrix  Riccati 
equation.  The  single-microphone  case  is  presented  first,  and  the  results  are  then  extended 
to  the  multiple-microphone  case. 

In  Chapter  4,  several  recursive/adaptive  algorithms  for  modeling  the  unwanted  noise  at 
a  single  microphone  are  developed.  We  assume  that  the  unwanted  noise  at  the  microphone 
is  the  output  of  an  all-pole  transfer  function  driven  by  a  white  process  and  develop  a 
recursive/adaptive  algorithm  for  estimating  the  parameters  of  this  transfer  function  based 
on  the  measurements  made  by  the  microphone.  An  adaptive  feedback  ANC  system  can  be 
obtained  by  combining  the  optimal  control  law  of  Chapter  3  with  these  recursive/adaptive 
algorithms  for  modeling  the  unwanted  noise.  Although  we  focus  specifically  on  developing 
estimation  algorithms  for  modeling  the  unwanted  noise  in  the  ANC  problem,  the  algorithms 
presented  in  this  chapter  can  be  applied  to  the  more  general  problem  of  identifying  non¬ 
stationary  autoregressive  (AR)  processes  embedded  in  white  noise. 

Chapter  5  develops  a  distributed  optimal  feedback  controller  for  minimizing  the  total 
acoustic  energy  in  an  enclosure.  Our  approach  is  to  model  the  residual  field  as  the  sum  of  the 
outputs  of  two  infinite-dimensional  linear  systems.  The  unwanted  noise  field  is  modeled  as 
the  output  of  an  infinite-dimensional  linear  system  driven  by  a  white  process.  Similarly,  the 
canceling  field  is  modeled  as  the  output  of  another  infinite-dimensional  linear  system  driven 
by  the  control  signal.  Finally,  the  residual  field  is  modeled  as  the  sum  of  the  outputs  of  these 
two  linear  systems.  The  inputs  to  the  controller  are  the  outputs  of  ordinary  microphones, 
and  the  control  signal  drives  an  ordinary  loudspeaker.  We  show  that  the  control  signal 
that  minimizes  a  certain  class  of  performance  criteria  is  a  linear  feedback  of  the  estimated 
states  of  these  two  linear  systems.  These  state  estimates  can  be  computed  using  a  Kalman 
filter,  and  the  feedback  gain  operator  can  be  computed  by  iterating  an  operator  Riccati 
equation.  This  is  the  first  mathematically  rigorous  formulation  of  the  distributed  ANC 
problem  known  to  us. 
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Lastly,  in  Chapter  6  we  summarize  the  major  contributions  of  this  thesis.  We  also 
suggest  several  potentially  important  directions  for  future  research. 
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Chapter  2 


Previous  Work  on  Active  Noise 
Cancellation 


2.1  Introduction 

This  chapter  motivates  our  study  of  active  noise  cancellation  and  presents  a  summary  of 
the  previous  work  done  in  this  area.  We  try  to  distinguish  between  the  acoustic  objectives 
of  different  ANC  systems  and  the  control  strategies  used  to  achieve  these  objectives.  To 
appreciate  the  advantages  and  the  limitations  of  active  noise  cancellation,  it  is  necessary 
to  understand  the  relevant  acoustical  principles  as  well  as  the  relevant  control  strategies. 
Hence,  this  chapter  includes  a  Section  on  acoustical  principles  behind  ANC  and  a  Section 
on  control  strategies  used  for  ANC.  The  relationship  between  active  noise  cancellation  and 
signal  estimation  is  also  studied  in  this  chapter. 

This  chapter  is  organized  as  follows.  Section  2.2  describes  the  physical  basis  for  active 
noise  cancellation,  with  emphasis  on  three  dimensional  sound  fields.  The  different  sound 
fields  that  are  generated  by  ANC  systems  designed  to  achieve  different  acoustic  objectives 
are  studied  in  this  Section  using  Modal  Analysis  techniques.  In  Section  2.3,  we  review  several 
commonly  used  control  strategies  for  ANC.  Many  examples  of  existing  ANC  systems  are 
presented  to  illustrate  the  use  of  these  control  strategies.  In  Section  2.4  ,  the  performance 
of  two  specific  control  strategies  for  ANC  is  analyzed  in  detail.  Finally  in  Section  2.5,  the 
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relationship  between  active  noise  cancellation  and  signal  estimation  is  explored.  A  brief 
summary  of  this  chapter  is  presented  in  Section  2.6. 

2.2  Acoustical  Principles 

The  fundamental  problem  in  active  noise  cancellation  is  to  generate  an  acoustic  field  that 
interferes  destructively  with  the  undesired  noise  [10,  13,  23,  40,  50,  59].  The  undesired  field 
is  usually  referred  to  <ls  the  “primary”  field,  and  the  interfering  field  is  referred  to  as  the 
“secondary”  field.  The  question  therefore  arises  as  to  whether  two  different  sources  of  sound 
can  generate  the  same  acoustic  field?  If  the  answer  to  this  question  is  yes  and  one  source  is 
under  our  control,  a  simple  change  of  sign  of  the  controlled  source  will  cancel  the  primary 
field  everywhere. 

We  shall  give  a  simple  example  of  two  different  sources  that  generate  the  same  acoustic 
field.  Let  q(x,t)  denote  a  sound  source  generating  the  pressure  field  where  x  is 

the  spatial  variable;  t  is  the  temporal  variable;  and  the  underline  denotes  a  vector-valued 
quantity.  Furthermore,  assume  that  the  source  is  only  non-zero  over  some  bounded  spatial 
domain  F.  In  this  case,  the  wave  equation  for  pressure  is 

=  (2.1) 

where  is  the  three  dimensional  Laplacian,  and  c  is  speed  of  sound.  The  pressure  field 
p{x,  t)  outside  the  region  of  support  of  the  source  will  remain  unchanged  if  the  source  field 
is  supplemented  by  j^-^q{x,t)  —  V^q{x,t).  To  see  this,  note  that  the  resulting  pressure 
field  in  this  case,  p'{x,t),  must  satisfy 

^■^p'ix,t)  -  VV(®,0  =  (2-2) 

and  eq.  (2.2)  can  be  rewritten  as 

^^{p'{x,t)  -  q{x,t))  -  V^{p'{x,t)  -  q{x,t))  =  q{x,t).  (2.3) 

Comparing  (2.3)  to  (2.1)  and  recalling  that  q(x,t)  is  zero  outside  T,  we  conclude  that 
p'(x,  t)  and  p{x,  t)  are  equal  outside  F.  Therefore,  we  have  found  an  example  of  two  distinct 
sources  that  produce  the  same  acoustic  pressure  field. 
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More  generally,  Kempton  [31]  has  shown  that  the  sound  field  generated  by  any  source 
can  be  reproduced  by  an  appropriate  infinite  series  of  source  singularities  positioned  at 
any  desired  point.  Hence,  complete  cancellation  of  the  primary  field  can  be  achieved  by 
arranging  these  source  singularities  to  produce  the  negative  of  the  primary  field. 

A  more  obvious  way  to  cancel  the  effect  of  an  acoustic  source  is  to  use  Kirchhoff’s 
theorem.  This  theorem  provides  a  formula  by  which  the  entire  effect  of  sources  inside  a 
closed  boundary  can  be  duplicated  by  sources  on  this  boundary.  Hence,  if  the  negative  of 
these  sources  is  used  on  the  boundary,  the  primary  field  in  the  exterior  of  the  boundary 
is  completely  canceled.  Jessel  and  Angevine  [28,  29]  have  proposed  such  an  ANC  system 
in  theory.  Unfortunately,  the  secondary  source  in  this  case  is  distributed,  and  a  discrete 
approximation  to  this  distributed  source  has  not  been  developed  yet.  Inspired  by  this  theory, 
several  researchers  have  built  experimental  ANC  systems  with  limited  success  [4,  32]. 

With  active  noise  cancellation,  the  combined  radiated  output  of  the  primary  and  sec¬ 
ondary  sources  is  typically  much  less  than  the  radiated  output  of  the  primary  source  alone. 
The  reduction  in  the  combined  radiated  power  is  achieved  by  a  reduction  in  the  radiation 
impedance  seen  by  the  primary  source  and/or  absorption  of  energy  by  the  secondary  source 
[17,  46]. 

Because  of  the  practical  difficulties  associated  with  cancellation  of  an  entire  noise  field, 
most  researchers  have  focused  on  developing  ANC  systems  with  very  limited  acoustic  ob¬ 
jectives.  The  most  commonly  used  acoustic  objective  is  attenuation  of  pressure  at  a  finite 
number  of  spatial  points. 

Modal  analysis  is  one  of  the  very  few  techniques  avaiable  for  studying  the  global  effects 
of  ANC  in  an  enclosed  volume. 

Modal  Analysis  for  ANC 

Modal  analysis  is  a  practical  approach  for  analyzing  the  sound  field  generated  by  an  acoustic 
system  operating  in  an  enclosure.  The  starting  point  is  to  assume  that  the  sound  field  in  the 
enclosure  has  a  periodic  time  dependence  .  Let  P(x,  t)  refer  to  the  pressure  at  location 
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X  and  at  time  t.  In  this  case,  P{x_,t)  can  be  expressed  as 

P{x,  t)  =  p{x,  (2.4) 

where  p(x,  fi)  is  not  a  function  of  time.  It  is  further  assumed  that  p(x,  fl)  can  be  expressed 
in  terms  of  a  finite  number  of  the  normal  modes  of  the  enclosure 

N 

p(x,ft)  =  ^  $„(x)a,j(n),  (2.5) 

n=0 

where  $n(2.)  is  the  the  n-th  normal  mode  of  the  enclosure,  and  an(n)  is  the  complex 
amplitude  of  the  n-th  mode. 

Next,  it  is  assumed  that  the  sound  field  is  linear  so  that  the  complex  amplitude  a„(fi) 
can  be  expressed  as  the  sum  of  the  contributions  from  the  primary  source  and  M  secondary 
sources,  i.e. 

M 

c„(Q)  =  aP(fi)  -t-  Yi  5„„,(fi)g^(n),  (2.6) 

m=l 

where  is  the  n-th  modal  amplitude  produced  by  the  primary  source;  is  the  net 

strength  of  the  m-th  secondary  source;  and  specifies  the  degree  to  which  the  m-th 

secondary  source  contributes  to  the  n-th  mode. 

Equation  (2.6)  can  be  rewritten  in  vector  form  as 

a  =  flP  -1-  Bq^,  (2.7) 

where  a  is  the  iV  x  1  vector 

ai(fi) 

Ct2(fi) 

£=  ;  (2.8) 

flAr(fl) 
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is  the  N  X  1  vector 


a?(n) 


is  the  M  x  1  vector 


qm 


(2.10) 


L  qhi^)  J 

and  B  is  the  N  x  M  matrix  of  complex  modal  coupling  coefficients 


Bii(n)  ... 


(2.11) 


Br^i{Q)  . . .  Bnm{^)  J 

The  total  time-averaged  acoustic  potential  energy  in  the  enclosure  is  denoted  by  Ep  and 
is  given  by 

Ep  -  bfe  fi)pdx,  (2.12) 

where  c  is  the  speed  of  sound,  and  p  is  the  density  of  fluid  in  the  enclosure.  Recalling  the 
orthonormality  of  'J'n’s,  we  see  that 


1  n 
4^^-  - 


(2.13) 


(2.14) 


Using  the  above  formulation,  one  can  find  the  secondary  source  strength  that  min¬ 
imizes  the  acoustic  potential  energy  Ep.  Since  Ep  is  a  quadratic  function  of  the  source 


strength  vector  ,  the  following  source  strength  vector  minimizes  the  total  acoustic  poten¬ 
tial  energy  in  the  enclosure 

£*  =  (2.15) 

and  the  resulting  minimum  value  of  the  potential  energy  is 

Ej,  =  ^[(aP)V  -  iaPfB[B^B]-^B”aP].  (2.16) 

A  similar  formulation  can  be  used  to  express  the  total  acoustic  kinetic  energy  in  the 
enclosure,  Ek,  &s  a.  quadratic  function  of  the  source  strength  vector  q‘.  Furthermore,  one 
can  find  the  secondary  source  strength  that  minimizes  Ek  or  (Ek  +  Ep).  Note  that  (Ek  +  Ep) 
is  the  total  acoustic  energy  in  the  enclosure. 

Nelson  et  al.  [36]  considered  the  problem  of  minimizing  the  acoustic  potential  energy  in 
a  rectangular  enclosure.  The  three  dimensional  sound  field  in  the  enclosure  was  expressed 
as  the  sum  of  the  contributions  of  7000  modes.  Given  the  position  and  the  strength  of  a 
pure  tone  primary  source  (a  point  source),  the  authors  used  equation  (2.15)  to  find  the 
amplitude  and  phase  of  a  number  of  secondary  point  sources  that  would  minimize  the  total 
acoustic  potential  energy  in  the  enclosure. 

At  frequencies  above  the  Schroeder  cut-off  frequency  of  the  enclosure,  the  authors  found 
that  substantial  reductions  in  acoustic  potential  energy  are  not  obtainable  unless  the  can¬ 
celing  sources  are  separated  from  the  primary  source  by  a  distance  of  no  more  than  half  the 
wave  length  [9]. 

At  frequencies  below  the  Schroeder  cut-off  frequency  of  the  enclosure,  the  system  exhib¬ 
ited  several  interesting  features  [19].  The  authors  found  that  appreciable  reductions  in  the 
overall  potential  energy  can  be  achieved  by  introduction  of  a  small  number  of  secondary 
sources  spaced  greater  than  half  a  wavelength  from  the  primary  source,  provided  that  the 
enclosure  is  being  excited  at  or  near  to  acoustic  resonance  (where  a  single  mode  dominates 
the  response).  In  the  low  frequency  case,  the  location  of  the  canceling  sources  was  found 
to  have  a  great  influence  on  the  performance  of  the  system.  It  was  found  that  at  excitation 
frequencies  where  only  a  single  mode  dominates  the  response  of  the  enclosure,  very  large 
reductions  in  potential  energy  can  be  achieved.  On  the  other  hand,  at  excitation  frequencies 
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where  many  modes  contribute  to  the  response  of  the  enclosure,  a  few  secondary  sources  are 
unable  to  control  these  modes  without  increasing  the  excitation  of  a  number  of  other  modes; 
hence,  little  reduction  in  the  total  acoustic  potential  energy  is  achieved. 

David  [15]  used  a  secondary  point  source  to  drive  the  acoustic  pressure  to  zero  at  a 
single  point  in  a  rectangular  enclosure  and  studied  the  result  of  this  action  on  the  rest  of 
the  sound  field  in  the  enclosure.  He  found  that  a  “zone  of  quiet”  around  the  cancellation 
point,  within  which  the  sound  pressure  level  is  reduced  by  more  than  15dB,  has  a  diameter 
of  about  one  tenth  of  the  wavelength  of  the  excitation  frequency.  The  mean-square  pressure 
away  from  the  point  of  cancellation  was  found  to  increase  if  the  transfer  impedance  between 
the  secondary  source  and  the  cancellation  point  was  very  small  at  the  excitation  frequency. 
This  led  the  authors  to  conclude  that  the  secondary  sources  should  be  placed  close  to  the 
cancellation  point  to  avoid  increasing  the  mean-square  pressure  away  from  the  cancellation 
point. 

Similar  studies  have  been  performed  for  the  free  field  case  [37].  In  these  studies,  the 
primary  source  and  the  secondary  sources  are  assumed  to  be  in  a  free  field  environment 
with  the  primary  source  radiating  at  a  single  known  frequency.  The  secondary  source 
strengths  that  minimize  the  power  output  of  the  combination  of  primary  and  secondary 
sources  is  then  found,  and  the  minimum  value  of  the  power  output  is  calculated.  It  is  found 
that  significant  reduction  in  the  combined  power  output  may  only  be  achieved  if  secondary 
sources  are  placed  within  a  distance  of  half  a  wavelength  of  the  primary  source. 

It  is  important  to  note  that  in  all  these  experiments  the  secondary  source  strengths 
were  calculated  based  on  the  exact  knowledge  of  the  strength  of  the  primary  source  for 
all  time.  Therefore,  these  results  can  only  be  used  as  rough  guidelines  for  predicting  the 
global  behavior  of  practical  ANC  systems  in  which  the  source  strength  is  unknown.  In  a 
practical  ANC  system,  the  source  strength  must  be  estimated  causally  based  on  the  available 
measurements. 

A  typical  ANC  system  utilizes  several  microphones  to  monitor  the  attenuated  field 
and/or  the  primary  field.  The  outputs  of  these  microphones  are  used  ais  inputs  to  some 
sort  of  electronic  controller.  The  controller  is  designed  to  drive  the  secondary  sources  in 
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such  a  way  that  the  desired  acoustic  objective  is  achieved.  The  desired  acoustic  objective 
of  almost  all  existing  ANC  systems  is  to  attenuate  the  sound  pressure  at  a  finite  number  of 
spatial  points.  In  the  next  Section,  we  will  concentrate  on  the  algorithms  that  are  used  by 
the  controller  to  achieve  this  specific  acoustic  objective. 


2.3  Control  Strategies  for  ANC 

A  recent  bibliography  of  references  for  ANC  contains  over  3450  entries  [27];  therefore,  it  is 
not  possible  to  cover  each  reference  individually.  However,  most  of  these  references  describe 
particular  applications  of  ANC  rather  than  new  concepts.  In  this  Section  we  concentrate 
on  those  control  strategies  that  are  most  commonly  used  in  active  noise  cancellation,  and 
selected  references  are  only  cited  to  illustrate  the  use  of  these  control  strategies. 

Most  existing  active  noise  cancellation  controllers  can  be  categorized  as  pointwise  feed¬ 
back  controllers  or  as  pointwise  feedforward  controllers.  The  specific  pointwise  feedback 
controller  and  pointwise  feedforward  controller  that  we  focus  on  in  this  thesis  are  presented 
next. 


2.3.1  Feedback  Controllers  for  Active  Noise  Cancellation 

A  generic  feedback  controller  of  the  type  we  focus  on  in  this  thesis  is  depicted  in  Fig.  2.1, 
where  the  input  to  the  controller  e{t)  is  the  sum  of  the  plant  output  c(t)  and  the  stochastic 
disturbance  n{t).  The  goal  is  to  choose  the  control  signal  r{t),  based  on  observations  of 
{e(r) :  r  <  t},  so  that  the  residual  signal  at  the  plant  output  e{t)  follows  a  desired  trajectory. 
Note  that  n{t),  e{t)  or  r{t)  can  be  scalar  valued  or  vector  valued.  The  important  assumption 
here  is  that  the  control  signal  is  generated  based  on  the  measurements  of  the  residual  signal 
at  the  plant  output,  and  not  based  on  direct  measurements  of  the  disturbance.  In  this  case, 
all  the  inputs  to  the  controller  contain  a  part  due  to  the  plant  output. 

An  illustrative  example  of  a  feedback  ANC  system  is  the  noise  canceling  headphone 
developed  by  Wheeler  [51].  Wheeler  used  an  analog  feedback  system  to  reduce  the  sound 
pressure  fluctuations  in  a  headphone  close  to  the  ear  of  the  listener.  The  noise  cancellation 
environment  and  the  corresponding  block  diagram  are  depicted  in  Fig.  2.2  and  Fig.  2.3, 
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n(t) 


Figure  2.1:  \3eneric  feedback  controller. 


respectively. 

Throughout  this  chapter  we  rely  on  the  following  model  of  the  microphone.  The  output 
of  the  microphone  is  assumed  to  be  the  value  of  the  acoustic  pressure  field  at  the  location 
of  the  microphone,  i.e.  the  microphone  is  assumed  to  measure  the  acoustic  pressure  at  the 
location  of  the  microphone  without  any  distortion  [6]. 

In  Fig.  2.3,  plant  G{s)  represents  the  overall  transfer  function  from  the  canceling  loud¬ 
speaker  input  r{t)  to  the  value  of  the  canceling  field  at  the  location  of  the  microphone 
c{t)  and  incorporates  the  transfer  functions  of  the  loudspeaJcer  and  the  propagation  path 
between  the  loudspeaker  and  the  microphone.  The  microphone  output  e{t)  is  the  sum  of 
the  unwanted  noise  at  the  location  of  the  microphone  n{t)  and  the  canceling  signal  at  the 
location  of  the  microphone  c(t).  The  goal  of  the  controller  is  to  minimize  the  sound  pres¬ 
sure  fluctuations  as  measured  by  the  microphone,  i.e.  the  controller  is  designed  to  minimize 
Comparing  Fig.  2.3  with  Fig.  2.1,  we  see  that  the  noise  canceling  headphone 
developed  by  Wheeler  is  a  scalar  analog  feedback  controller. 

Note  that  this  system  attenuates  the  noise  power  at  the  microphone,  and  it  is  assumed 
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that  attenuating  the  noise  at  the  microphone  will  result  in  noise  attenuation  at  the  ear  of 
the  listener.  In  noise  canceling  headphones  with  low  frequency  noise,  this  assumption  is 
quite  reasonable  because  the  ear  of  the  listener  is  only  about  3  centimeters  away  form  the 
microphone  [15]. 


Figure  2.2:  Noise  canceling  headphone  proposed  by  Wheeler. 


In  the  system  of  Fig.  2.3,  the  noise  field  at  the  microphone  without  the  noise  cancellation 
system  is  n(t).  With  the  noise  cancellation  system,  the  noise  field  at  the  microphone  is  e(t). 
The  power  spectra  of  these  two  signals  are  related  by  the  closed-loop  transfer  function  T{s) 
that  is  defined  as 


T{s) 


1 

l-H{s)G{s)  ’ 


where  H (s)  is  the  transfer  function  of  the  controller.  Specifically, 


(2.17) 


Peeim  =  |r(jfi)p  Pnnim.  (2.18) 

where  Pee{e^^)  is  the  power  spectrum  of  e(t),  and  PnnU^)  is  the  power  spectrum  of  n{t) 
with  fl  being  the  continuous- time  frequency  variable.  Looking  at  eq.  (2.18),  we  see  that  to 
achieve  attenuation  over  a  certain  frequency  band,  |r(jfl)|  must  be  made  small  over  this 
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r{t) 


J  Controller 

1  His) 

Figure  2.3:  Block  diagram  of  the  noise  canceling  headphone  proposed  by  Wheeler. 

band.  Typically,  the  analog  controller  is  designed  so  that  the  closed-loop  system  achieves 
moderate  noise  attenuation  (15dB)  over  a  relatively  wide  range  of  frequencies  (50Hz-500Hz) 
[12,  51].  Note  that  due  to  its  analog  nature,  the  controller  in  Fig.  2.2  cannot  be  changed 
once  it  is  built. 

A  Discrete-Time  Feedback  Controller 

A  discrete-time  version  of  the  system  in  Fig.  2.2  has  been  proposed  by  Graupe  [26].  Since 
the  controller  proposed  by  Graupe  is  a  discrete-time  controller  we  assume  throughout  this 
Section  that  the  noise  to  be  canceled,  the  canceling  signal,  and  all  other  intermediate  signals 
are  appropriately  sampled.  Consequently  in  this  Section,  “t”  represents  the  normalized 
sampling  time.  Note  that  up  to  this  point,  we  have  been  using  “f’  as  the  continuous- time 
variable.  Occasionally  this  notation  can  be  ambiguous;  however,  we  think  that  in  most 
cases  the  reader  will  be  able  to  resolve  this  ambiguity  based  on  the  context.  Whenever  this 
ambiguity  cannot  be  resolved  based  on  the  context,  we  explicitly  state  whether  “t”  denotes 
the  continuous-time  variable  or  the  discrete-time  variable. 

The  block  diagram  for  the  system  of  Graupe  is  depicted  in  Fig.  2.4,  where  G{z)  is 
the  overall  transfer  function  from  the  canceling  loudspeaker  input  r{i)  to  the  value  of  the 
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canceling  field  at  the  location  of  the  microphone  c{t)  and  incorporates  the  transfer  function 
of  the  loudspeaker  and  the  propagation  path  between  the  loudspeaker  and  the  microphone. 
The  microphone  output  e(t)  is  the  sum  of  the  unwanted  noise  n{t)  and  the  canceling  signal 
c{t).  n{t)  is  modeled  as  the  output  of  the  transfer  function  $(z)  driven  by  the  white  noise 
process  w{t).  The  goal  of  the  controller  is  to  minimize  the  sound  pressure  at  the  microphone, 
i.e.  the  controller  is  designed  to  minimize  E{e'^{t)}. 


Figure  2.4:  Block  diagram  of  the  ANC  system  proposed  by  Graupe. 

A  modified  version  of  the  Box- Jenkins  controller  [8]  is  used  to  continuously  adjust  ff(z) 
such  that  E{e{t)'^)  is  minimized.  This  system  uses  the  noise  model  $(<;)  to  adapt  the 
controller  to  the  specific  type  of  noise  that  is  being  canceled.  Based  on  computer  simulations, 
the  authors  claim  that  this  system  achieves  54dB  of  noise  cancellation  in  the  context  of 
compressor  noise. 

Incorporation  of  a  stochastic  model  for  the  unwanted  noise  in  the  system  proposed  by 
Graupe  is  one  reason  for  the  much  better  performance  of  this  system  than  the  performance 
of  the  system  proposed  by  Wheeler.  The  digital  system  of  Graupe  uses  the  model  of  the 
unwanted  noise  to  determine  those  frequencies  at  which  the  noise  signal  has  significant 
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energy,  and  it  adjusts  H{z)  such  that  the  resulting  closed-loop  transfer  function,  = 

small  over  these  frequencies.  Since  the  compressor  noise  is  extremely 
narrow  band,  the  resulting  closed-loop  system  essentially  acts  as  a  notch  filter.  On  the  other 
hand,  the  analog  system  proposed  by  Wheeler  is  fixed  and  makes  no  use  of  the  spectral 
characteristics  of  the  unwanted  noise;  therefore,  it  achieves  modest  noise  attenuation  over 
a  wide  range  of  frequencies. 

The  second  reason  for  the  excellent  results  obtained  by  Graupe  is  that  his  simulations 
were  performed  with  a  G{z)  that  did  not  contain  any  delay.  Moreover,  the  conroller  proposed 
by  Graupe  is  derived  based  on  the  explicit  assumption  that  G{z)  contains  no  delay.  This 
severely  limits  the  applications  in  which  this  algorithm  can  be  used,  since  there  are  many 
scenarios  in  which  the  propagation  delay  between  the  loudspeaker  and  the  error  microphone 
is  of  the  order  of  few  sample  intervals. 

2.3.2  Feedforward  Controllers  for  Active  Noise  Cancellation 

A  generic  feedforward  controller  of  the  type  we  focus  on  in  this  thesis  is  depicted  in  Fig. 
2.5,  where  q{t)  is  the  input  of  the  controller,  and  n{i)  is  the  stochastic  disturbance.  The 
goal  is  to  choose  the  controller  such  that  the  residual  signal  at  the  plant  output  e{t)  is  as 
close  as  possible  to  a  desired  trajectory.  In  this  case,  the  input  of  the  controller  contains 
no  part  due  to  the  plant  output  c(t),  i.e.  there  is  no  feedback  from  the  plant  output  to  the 
inputs  of  the  controller.  This  is  in  contrast  to  the  feedback  scheme  of  the  previous  Section 
where  the  input  of  the  controller  was  the  sum  of  the  plant  output  and  the  disturbance. 


n(t) 


Figure  2.5:  Generic  feedforward  controller. 
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Single-Channel  Feedforward  ANC  Systems 

Most  of  the  early  work  in  ANC  was  related  to  attenuation  of  duct  noise  at  low  frequencies 
using  feedforward  controllers.  A  very  thorough  analysis  of  the  duct  problem  was  carried 
out  by  Swinbanks  in  1973  [49].  We  shall  present  his  work  in  some  detail  next  because  it 
illustrates  many  of  the  issues  involved  in  using  feedforward  controllers  for  ANC.  The  system 
developed  by  Swinbanks  is  a  continuous-time  system;  therefore,  all  the  signals  involved  in 
this  Section  are  continuous-time  signals. 

Figure  2.6  is  a  simplified  diagram  of  the  system  proposed  by  Swinbanks.  The  primary 
wave,  propagating  in  the  down  stream  direction  of  an  infinitely  long  duct  with  square 
cross  Section,  is  detected  by  a  microphone  which  supplies  the  information  necessary  for  the 
operation  of  the  controller.  The  controller  generates  a  secondary  field  in  zone  3  which  is  the 
negative  of  the  primary  field  in  this  zone.  For  the  microphone  to  only  detect  the  primary 
field,  the  secondary  field  must  not  propagate  upstream  towards  zone  1.  Moreover,  if  there 
were  any  upstream  obstructions  in  zone  1,  the  secondary  field  propagating  upstream  would 
get  reflected  back  to  zone  3,  and  one  would  have  to  deal  with  the  problem  of  controlling 
both  the  primary  wave  and  the  reflected  wave  in  zone  3. 


Zone  2 


Figure  2.6:  Sketch  of  the  ANC  system  proposed  by  Swinbanks. 


Swinbanks  showed  that  by  using  a  single  ring  sources,  consisting  of  four  point  sources,  it 
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is  possible  to  generate  an  output  consisting  only  of  a  propagating  plane  wave  for  frequencies 
up  to  2.1  times  the  cut-oflf  frequency  of  the  duct,  see  Fig.  2.7.  More  precisely,  for  these 
frequencies,  the  non-plane  wave  modes  in  the  duct  will  decay  exponentially  along  the  long 
axis  of  the  duct  away  from  the  ring  source.  Swinbanks  also  showed  that  by  combining  the 
effects  of  two  such  ring  sources,  it  is  possible  to  generate  an  approximately  unidirectional 
plane  wave  in  the  duct  up  to  2.1  times  the  cut-off  frequency  of  the  duct,  see  Fig.  2.7. 
Specifically,  if  the  source  strength  of  the  upstream  ring  is  mi(t)  and  the  strength  of  the 
downstream  ring  is  m2{t),  there  will  be  approximately  zero  output  in  the  upstream  direction 
provided  that 

mi{t)  = -m2{t  -  To),  (2.19) 

where  tq  is  an  appropriate  amount  of  delay  which  depends  on  the  separation  between  the 
two  rings  and  the  speed  of  sound  in  the  duct. 

As  illustrated  in  Fig.  2.7,  the  controller  in  this  case  consists  of  a  delay  element 
followed  by  an  analog  filter  whose  purpose  is  to  equalize  the  combined  frequency  response 
of  the  secondary  sources.  Let  p{x,t)  denote  the  pressure  field  at  location  x  (measured  along 
the  long  ajcis  of  the  duct)  and  at  time  t.  Over  the  frequency  range  of  interest,  the  equalizer 
tries  to  make  p(b,t)  —  m2{t  —  T2)  for  some  delay  T2- 

The  operation  of  the  ANC  system  of  Swinbanks  can  be  summarized  as  follows.  At 
time  t,  the  controller  predicts  p{b,  t  -H  T2)  and  drives  the  equalized  secondary  sources  with 
the  negative  of  this  predicted  value.  If  the  prediction  is  perfect,  the  primary  field  and  the 
secondary  field  will  interfere  destructively  at  a;  =  6,  and  the  result  will  be  complete  silence 
at  this  point.  However,  every  other  point  in  zone  3  wiU  also  be  completely  silent,  since  the 
residual  field  in  zone  3  is  a  plane  wave  propagating  to  the  right. 

Implementation  of  Swinbanks  method,  with  minor  modifications,  resulted  in  a  flurry  of 
publications  (e.g.,  [7,  11,  41,  44]).  These  implementations  showed  that  the  main  drawback 
of  Swinbanks’  method  is  that  very  precise  analog  electronics  are  needed  to  implement  the 
feedforward  controller.  The  required  degree  of  precision  in  the  magnitude  and  the  phase  of 
the  controller  was  found  to  be  impossible  to  maintain  except  over  a  very  limited  range  of 
frequencies.  In  short,  Swinbanks’  method  was  found  to  lack  robustness. 
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Figure  2.7:  Block  diagram  of  the  ANC  system  proposed  by  Swinbanks. 

The  most  interesting  feature  of  the  system  proposed  by  Swinbanks  is  that  it  achieves 
complete  noise  cancellation  at  every  point  in  the  downstream  Section  of  the  duct,  although 
the  controller  is  designed  to  cancel  the  noise  at  a  single  point.  This  is  a  direct  consequence 
of  the  unidirectionality  of  the  primary  and  secondary  fields.  In  most  applications  of  ANC, 
the  sound  fields  involved  are  diffused;  hence,  the  method  of  Swinbanks  often  cannot  be 
used. 

In  1981  Burgess  proposed  an  adaptive  discrete-time  feedforward  controller  for  ANC  [10]. 
The  system  proposed  by  Burgess  is  quite  robust  and  does  not  rely  on  the  unidirectionality  of 
the  sound  fields  involved.  Therefore,  this  system  can  be  used  in  a  wide  range  of  applications. 
The  acoustic  goal  of  this  system  is  to  attenuate  the  noise  at  the  location  of  a  microphone 
based  on  the  information  provided  by  a  reference  signal. 

As  depicted  in  Fig.  2.8,  the  system  proposed  by  Burgess  consists  of  an  error  microphone 
and  a  canceling  loudspeaker.  The  goal  is  to  to  attenuate  the  noise  at  the  location  of  the  error 
microphone.  The  reference  signal  q{t)  provides  information  about  the  unwanted  noise  at  the 
error  microphone  and  is  the  input  of  the  controller  H.  If  the  noise  is  generated  by  a  pump 
or  an  engine,  the  reference  signal  is  typically  the  output  of  a  non-acoustical  sensor  attached 


33 


to  the  pump  or  the  engine  (e.g.,  a  tachometer  on  the  engine  fly  wheel).  The  controller  is 
a  discrete-time  finite  impulse  response  system  whose  coefficients  are  continuously  adjusted 
to  minimize  the  power  of  the  residual  signal  e(t)  at  the  error  microphone. 


Primary  Field 


Error  Microphone 


Reference  Signal 


H 

(FIR  System) 


Feedforward  Cmtroller 


Cancelling  Loudspeaker 


Figure  2.8:  ANC  system  proposed  by  Burgess. 


The  block  diagram  for  the  system  in  Fig.  2.8  is  depicted  in  Fig.  2.9,  where  q{i)  is 
the  reference  signal;  H  is  the  FIR  controller  whose  coefficients  at  time  “t”  are 
and  G{z)  is  the  transfer  function  from  the  input  of  the  loudspeaker  r{t)  to  the  value  of 
the  canceling  field  at  the  location  of  the  error  microphone  c{t).  The  output  of  the  error 
microphone  e{t)  is  the  sum  of  the  unwanted  noise  at  the  location  of  the  microphone  n{t)  and 
the  value  of  the  canceling  field  at  the  location  of  the  microphone  c{t).  This  block  diagram 
is  derived  under  the  assumption  that  there  is  no  feedback  from  the  output  of  the  canceling 
loudspeaker  to  the  reference  signal. 


Feedforward  Controller 


Figure  2.9:  Block  diagram  of  the  ANC  system  proposed  by  Burgess. 

Assuming  that  G{z)  is  known,  an  LMS-type  algorithm  [52]  can  be  derived  for  adjusting 
the  coefficients  such  that  the  mean-square  power  of  the  residual  signal  e(<)  is 
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minimized.  The  update  equations  for  the  coefficients  of  the  controller  are: 

h{i,t  +  1)  =  h{i,t) -2  IJ,  e{t)  v{t  -  i)  ,  (2.20) 

where  u(t)  is  defined  as  v{t)  q{t)  * g{t),  and  fi  is  the  step-size  of  the  algorithm.  Note  that 
v{t)  can  be  computed  by  “filtering”  the  known  signal  q{t)  with  the  known  impulse  response 
g{t)]  hence,  this  algorithm  is  sometimes  referred  to  as  the  “filtered- x”  algorithm  [53]. 

Multiple-Channel  Feedforward  ANC  Systems 

A  multiple- channel  version  of  the  algorithm  of  Burgess  has  been  developed  by  Elliot  et. 

al  [21,  18].  This  algorithm  is  a  discrete-time  algorithm;  therefore,  all  signals  involved  are 

assumed  to  be  appropriately  sampled.  As  depicted  in  Fig.  2.10,  this  system  consists  of  L 

microphones  and  M  loudspeakers.  The  controller  in  this  case  is  a  bank  of  M  finite  impulse 

response  filters  where  each  filter  is  driven  by  the  reference  signal  q{t).  The  goal 

is  to  adjust  the  coefficients  of  these  filters  such  that  a  cost  function  involving  the  mean- 

square  sum  of  the  microphone  outputs  and  the  mean-square  sum  of  the  control  signals  is 

minimized.  Specifically,  given  the  matrix  of  transfer  functions  from  the  loudspeaker  inputs 

to  the  microphone  outputs,  the  authors  develop  an  LMS-type  algorithm  for  adjusting  the 

coefficients  of  these  M  filters  such  that 

{L  M  ] 

^e?(0  +  a^r2(t)  | 

is  minimized,  where  a  is  a  positive  weighting  coefficient  for  the  control  effort.  Note  that 
the  cost  function  used  here  is  a  bit  more  general  than  the  one  used  by  Burgess  [10]. 

The  influence  of  the  loudspeaker  transfer  function  and  acoustic  delay  on  the  performance 
of  this  algorithm  has  been  thoroughly  analyzed  in  a  number  of  papers  (e.g.,  [20,  47,  48]). 

2.4  Performance  Analysis  for  Two  ANC  Controllers 

The  performance  of  the  continuous- time  single-channel  feedback  ANC  system  and  the  per¬ 
formance  of  the  continuous- time  single-channel  feedforward  ANC  system  are  analyzed  in 
this  Section. 
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Figure  2.10:  Multiple-channel  feedforward  ANC  system. 

2.4.1  Optimal  Single-Channel  Feedback  Controller 

Further  insight  about  the  continuous-time  feedback  ANC  system  of  Wheeler  [51]  can  be 
gained  by  finding  the  optimal  feedback  controller  in  this  case.  Specifically  in  Fig.  2.3,  we 
would  like  to  find  the  causal  controller  H{s)  that  minimizes  assuming  that  the 

noise  n{t)  is  stationary  with  known  power  spectrum  Pnn{j^)  and  assuming  that  the  plant 
G(s)  is  known.  Our  strategy  is  to  first  find  the  minimizing  controller  within  a  restricted 
class  of  controllers,  and  then  show  that  no  other  controller  can  outperform  this  controller. 
We  shall  make  the  further  assumption  that  G(s)  is  stable.  This  is  a  reasonable  assumption, 
since  the  physical  system  that  G(s)  corresponds  to  is  passive. 

Let  us  start  by  looking  at  controllers  of  the  form  depicted  in  Fig.  2.11.  In  this  figure, 
G{s)  (the  transfer  function  of  the  plant)  is  fixed,  and  Q{s)  is  a  causal  and  stable  system  that 
can  be  varied  subject  to  the  restriction  that  the  resulting  closed-loop  system  must  be  stable. 
Let  us  take  an  arbitrary  controller  within  this  class  and  assume  that  this  controller  results 
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from  some  <3(s),  then  the  closed-loop  transfer  function  corresponding  to  this  controller  is 


1 

l-Gis)H{s) 
_ 1 _ 

=  1  -l-  G{s)Q{s). 


(2.21) 

(2.22) 

(2.23) 


It  is  important  to  note  that  since  both  Q(5)  and  G(s)  are  stable,  there  is  no  unstable 
pole/zero  cancellation  in  the  closed-loop  system. 

With  this  controller,  the  power  spectrum  of  the  residual  signal  e{t)  can  be  expressed  in 
terms  of  the  power  spectrum  of  the  unwanted  noise  n(t)  as 

Pee(ifi)  =  |1  +  Gijn)  Qijnf  Pnnim-  (2-24) 


Therefore,  the  first  step  is  to  find  the  causal  and  stable  Q(6)  in  eq.  (2.24)  that  minimizes 
the  power  of  the  residual  signal  e{t)  for  a  given  Pnn(s). 


H(s) 


Figure  2.11:  Restricted  cIclss  of  controllers. 


Figure  2.12:  Equivalent  block  diagram  of  the  the  closed-loop  system. 

A  convenient  way  to  find  the  minimizing  choice  for  Q{s)  is  to  represent  equation  (2.23)  in 
block  diagram  form  as  depicted  in  Fig.  2.12.  Looking  at  Fig.  2.12,  it  is  clear  that  the  choice 
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for  Q(5)  that  minimizes  E{€^{t)}  is  the  causal  Wiener  filter  [11]  that  produces  the  linear 
least-squares  estimate  (ELSE)  of  —n(t)  based  on  {p(r)  :  t  <t}.  Let  us  denote  this  particular 
choice  for  Q{s)  by  Qwin(s),  the  corresponding  controller  by  ,  and 

the  resulting  E{e^{t)}  by  Eujin-  From  eq.  (2.23),  we  see  that  the  closed-loop  system  in  this 
case  is  also  stable,  since  G{s)  and  Qwin{s)  are  both  stable. 

It  remains  to  be  shown  that  no  other  causal  controller  can  outperform  We  will 

prove  this  by  contradiction.  To  this  end,  assume  that  such  a  controller  IE^{s)  exists,  and 
the  resulting  value  of  E{e^{t)}  is  E*,  which  is  less  than  E^in-  It  can  be  easily  shown  that 
this  leads  to  a  contradiction  by  choosing 


Qis)  = 


H*{s) 

1  -  H*{s)G{s) 


(2.25) 


in  Fig.  2.12  and  noting  that  the  resulting  value  for  E{e^(t)}  is  E*,  which  is  assumed  to 
be  less  than  E^in.  This  is  a  contradiction,  since  E^in  is  obtained  by  minimizing  E{e'^{t)} 
over  all  possible  causal  choices  for  Q{s).  Note  that  Q(s)  in  eq.  (2.25)  is  causal,  since  it  can 
be  realized  as  feedback  interconnection  of  jff^(s)  and  C?(s).  Using  the  uniqueness  of  the 
linear  least-squares  estimate,  it  can  be  shown  that  the  canceling  signal  c(t)  that  minimizes 
E{€^(t)}  is  unique. 

This  analysis  shows  that  with  the  optimal  controller,  the  canceling  signal  c(i)  equals  the 
ELSE  of  —n(t)  based  on  {p(r)  :  r  <  t),  where  p(t)  =  n{t)*g{t).  This  analysis  can  be  taken 
one  step  further.  To  this  end,  let  us  express  G{s)  as 


G(s)  —  Gall{s)Gmin{^), 


(2.26) 


where  Gaii{s)  is  an  all-pass  system  function,  and  Gmin(s)  is  a  minimum-phase  system  func¬ 
tion.  We  also  define  p'{t)  as 

p'it)  =  n{t)*gaii{t).  (2.27) 

The  ELSE  of  n{t)  based  on  {p(r)  :  r  <  t}  is  the  same  as  the  ELSE  of  n(t)  based  on 
{p'{t)  :r  <t},  since  p{t)  =  p'(t)  *  gmin{i)  and  Gmin{s)  is  minimum-phase.  Recalling  that 
the  minimizing  choice  for  c(t)  is  the  ELSE  of  —n{t)  based  on  {p(r) :  t  <  t},  we  see  that  the 
minimizing  choice  for  c{t)  can  be  alternatively  expressed  as  the  ELSE  of  —n{t)  based  on 
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{p'(r)  :  T  <  t].  Hence,  the  optimal  performance  of  the  noise  cancellation  system  depends 
on  the  joint  second-order  statistics  of  {n(t),p'(t)},  which  can  be  calculate  form  the  power 
spectrum  of  n{t),  since  p'{t)  =  n{t)  *  gaiiit)  and  gaii{t)  is  known.  Clearly,  noise  attenuation 
will  be  high  whenever  p'{t)  is  highly  correlated  with  n{t). 

If  the  phase  of  Ga:i(s)  is  approximately  linear  over  the  frequencies  for  which  n(i)  has 
significant  energy,  the  problem  will  reduce  to  predicting  -n(t)  based  on  {n(r)  :  r  <  t-  tq), 
where  tq  is  the  slope  of  the  linear  approximation  to  the  phase  of  Gaii(s). 

2.4.2  Optimal  Single-Channel  Feedforward  Controller 

Further  insight  about  the  feedforward  controller  of  Burgess  can  be  gained  by  looking  at  the 
system  in  continuous-time.  The  block  diagram  for  the  system  in  continuous-time  is  depicted 
in  Fig.  2.13,  where  g(t)  is  the  reference  signal;  ff(s)  is  a  continuous-time  LTI  controller;  and 
G(s)  is  the  transfer  function  from  the  loudspeaker  input  r(i)  to  the  value  of  the  canceling 
field  at  the  location  of  the  error  microphone.  The  output  of  the  error  microphone  e(f)  is 
the  sum  of  the  unwanted  noise  n(i)  and  the  canceling  signal  c(t). 


n(t) 

q(t)  j  Controller  ■  r(t) 

G(s) 

1 _ 

"1  h(S)  ; 

Figure  2.13:  Block  diagram  of  the  feedforward  ANC  system  in  continuous  time. 

Let  us  express  G(s)  as 

G(s)  =  Galt(s)Gmm(^)f  (2.28) 

where  Gaii{s)  is  an  all-pass  system  function  and  Gmin(s)  is  a  minimum-phase  system  func¬ 
tion.  Assume  that  n(t)  and  q(i)  are  jointly  stationary.  In  Fig.  2.13,  G(s)  and  if(s)  can  be 
interchanged  without  affecting  c(t)  or  e(t).  The  resulting  system  is  depicted  in  Fig  2.14. 

In  Fig.  2.14,  the  choice  for  If(s)  that  minimizes  F{e^(t)}  is  the  causal  Wiener  filter 
[57]  that  produces  the  linear  least-squares  estimate  of  —n(t)  based  on  {p(r)  :  t  <  t).  The 
minimum  value  of  E{e^(i)}  can  be  calculated  based  on  the  joint  second-order  statistics  of 
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{7i(f),p(i)}  ,  which  in  turn  can  be  calculated  based  on  the  joint  second-order  statistics  of 
{Ti{t),q{t)},  since  p{t)  =  q{t)  *  g{t)  and  g{t)  is  known. 


G(s) 


Figure  2.14:  Alternate  block  diagram  for  the  ANC  system  proposed  by  Burgess. 

This  analysis  can  be  taken  one  step  further  by  observing  that  the  linear  least-squares 
estimate  of  n(t)  based  on  {p(r)  :  r  <  t}  is  equivalent  to  the  linear  least-squares  estimate  of 
n{t)  based  on  {p'\r)  :  r  <  t},  since  p{t)  =  p"{t)  *  gm{n{s)  and  Gmin{s)  is  minimum  phase. 
Therefore,  the  noise  attenuation  will  be  high  whenever  p"{t)  is  highly  correlated  with  n{i). 

If  the  phase  of  Gaii{s)  is  approximately  linear  over  those  frequencies  for  which  q{t)  has 
significant  energy,  the  problem  will  reduce  to  predicting  —n{t)  based  on  {q(r)  :  r  <  t  -  ro}, 
where  tq  is  the  slope  of  the  linear  approximation  to  the  phase  of  Gaii{s).  In  this  case,  the 
noise  attenuation  is  substantial  whenever  q{t)  is  highly  correlated  with  the  future  values  of 
n{t).  Note  that  this  approximation  is  quite  reasonable  whenever  the  noise  to  be  canceled  is 
relatively  narrow  band. 

A  major  shortcoming  of  the  ANC  system  developed  by  Burgess  is  discussed  next.  The 
above  analysis  implies  that  the  system  of  Burgess  works  well  whenever  the  reference  signal 
is  highly  correlated  with  the  noise  at  the  error  microphone,  and  there  is  no  feedback  from 
the  canceling  loudspeaker  to  the  reference  signal.  However,  the  formulation  of  Burgess 
does  not  tell  us  how  to  obtain  such  a  reference  signal.  The  author  suggests  the  use  of 
non- acoustical  sensors  attached  to  the  noise  source.  In  many  applications,  the  source  of 
the  noise  is  inaccessible  or  unknown;  hence,  it  is  often  not  possible  to  acquire  the  reference 
signal  in  this  way. 

A  common  practice  advocated  by  many  researchers  is  to  use  the  output  of  a  second 
microphone  (called  the  reference  microphone)  as  the  reference  signal.  Unfortunately,  there 
is  typically  some  feedback  from  the  canceling  loudspeaker  to  the  output  of  the  reference 
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microphone.  In  the  presence  of  this  feedback,  the  algorithm  proposed  by  Burgess  no  longer 
minimizes  £{e^(t)};  moreover,  the  overall  system  has  been  observed  to  go  unstable  in  the 
presence  of  this  feedback.  To  decrease  the  feedback  from  the  canceling  loudspeaker  to  the 
reference  microphone,  the  reference  microphone  is  usually  placed  far  away  form  the  canceling 
loudspeaker  and  the  error  microphone.  However,  this  decreases  the  correlation  between  the 
noise  at  the  error  microphone  and  the  noise  at  the  reference  microphone,  resulting  in  low 
levels  of  noise  attenuation. 

The  multiple-channel  feedforward  ANC  algorithm  proposed  by  Elliot  et.  al  [21,  18]  suf¬ 
fers  from  the  same  shortcoming  that  the  single- channel  algorithm  of  Burgess  does.  Specifi¬ 
cally,  there  is  the  problem  of  how  to  acquire  a  reference  signal  that  is  highly  correlated  with 
the  unwanted  noise  and  at  the  same  time  avoid  acoustic  feedback  from  the  output  of  the 
canceling  loudspeakers  to  the  reference  signal. 

To  circumvent  the  acoustic  feedback  problem  associated  with  the  use  of  a  reference 
microphone,  Eriksson  et.  al.  [23]  proposed  an  ANC  system  that  takes  the  effect  of  the  feed¬ 
back  into  account  in  the  design  of  the  adaptive  controller.  Specifically,  Eriksson  proposed 
canceling  the  effect  of  the  feedback  at  the  output  of  the  reference  microphone,  i.e.  the  input 
of  the  controller  in  this  case  is  the  sum  of  the  output  of  the  reference  microphone  and  the 
negative  of  the  estimate  of  the  feedback.  The  resulting  system  uses  two  LMS  algorithms 
and  is  quite  complicated.  More  importantly,  instability  can  result  if  the  parameters  of  the 
adaptive  algorithm  are  arbitrarily  set,  and  the  authors  offer  no  procedure  for  properly  set¬ 
ting  these  parameters.  Furthermore,  this  system  still  works  poorly  whenever  the  noise  at 
the  reference  microphone  is  not  highly  correlated  with  the  noise  at  the  error  microphone. 

2.5  Relationship  between  Signal  Estimation  and  Active  Noise 
Cancellation 

The  relationship  between  signal  estimation  and  active  noise  cancellation  is  explored  in  this 
Section. 

Signal  estimation  refers  to  estimating  one  stochastic  process  from  observations  of  another 
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stochastic  process.  The  goal  in  signal  estimation  is  to  generate,  based  on  an  observed  signal, 
an  output  signal  that  is  as  close  as  possible  to  some  desired  signal.  In  a  large  class  of 
signal  estimation  problems,  the  desired  signal  and  the  output  signal  are  both  numerical 
representation  of  some  underlying  physical  signal.  For  example,  in  speech  enhancement  the 
desired  signal  is  a  numerical  representation  of  the  pressure  fluctuations  produced  by  the 
vocal  cords  of  the  speaker,  and  the  output  of  the  speech  enhancement  system  is  a  numerical 
estimate  of  these  pressure  fluctuations.  In  this  class  of  problems,  all  the  signals  involved 
(including  the  desired  signal)  are  numerical;  hence,  there  is  no  need  to  consider  the  physical 
signals  that  these  numerical  signals  represent.  Moreover  in  this  class  of  problems,  the 
process  by  which  the  output  signal  is  obtained  is  purely  numerical.  In  the  reminder  of  this 
Section  the  term  “signal  estimation”  refers  exclusively  to  this  type  of  numerical  estimation. 

In  active  noise  cancellation,  the  desired  signal  is  a  physical  signal  and  not  a  numerical 
signal;  therefore,  the  output  signal  in  active  noise  cancellation  is  also  a  physical  signal. 
Furthermore,  the  output  signal  in  ANC  is  generated  as  the  sum  of  a  canceling  signal  and 
one  of  the  inputs.  Since  the  output  signal  is  a  physical  signal,  the  canceling  signal  and  the 
input  signal  to  which  the  canceling  signal  is  added  are  both  physical  signals  as  well.  Hence, 
every  active  noise  cancellation  system  must  have  a  transducer  for  physically  generating  the 
canceling  signal.  If  the  input  of  this  transducer  is  significantly  different  from  its  output,  the 
design  of  the  ANC  system  must  take  the  effect  of  the  transducer  into  account. 

For  example,  in  acoustic  active  noise  cancellation,  the  canceling  signal  and  the  unwanted 
noise  are  both  physical  signals  corresponding  to  the  primary  pressure  field  and  the  secondary 
pressure  field,  respectively.  The  canceling  signal  is  generated  by  a  loudspeaker  and  is  added 
physically  to  the  unwanted  noise  through  the  interference  of  the  primary  field  and  the 
secondary  field.  Note  that  the  process  by  which  the  output  signal  is  obtained  is  not  purely 
numerical,  since  the  output  signal  is  obtained  through  physical  interference  of  two  acoustic 
fields. 

The  relationship  between  signal  estimation,  in  which  the  desired  signal  is  a  numerical 
signal,  and  active  noise  cancellation,  in  which  the  desired  signal  is  a  physical  signal,  is 
explored  in  this  Section.  In  Section  2.5.1,  a  very  general  class  of  single- input/single-output 
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signal  estimation  systems  is  compared  to  the  class  of  single-channel  feedback  ANC  systems. 
In  Section  2.5.2,  the  class  of  single-channel  feedforward  ANC  systems  is  compared  to  a 
special  class  of  two-input /single-output  signal  estimation  systems. 


2.5.1  Relationship  between  Single-Channel  Signal  Estimation  and  Single- 
Channel  ANC 

The  relationship  between  single- channel  signal  estimation  and  single-channel  feedback  active 
noise  cancellation  is  studied  in  this  Section. 


H 

S 

(b) 

Figure  2.15:  General  single-channel  signal  estimation  system. 


Let  us  consider  a  completely  general  single-input/single-output  signal  estimation  system. 
The  estimator  can  be  causal,  non-causal,  linear,  or  non-linear.  As  depicted  in  Fig.  2.15(a), 


this  system  is  a  single-input /single-output  system  with  input  Xs{t)  and  output  es(t).  The 
objective  is  to  find  the  system  Hs  that  makes  es{t)  as  close  as  possible  to  some  desired 
signal  ds{t).  Alternatively,  the  single-channel  signal  estimation  system  can  be  expressed 
in  block  diagram  form  as  depicted  in  Fig.  2.15(b),  where  =  Hs  —  I  with  I  being  the 
identity  system.  Referring  to  Fig.  2.15(b),  the  single- channel  signal  estimation  problem  can 
be  stated  as  follows: 

Find  the  single-input/single-output  system  H's  that  minimizes 

J{{[xs{T)  +  H's(xs(T)),ds(T)]:  0<r<T}),  (2.29) 

where  J(.)  is  a  pre-specified  cost  function  that  indicates  how  close  the  output  es{t)  is  to  the 
desired  signal  ds{t)  over  the  time  interval  0  <  t  <  T,  e.g. 

J  ({ [x,(r)  -h  H'ixsir)),  ^.(r)]  :  0  <  r  <  T})  =  (2.30) 

T 

i;|[a;5(r)-l- /f'(xs(r))  -  <is(r)]^}dr.  (2.31) 

As  depicted  in  Fig.  2.16(a),  a  single-channel  feedback  ANC  system  is  a  single-input/single- 
output  feedback  system  with  input  Xa(t)  and  output  ea{t).  The  output  €a(<)  is  generated 
as  the  sum  of  the  canceling  signal  Ca(t)  and  the  input  signal  Xa{t).  Moreover,  the  canceling 
signal  Caft)  is  generated  as  the  output  of  system  G  with  input  r{t),  and  it  is  assumed  that 
G  is  not  an  identity  system.  In  this  system,  Xa{t)  and  Ca(t)  are  both  physical  signals,  and 
system  G  represents  the  effect  of  the  transducer  that  physically  generates  the  canceling 
signal. 

The  objective  of  the  single-channel  feedback  ANC  system  is  to  make  ea{t)  as  close  as 
possible  to  some  desired  signal  da{t).  Specifically,  the  single- channel  feedback  ANC  problem 
is  to  find  the  system  Ha  in  Fig.  2.16(a)  that  makes  ea{t)  as  close  as  possible  to  da(t). 
Following  the  same  steps  used  for  the  derivation  of  the  optimal  feedback  ANC  system 
in  Section  2.4.1,  it  can  be  shown  that  the  single- channel  feedback  ANC  problem  can  be 
equivalently  stated  as  finding  the  system  H'^  in  Fig.  2.16(b)  that  makes  ea(0  as  close  as 
possible  to  the  desired  signal  da(t).  Referring  to  Fig.  2.16(b),  the  single-channel  ANC 
problem  can  be  stated  as: 
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Find  the  single-input/single-output  system  that  minimizes 

J  ({ K{r)  +  {GoH',)Mt)),  4(t)]  :  0  <  r  <  T})  ,  (2.32) 

where  “o”  is  the  composition  operator. 

Comparing  the  single-channel  signal  estimation  problem  to  the  single-channel  feedback 
ANC  problem,  we  see  that  the  ANC  problem  can  be  viewed  as  a  constrained  version  of 
the  signal  estimation  problem.  Specifically,  the  single- channel  feedback  ANC  problem  can 
be  formulated  as  finding  that  minimizes  the  expression  in  eq.  (2.29)  subject  to  the 
constraint  that  must  be  of  the  following  form 

K  =  GoK  ,  (2.33) 

for  some  system  and  a  given  system  G. 

We  now  compare  the  single-channel  signal  estimation  system  to  the  single- channel  feed¬ 
back  ANC  system  in  terms  of  their  input/output  characteristics.  To  this  end,  we  assume 
that  both  systems  are  driven  by  the  same  input  and  that  the  desired  signal  for  both  systems 
is  the  same.  We  also  assume  that  each  system  is  designed  so  that  the  the  distance  between 
its  output  and  the  desired  signal  is  minimized. 

Under  these  assumptions,  the  performance  of  the  ANC  system  is  never  better  than  the 
performance  of  the  corresponding  signal  estimation  system,  since  the  optimal  ANC  system 
is  obtained  by  minimizing  J(.)  subject  to  a  constraint,  while  the  optimal  signal  estimation 
system  is  obtained  by  minimizing  J(.)  without  any  constraints. 

Again  under  these  assumptions,  if  the  minimizing  in  eq.  (2.29)  happens  to  equal 
{GoQ)  for  some  system  Q,  the  optimal  ANC  system  will  be  identical  to  the  optimal  signal 
estimation  system.  Specifically,  in  this  case,  the  minimizing  if'  in  the  ANC  problem  is 
system  Q.  For  example,  if  G  is  invertible,  the  minimizing  Hg  can  be  expressed  as 

H’  =  loH'g  (2.34) 

=  Go{G-^oH'g) ,  (2.35) 

where  I  is  the  identity  system.  This  implies  that  the  solution  to  the  single-channel  feedback 
ANC  problem  can  be  expressed  in  terms  of  the  solution  to  the  corresponding  single- channel 
signal  estimation  problem  with  if'  =  G~^oHg. 
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The  above  discussion  implies  that  unless  the  solution  to  the  signal  estimation  problem 
is  {GoH'^)  with  iT'  begin  the  solution  to  the  corresponding  ANC  problem,  the  optimal 
feedback  ANC  system  will  be  different  from  the  optimal  signal  estimation  system. 

In  general,  a  procedure  for  solving  the  single-channel  signal  estimation  problem  cannot 
be  used  to  solve  the  single-channel  feedback  ANC  problem,  since  a  procedure  for  solving  an 
un-constrained  minimization  problem  cannot  be  used  to  solve  a  constrained  minimization 
problem. 

If  the  minimizing  in  the  ANC  problem  can  be  commuted  with  G,  the  single-channel 
feedback  ANC  problem  can  be  solved  by  solving  a  related  two-input /single-output  signal 
estimation  problem.  The  ANC  system  in  this  case  can  be  equivalently  expressed  in  block 
diagram  form  as  depicted  in  Fig.  2.17  where  qa{t)  =  G{xa{t).  Referring  to  Fig.  2.17, 
the  single- channel  ANC  problem  in  this  case  can  be  restated  as  finding  the  system  H'^ 
that  makes  the  output  ea(t)  as  close  as  possible  to  the  desired  signal.  This  class  of  two- 
input/single-output  signal  estimation  systems  is  studied  in  the  next  Section. 
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Figure  2.17:  Equivalent  block  diagram  for  the  single-channel  feedback  ANC  system,  pro¬ 
vided  that  G  and  /f'  can  be  commuted. 


2.5.2  Relationship  between  Feedforward  Signal  Estimation  and  Feedfor¬ 
ward  ANC 

One  way  to  do  signal  estimation  is  to  first  generate  a  canceling  signal,  and  then  form  the 
estimate  as  the  sum  of  this  canceling  signal  and  one  of  the  observations.  We  shall  refers  to 
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this  approach  as  feedforward  signal  estimation. 


Figure  2.18:  Feedforward  signal  estimation  system. 

A  generic  block  diagram  for  a  single-channel  feedforward  signal  estimation  system  is 
depicted  in  Fig.  2.18.  As  depicted  in  Fig.  2.18,  this  system  is  a  two-input/single-output 
system  in  which  the  output  es(i)  is  the  sum  of  the  canceling  signal  Cs(t)  and  the  input  signal 
Xs(t).  Moreover,  the  canceling  signal  is  generated  solely  based  on  one  of  the  inputs,  namely 
9s(t). 

The  adaptive  noise  canceling  system  proposed  by  B.  Widrow  [52]  has  the  same  structure 
as  the  system  in  Fig.  2.18,  except  that  in  Widrow’s  system  is  continuously  adjusted. 

The  objective  of  the  feedforward  signal  estimation  system  is  to  make  es(t)  as  close  as 
possible  to  some  desired  signal  ds(t).  Specifically,  the  feedforward  signal  estimation  problem 
can  be  stated  as  follows: 

Find  the  single-input/single- output  system  Hg  that  minimizes 

J{{[xg{T)  +  Hsiqg{r)),dg{T)]:  0<r<T}).  (2.36) 

Next,  we  consider  the  single- channel  feedforward  ANC  system.  As  depicted  in  Fig. 
2.19,  a  single-channel  feedforward  ANC  system  is  a  two-input/single-output  system  where 
the  output  ea{t)  is  generated  as  the  sum  of  the  canceling  signal  Ca{t)  and  the  input  signal 
Xa{t).  Moreover,  the  canceling  signal  Ca{t)  is  generated  as  the  output  of  system  G,  and  it  is 
assumed  that  G  is  not  an  identity  system.  In  this  system,  laCO  and  Ca{t)  are  both  physical 
signals,  and  system  G  represents  the  effect  of  the  transducer  that  physically  generates  the 
canceling  signal. 
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Figure  2.19:  Feedforward  ANC  system. 


The  objective  of  the  single-channel  feedforward  ANC  system  is  to  make  ea[t)  as  close 
as  possible  to  some  desired  signal  da{i).  Specifically,  the  feedforward  ANC  problem  can  be 
stated  as  follows: 

Find  the  single-input/ single-output  system  Ha  that  minimizes 

J{{[xa{r)F{GoHa){qa{r)),da{r)]  :  0  <  r  <  T}) .  (2.37) 


Comparing  the  feedforward  signal  estimation  problem  to  the  feedforward  ANC  problem, 
we  see  that  the  feedforward  ANC  problem  can  be  viewed  as  a  constrained  version  of  the 
feedforward  signal  estimation  problem.  Specifically,  the  feedforward  ANC  problem  can  be 
stated  as  finding  Hs  that  minimizes  the  expression  in  eq.  (2.36)  subject  to  the  constraint 
that  Hs  must  be  of  the  following  form 

Hs  =  GoHa,  (2.38) 

for  a  given  system  G  and  some  system  Ha- 

We  now  compare  the  single-channel  feedforward  signal  estimation  system  to  the  single¬ 
channel  feedforward  ANC  system  in  terms  of  their  input/output  characteristics.  To  this 
end,  we  assume  that  both  systems  are  driven  by  the  same  inputs  and  that  the  desired  signal 
for  both  systems  is  the  same.  We  also  assume  that  each  system  is  designed  so  that  the  the 
distance  between  its  output  and  the  desired  signal  is  minimized. 
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Under  these  assumptions,  the  performance  of  the  ANC  system  is  never  better  than  the 
performance  of  the  corresponding  signal  estimation  system,  since  the  optimal  ANC  system 
is  obtained  by  minimizing  J{.)  subject  to  a  constraint,  and  the  optimal  signal  estimation 
system  is  obtained  by  minimizing  J(.)  without  any  constraints. 

Again  under  these  assumptions,  if  the  minimizing  happens  to  equal  (GoQ)  for  some 
system  Q,  the  optimal  ANC  system  will  be  identical  to  the  optimal  signal  estimation  system. 
Specifically  in  this  case,  the  minimizing  in  the  ANC  problem  is  system  Q.  For  example, 
if  G  is  invertible,  the  minimizing  Hs  can  be  expressed  as 

Hs  =  loHs  (2.39) 

=  Go{G-^oHs).  (2.40) 

This  implies  that  the  solution  to  the  feedforward  ANC  problem  can  be  expressed  in  terms  of 
the  solution  to  the  corresponding  feedforward  signal  estimation  problem  with  Ha  —  G~^oHs- 
The  above  discussion  implies  that  unless  the  solution  to  the  signal  estimation  problem  is 
(GoHa)  with  H'^  being  the  solution  to  the  corresponding  ANC  problem,  the  optimal  ANC 
system  will  be  different  from  the  optimal  signal  estimation  system. 

In  general,  a  procedure  for  solving  the  feedforward  signal  estimation  problem  can  not  be 
used  to  solve  the  feedforward  ANC  problem,  since  a  procedure  for  solving  an  un-constrained 
minimization  problem  cannot  be  used  to  solve  a  constrained  minimization  problem. 

If  the  minimizing  Ha  in  the  ANC  system  can  be  commuted  with  G,  the  feedforward 
ANC  problem  can  be  solved  by  solving  a  related  feedforward  signal  estimation  problem. 
Specifically  in  this  case,  the  minimizing  Ha  is  the  solution  to  the  signal  estimation  problem 
with  Xs{t)  =  Xa{t)  and  qs{t)  =  G{qa{t)). 

2.6  Summary 

Acoustical  principles  behind  ANC  and  several  commonly  used  control  strategies  for  ANC 
were  reviewed  in  this  chapter.  The  relationship  between  active  noise  cancellation  and  signal 
estimation  was  also  explored. 

Numerical  studies  were  presented  to  demonstrate  the  feasibility  of  minimizing  the  total 
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acoustic  potential  energy  in  an  enclosure.  These  studies  showed  that  at  low  frequencies, 
a  few  secondary  sources  can  be  used  to  obtain  substantial  reduction  in  the  total  acoustic 
potential  energy  in  an  enclosure  provided  that  two  or  three  modes  dominate  the  response 
of  the  enclosure.  In  another  study,  it  was  found  that  driving  the  acoustic  pressure  to  zero 
at  a  single  point  in  an  enclosure  produces  a  zone  of  quiet  around  this  point  which  has  a 
diameter  of  about  one  tenth  of  the  wavelength  of  the  excitation  frequency.  It  is  important 
to  recall  that  the  secondary  source  strengths  in  all  these  studies  were  calculated  based  on 
the  exact  knowledge  of  the  primary  source  strength  for  all  time.  However  in  a  practical 
ANC  system,  the  primary  source  strength  is  unknown  and  must  be  estimated  based  on  the 
available  measurements. 

Existing  control  strategies  for  ANC  were  characterized  as  pointwise  feedforward  con¬ 
trollers  or  as  pointwise  feedback  controllers. 

Feedforward  systems  achieve  noise  cancellation  by  exploiting  the  cross-correlation  be¬ 
tween  the  reference  signal  q(t)  and  the  unwanted  noise  n(t)  at  the  error  microphone.  For 
the  single- channel  feedforward  ANC  system,  we  showed  that  the  optimal,  in  the  minimum 
mean-square  sense,  canceling  signal  at  the  error  microphone  is  the  ELSE  of  —n{t)  based  on 

p"it)  =  q{t)*  gaii{i),  (2.41) 

where  is  the  aU-pass  part  of  the  transfer  function  from  the  input  of  the  canceling 

loudspeaker  to  the  value  of  the  canceling  field  at  the  location  of  the  error  microphone,  see 
equation  (2.28).  We  also  showed  that  this  formulation  can  be  used  to  compute  the  optimal 
performance  of  the  feedforward  system  based  on  the  second-order  statistics  of  {n{t),q(t)}. 

Feedback  systems  achieve  noise  cancellation  by  exploiting  the  auto-correlation  of  the 
noise  n(t)  at  the  error  microphone.  For  the  single-channel  feedback  ANC  system,  we  showed 
that  the  optimal  canceling  signal  at  the  error  microphone  is  the  ELSE  of  —n{t)  based  on 

p'{i)  =  n(t)*  gau{t),  (2.42) 

where  Gaiiis)  is  the  aU-pass  part  of  the  transfer  function  from  the  input  of  the  canceling 
loudspeaker  to  the  output  of  value  of  the  canceling  field  at  the  location  of  the  error  micro¬ 
phone,  see  equation  (2.26).  We  also  showed  that  this  formulation  can  be  used  to  compute 
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the  optimal  performance  of  the  feedback  system  based  on  the  second-order  statistics  of  n{t). 

Using  the  above  analysis  we  can  compare  the  optimal  performance  of  the  feedforward 
system  to  the  optimal  performance  of  the  feedback  system.  Based  on  this  analysis,  we 
expect  the  feedforward  system  to  achieve  low  levels  of  attenuation  whenever  the  reference 
signal  is  weakly  correlated  with  the  noise  at  the  error  microphone. 

The  main  drawback  of  the  feedforward  systems  is  the  problem  of  how  to  acquire  a 
reference  signal  that  is  highly  correlated  with  the  unwanted  noise  and  at  the  same  time 
avoid  acoustic  feedback  from  the  output  of  the  canceling  loudspeaker  to  the  reference  signal. 

The  main  drawback  of  analog  feedback  systems  is  that  these  systems  can  not  be  tuned  to 
take  advantage  of  the  statistical  characteristics  of  the  unwanted  noise.  The  main  drawback 
of  the  digital  feedback  system  of  Graupe  [26]  is  that  it  can  only  be  used  when  the  transfer 
function  from  the  canceling  loudspeaker  to  the  error  microphone  contains  no  delay.  This 
severely  limits  the  applications  in  which  this  system  can  be  used.  Furthermore,  a  multiple- 
microphone  version  of  the  system  proposed  by  Graupe  does  not  exist. 

In  Section  2.5,  we  explored  the  relationship  between  signal  estimation  and  active  noise 
cancellation.  The  systems  studied  in  this  Section  were  free  to  be  linear  or  non-linear  and  the 
cost  functions  involved  were  quite  general.  By  structuring  the  signal  estimation  problem  in  a 
particular  way,  we  were  able  to  illustrate  the  close  relationship  between  single- channel  signal 
estimation  and  single- channel  feedback  active  noise  cancellation  and  the  close  relationship 
between  feedforward  signal  estimation  and  and  feedforward  active  noise  cancellation.  In 
each  case,  we  showed  that  the  active  noise  cancellation  problem  can  be  viewed  as  a  con¬ 
strained  version  of  the  signal  estimation  problem. 

In  the  context  of  pointwise  active  noise  cancellation,  our  goal  is  to  develop  a  general 
framework  for  designing  feedback  active  noise  cancellation  systems  that  can  be  used  in  a 
wide  range  of  applications.  Our  formulation  can  be  used  in  single-microphone  as  well  as 
multiple-microphone  configuration,  and  the  resulting  algorithms  will  take  advantage  of  the 
the  statistical  characteristics  of  the  unwanted  noise. 

We  will  then  extend  the  above  framework  to  the  case  where  the  performance  criterion 
is  distributed.  The  performance  criterion  in  this  case  is  the  total  acoustic  energy  in  an 
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enclosure.  However,  the  measurements  are  still  made  using  ordinary  microphones,  and  the 
control  signal  is  still  used  to  drive  an  ordinary  loudspeaker.  Most  importantly,  the  control 
signals  are  generated  causally  based  only  on  the  outputs  of  the  microphones. 
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Chapter  3 


Pointwise  Optimal  Feedback 
Control  for  ANC 


3.1  Introduction 

In  this  chapter  we  develop  an  optimal  feedback  controller  for  attenuating  the  acoustic  pres¬ 
sure  at  the  locations  of  a  finite  number  of  error  microphones.  Our  controller  is  designed  to 
achieve  the  same  acoustic  objective  as  the  multiple-channel  feedforward  controller  in  [18]; 
however,  our  controller  generates  the  control  signal  based  on  the  outputs  of  the  error  micro¬ 
phones,  and  not  based  on  a  separate  reference  signal.  Therefore,  the  feedback  ANC  system 
developed  in  this  chapter  achieves  the  same  acoustic  objective  as  the  feedforward  ANC 
system  in  [18],  without  the  problems  associated  with  acquiring  an  appropriate  reference 
signal. 

An  important  feature  of  the  optimal  feedback  controller  is  that  it  takes  full  advantage 
of  the  statistics  of  the  noise.  LMS-type  feedforward  controllers  used  in  ANC  systems  take 
advantage  of  the  statistical  characteristics  of  the  noise  to  be  canceled.  However,  existing 
feedback  ANC  controllers  do  not  take  advantage  of  the  noise  statistics  (except  under  very 
restrictive  assumptions  [26]). 

We  present  a  procedure  for  designing  optimal  feedback  controllers  for  pointwise  ANC 
that  can  be  used  under  very  general  conditions.  The  error  microphones  are  placed  at  those 
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locations  where  noise  attenuation  is  desired,  and  the  canceling  loudspeaker  is  placed  in 
the  vicinity  of  these  microphones.  Given  a  specific  arrangement  of  the  error  microphones 
and  the  canceling  loudspeaker,  and  given  a  specific  statistical  model  for  the  noise  at  these 
microphones,  our  design  procedure  determines  the  optimal  feedback  controller.  Specifically, 
our  design  procedure  finds  the  controller  that  minimizes  a  performance  criterion  involving 
the  mean-square  sum  of  the  residual  signals  at  the  error  microphones  and  the  mean-square 
sum  of  the  control  signal. 

Our  approach  is  to  model  the  residual  signals  at  the  microphones  as  the  sum  of  the 
outputs  of  two  linear  system.  The  unwanted  noise  signals  at  the  microphones  are  modeled 
as  the  outputs  of  a  multiple-input /multiple- output  linear  system  driven  by  a  white  process. 
Similarly,  the  canceling  signals  at  the  microphones  are  modeled  as  the  outputs  of  a  single- 
input/multiple-output  linear  system  driven  by  the  control  signal.  Finally,  the  residual 
signcils  at  the  microphones  are  modeled  as  the  sum  of  the  outputs  of  these  two  linear  systems. 
We  show  that  the  control  signal  that  minimizes  a  certain  class  of  pointwise  performance 
criteria  is  a  linear  feedback  of  the  estimated  states  of  these  two  linear  systems.  These 
state  estimates  are  computed  using  a  Kalman  filter,  and  the  feedback  gain  is  computed  by 
iterating  a  Riccati  equation.  The  single-microphone  case  is  presented  first,  and  the  results 
are  then  extended  to  the  multiple-microphone  case. 

In  Section  3.4,  we  show  that  the  optimal  feedback  controller  developed  in  this  chapter 
can  be  thought  of  as  a  generalization  of  the  optimal  feedforward  controller.  Specifically, 
the  optimal  feedforward  controller  generates  the  canceling  signal  based  on  all  the  available 
measurement,  but  the  optimal  feedforward  controller  generates  the  canceling  signal  based 
on  a  proper  subset  of  the  available  measurements. 

The  control  algorithms  developed  in  this  chapter  are  discrete-time  algorithms.  Unless 
stated  otherwise,  we  assume  throughout  this  chapter  that  the  noise  to  be  attenuated,  the 
canceling  signal,  and  all  other  intermediate  signals  are  appropriately  sampled.  Unless  stated 
otherwise,  throughout  this  chapter  “t”  represents  the  normalized  sampling  time. 
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3.2  Single-Microphone  Optimal  Feedback  Controller  for  ANC 


A  single- microphone  optimal  feedback  controller  for  ANC  is  developed  in  this  Section.  The 
system  utilizes  one  error  microphone  and  one  canceling  loudspeaker.  The  microphone  is 
placed  at  the  point  where  noise  attenuation  is  desired,  and  the  loudspeaker  is  placed  as 
close  as  possible  to  this  microphone.  The  acoustic  objective  is  to  attenuate  the  pressure  at 
the  location  of  the  microphone.  The  loudspeaker  produces  a  secondary  field  that  interferes 
destructively  with  the  primary  field  at  the  microphone,  generating  a  zone  of  quiet  around 
the  microphone  as  reported  in  [15]. 

A  generic  single-microphone  feedback  ANC  system  is  depicted  in  Fig.  3.1  where  m{t)  is 
the  microphone  output,  and  r(t)  is  the  input  of  the  canceling  loudspeaker.  The  microphone 
measures  the  sum  of  the  unwanted  noise  n{t)  and  the  canceling  signal  c{t).  The  objective 
is  to  generate  r{t),  based  on  the  measurements  of  m(t),  so  that  the  power  of  the  residual 
signal  at  the  microphone  is  minimized. 

Throughout  this  chapter  we  rely  on  the  following  model  for  the  microphone.  The  output 
of  the  microphone  is  assumed  to  be  the  sum  of  the  acoustic  pressure  at  the  location  of 
the  microphone  and  the  microphone  measurement  noise.  According  to  this  model,  the 
only  distortion  introduced  by  the  microphone  is  this  additive  measurement  noise.  This 
microphone  model  is  quite  accurate  for  most  practical  pressure  microphones  at  frequencies 
below  lOOOHz  [6]. 

The  block  diagram  for  the  single- microphone  ANC  system  that  we  propose  and  develop 
in  this  Section  is  depicted  in  Fig.  3.2.  The  system  G{z)  represents  the  overall  transfer 
function  from  the  canceling  loudspeaker  input  r(t)  to  the  value  of  the  canceling  field  at  the 
location  of  the  microphone  c{t)  and  incorporates  the  transfer  function  for  the  loudspeaker 
and  the  propagation  path  between  the  loudspeaker  and  the  microphone.  The  microphone 
output  m{t)  is  the  sum  of  the  unwanted  noise  n{t),  the  canceling  signal  c(t),  and  the 
microphone  measurement  noise  u(t). 

Our  betsic  strategy  for  solving  the  ANC  problem  is  to  view  this  problem  as  a  stochastic 
optimal  control  problem  with  n(t)  being  the  stochastic  disturbance  and  G{z)  being  the  plant. 
From  this  point  of  view,  the  objective  is  to  keep  the  output  of  the  plant,  €{t)  =  c(t)  +  n{t), 
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as  close  to  zero  as  possible.  Note  that  e{t)  is  the  value  of  the  residual  acoustic  field  at  the 
location  of  the  microphone,  i.e.  the  sum  of  the  noise  field  at  the  location  of  the  microphone 
and  canceling  field  at  the  location  of  the  microphone. 


Noise  Source 


Figure  3.1:  Generic  single- mi  crop  hone  feedback  active  noise  cancellation  system. 

3.2.1  Model  Specification 

The  stochastic  optimal  control  formulation  requires  a  disturbance  model  and  a  plant  model. 
We  model  the  noise  n{t)  as  the  output  of  a  pole-zero  system  function  driven  by  a  white 
process: 

9  j 

^kn{t  -  k)  +  -  k)  j  <  q,  (3.1) 

k=l  k^l 

where  w{t)  is  a  zero  mean  unit  variance  white  Gaussian  process,  i.e.  n{t)  is  modeled  as  an 
ARMA  process.  The  plant  model,  corresponding  to  the  transfer  function  from  the  input 
of  the  canceling  loudspeaker  to  the  value  of  the  canceling  field  at  location  of  the  error 
microphone,  is  the  pole-zero  system  function 

Ylk=ibkZ-^  ^ 
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u(t) 


Figure  3.2:  Block  diagram  of  the  single-microphone  feedback  active  noise  cancellation  sys¬ 
tem. 
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We  represent  the  system  of  Fig.  3.2  in  state-space  form  in  three  steps.  First  we  will 
derive  a  state-space  representation  for  the  disturbance  model  eq.  (3.1),  and  then  we  will 
derive  a  state-space  representation  for  the  plant  model  eq.  (3.2).  These  two  representations 
are  finally  combined  to  get  the  overall  state-space  representation  for  the  system  in  Fig.  3.2. 

The  disturbance  model  can  be  expressed  in  state-space  form  as  [34,  3,  2]. 


xi{t)  =  Xd{t  -1)  +  Ld  w{t  -  1) 
n{t)  =  Hd  Xd{t), 


(3.3) 

(3.4) 


where  Xd{t)  is  the  {q  x  1)  state  vector  for  the  disturbance  model;  ^d  is  the  corresponding 
{q  X  q)  transition  matrix 


— ai  1  0 

-02  0 

-a,-!  0  ■  •  ■ 


-O',  0 


0 


1 

0 


Ld  is  the  {q  X  1)  vector 


Ld  =  [Pi  •••  /?,  O-.-O]^; 

and  Hd  is  the  (1  x  g)  unit  vector 

Hd  =  [lO-.-O]. 

Similarly,  the  plant  model  can  be  expressed  in  state-space  form  as  [34]: 

Xp{t)  =  Xp{t  -\)-\-Gp  r{t  -l)  +  Lp  u{t  -  1) 
c(t)  =  Hp  a;p(t). 


(3.5) 


(3.6) 

(3.7) 


where  Xp(t)  is  the  (n  x  1)  state  vector  for  the  plant  model;  u{t)  is  an  {n  x  1)  vector  of 
random  variables  such  that 


E  |•u(<)^i(T)^|  =  S{t  -  r)/„xu, 
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with  u{t)  independent  of  $p  is  the  (n  x  n)  transition  matrix  of  the  plant 


-ai  1  0 


-a2  0 


-On-l  0 


-an  0 


Gp  is  the  (n  X  1)  vector 


Lp  is  the  (n  x  n)  matrix 


Gp  —  [bi  ■■  -  bm  0  •  •  -0]^; 


and  Hp  is  the  (n  x  1)  unit  vector 


Lp  —  a-nlnxm 


Hp  =  [\ 


The  state-space  model  for  the  overall  system  is  now  obtained  by  combining  (3. 3-3.4) 
with  (3.6-3.7): 


x{t)  =  $  a;(t  —  1) -1- (j  7-(t  —  1) -H  i 


m{i)  —  Hx{t) -]r  v{t), 


w{t  —  1) 
u(t  —  1) 


(3.10) 


where  x{t)  is  the  state  vector  of  the  overall  system 


$  is  the  corresponding  (g  +  n)  x  {q  +  n)  transition  matrix 


(3.11) 


0 


0 


(3.12) 
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G  is  the  (9  +  J^)  X  1  vector 


0 


L  is  the  {q  +  n)  X  {q-\-  n)  matrix 


and  H  is  the  1  x  (g  +  n)  vector 


Ld  0 

0  Lp 


(3.13) 


H  =  [Hd  Hp]. 

The  measurement  noise  in  the  error  microphone  is  modeled  by  the  white  Gaussian  process 
v{t)  whose  variance  is  <7^,  and  v{t)  is  independent  of  w{t)  and  u{t). 


3.2.2  The  Optimal  Feedback  Controller 

The  goal  of  the  ANC  system  is  to  choose  {r(f)  :  t  =  0, . . ., iV  —  1}  so  that  the  following 
performance  measure  is  minimized: 


(3.14) 


To  be  causally  realizable,  the  choice  for  r(t)  must  be  a  function  of  the  available  measure¬ 
ments  at  time  “t”,  i.e.  it  must  be  a  function  of  {m(r)  :  r  <  t  —  1}.  The  first  summand  in 
(3.14)  represents  the  power  of  the  residual  signal  at  the  microphone,  and  the  second  sum¬ 
mand  represents  the  power  expended  by  the  control  signal.  Here,  p  is  a  weighting  factor 
which  is  chosen  to  achieve  a  trade  off  between  control  effort  and  noise  attenuation.  Recalling 
that 

e{t)  =  H  x{t),  (3.15) 


we  can  rewrite  J{N)  in  terms  of  the  state  vector  x{i)  as 

r7V-i  > 

J(A)  =  e\Y^  {x'^{t)H^Hx(t)  +  p  -h  x^{N)H'^Hx{N) 


(3.16) 


t=o 
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It  is  important  to  observe  that  (3.16)  is  the  expected  value  of  a  quadratic  function  of  the 
state  vector  x{t)  and  the  control  signal  r{t).  Following  standard  terminology  from  control 
literature,  we  shall  refer  to  this  minimization  problem  as  the  “Finite  Horizon  Optimal  Linear 
Quadratic  Control  Problem”. 

Since  the  performance  measure  J{N)  is  the  expected  value  of  a  quadratic  form  in  the 
state  and  the  control,  the  minimizing  choice  for  {r(t)  :  t  =  0, . . .,  iV  -  1}  is  a  linear  feedback 
of  the  optimal  state  estimates  [24].  The  optimal  state  estimates  can  be  computed  using 
a  Kalman  filter,  and  the  feedback  gain  can  be  computed  by  iterating  a  Riccati  difference 
equation.  To  this  end,  let  us  denote  by 

i(tlt)  =  £  {a;(t)|m(l), . . .,  m(t)}  (3-17) 

the  state  estimate  based  on  observations  up  to  time  “t”,  and  by 

P{t\t)  =  E  |[a:(tlt)  -  a:(t)][i(t|t)  -  a:(t)]^|m(l), . . . ,  m(t)|  (3.18) 

the  associated  error  covariance  matrix. 

Then  using  the  standard  Kalman  filter  formulation,  i(t|t)  and  Pit\t)  can  be  computed 


sequentially  in  time,  in  two  stages,  as  follows: 

Propagation  Equations: 

a:(t|t  —  1)  =  #  i(t  -  l|t  -  1)  +  G  ti(t  -  1)  (3.19) 

P{t\t-1)  =  $  P(t- 1|/- 1)  +  (3.20) 

Updating  Equations: 

x{t\t)  =  x{t\t  —  1)  +  K f{t)[m(t)  —  H  x{t\t  —  1)]  (3.21) 

Pit\t)  =  [I-Kj{t)H]P{t\t-l),  (3.22) 

where  K f{t)  is  the  Kalman  gain  given  by 


The  minimizing  choice  for  {r{t)  :  t  =  0, . . .,  iV  -  1}  is  then  given  by 
r{N-j)  =  -[GP{j-l)G  +  p]-^G^P{j-l)^  x{N-j\N-j)  j  =  l,...,N 

=  K,ij-1)  xiN-j\N-j),  (3.24) 
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where  the  (m  +  n)  x  (m  +  n)  matrix  P{j)  is  obtained  by  iterating  the  following  Riccati 
difference  equation: 

P{j  +  1)  =  )$  -  (3.25) 

with  the  initial  condition  P(0)  =  P^P.  The  resulting  minimum  value  of  the  performance 
index  J{N)  can  also  be  computed  based  on  P{N)  [35].  The  block  diagram  for  the  resulting 
ANC  system  is  depicted  in  Figure  3.3. 


. . . . . . 


Optimal  Feedback  Controller 


Figure  3.3:  Block  diagram  of  our  single- microphone  feedback  ANC  system. 


Infinite  Horizon  Case 

If  we  consider  the  above  problem  with  both  time-invariant  system  matrices  and  with  con¬ 
stant  weighting  factor  p,  and  if  we  allow  N  to  approach  infinity,  a  new  performance  index 
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can  be  defined  as 


J=  lim  ^JiN). 

In  this  case,  we  ask  for  the  choice  of  {r(t)  :  <  =  0, 1, . . .}  that  minimizes  J.  This  is  a  well 
posed  problem  provided  that  [$,G]  is  stabilizable,  and  is  detectable.  In  this  case, 

P{j)  converges  to  a  constant  matrix  P(oo)  which  is  the  unique  positive  definite  solution  of 
the  following  algebraic  Riccati  equation 

P{oo)  =  $^P(oo)$  -  #^P(oo)G[G'^P(oo)G  +  p]"^G^P(oo)$  +  H'^H.  (3.26) 

In  this  case,  the  minimizing  choice  for  r(t)  is  given  by 

r(t)  =  -[G^P(oo)G  +  p]-^G^P(oc)^  x(tjt)  (3.27) 

=  Kc(oc)  (3.28) 

Note  that  the  control  gain,  Kc(oo),  is  time-invariant. 

Similarly,  the  error  covariance  matrix,  P(t\t  —  1)  in  the  Kalman  filter,  converges  to 
a  constant  matrix  P(oo)  which  is  the  unique  positive  definite  solution  of  the  following 
algebraic  Riccati  equation 

P(oo)  =  $F(oo)^^  -  ^P(oo)H^[HP(oo)H^  +  H P(oo)^'^  +  LL^ , 

provided  that  [^,T]  is  stabilizable  and  [$,ir]  is  detectable.  The  Kalman  filter  equations 
reduce  to: 

Propagation  Equation: 

x(t\t  —  1)  =  $  x(t  —  l|t  —  1)  -f  G  u(t  —  1)  (3.29) 

Updating  Equation: 

x(t\t)  =  x(t\t  —  1)  +  K /(oo)[m(t)  -  H  x(t\t  —  1)],  (3.30) 

where  K f(oo)  is  the  Kalman  gain  given  by 

Note  that  the  Kalman  filter  in  this  case  is  linear  and  time-invariant.  Recalling  that  the 
control  gain  Kc(oo)  is  also  time-invariant,  we  conclude  that  the  controller  becomes  a  linear 
time-invariant  system  in  this  case. 
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3.2.3  Frequency-Domain  Perspective 


Our  formulation  of  the  ANC  problem  so  far  has  been  entirely  in  the  time  domain.  Addi¬ 
tional  insight  can  be  gained  by  looking  at  the  closed-loop  ANC  system  in  the  frequency 
domain.  For  this  purpose,  we  shall  concentrate  on  the  infinite  horizon  case  in  which  the 
overall  system  is  linear  and  time-invariant,  and  all  the  signals  involved  are  jointly  station¬ 
ary.  Furthermore,  we  assume  that  the  objective  is  to  minimize  the  pressure  fluctuations  at 
the  error  microphone,  i.e.  we  assume  that  the  weighting  factor  for  control  effort,  p  in  eq. 
(3.14),  is  zero. 

Let  H{z)  denote  the  transfer  function  of  the  controller  in  Fig.  3.3.  Since  the  measure¬ 
ment  noise  v(t)  is  white  and  G{z)  contains  at  least  one  delay,  the  causal  controller  that 
minimizes  E{m^(t)}  is  the  same  as  the  causal  controller  that  minimizes  E{e^{t)}  (note 
that  m{t)  =  e{t)  -f  v{t)).  Therefore,  we  can  characterize  the  optimal  feedback  controller  as 
the  controller  that  minimizes  E{m\t)}  in  Fig.  3.3.  Referring  to  this  figure,  we  define  the 
closed-loop  transfer  function  T{z)  as 


,  and  we  define  z(t)  as 


T{z) 


1 

1  -  Hiz)  Giz) 


z(t)  =  n(i)  +  v(t}. 


(3.32) 


(3.33) 


The  power  spectrum  of  the  microphone  output  can  be  expressed  in  terms  of 

the  power  spectrum  of  z(t),  i.e. 

|2 


Pmmieh  =  |r(e^“)| 


(3.34) 


Equations  (3.32)-(3.34)  show  how  the  controller  affects  the  power  spectrum  of 

the  microphone  output.  To  minimize  the  variance  of  m{t)  or  equivalently  to  minimize 
the  variance  of  e(t),  we  need  to  make  |r(e-'‘^)|  small  over  those  frequencies  for  which  the 
unwanted  noise  has  significant  power.  This  is  precisely  what  the  optimal  controller  does. 
Specifically,  the  optimal  controller  uses  the  noise  model  to  determine  those  frequencies 
for  which  the  unwanted  noise  signal  has  significant  power,  and  it  chooses  H{z)  so  that 
the  magnitude  of  the  resulting  closed-loop  transfer  function,  |T(e^“)|,  is  small  over  these 
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frequencies.  From  this  point  of  view,  it  is  clear  that  the  more  narrow-band  the  spectrum  of 
n(t),  the  higher  the  attenuation  level  achieved  by  the  feedback  controller.  An  example  of 
this  will  be  presented  in  the  next  section. 

Having  a  good  model  for  the  unwanted  noise  is  crucial  for  achieving  high  attenuation 
levels.  In  the  absence  of  an  accurate  model  for  the  the  unwanted  noise  (i.e.  without 
knowing  at  which  frequencies  the  power  of  the  unwanted  noise  is  concentrated),  one  is 
forced  to  assume  that  the  spectrum  of  the  noise  is  flat  over  a  wide  range  of  frequencies. 
This  results  in  an  ANC  system  with  low  attenuation  over  a  wide  range  of  frequencies  (e.g., 
noise  canceling  headphones  developed  by  Wheeler  [51]).  On  the  other  hand,  in  cases  for 
which  the  unwanted  noise  is  relatively  narrow  band  and  the  noise  is  well  modeled,  the  active 
noise  cancellation  system  can  be  designed  to  achieve  very  high  attenuation  levels.  This  is 
done  by  choosing  the  controller  such  that  the  closed-loop  system  has  high  attenuation  over 
the  narrow  band  of  frequencies  where  the  power  of  the  noise  is  concentrated. 

Analysis  of  the  feedback  controller  in  chapter  2  can  also  be  used  to  explain  why  relatively 
narrow  band  noise  is  highly  attenuated  by  the  optimal  feedback  controller.  Recall  that  this 
analysis  was  done  in  continuous  time;  hence,  “f’  represent  the  continuous  time  variable  in 
this  paragraph.  Referring  to  Fig.  2.3  of  Chapter  2,  it  was  shown  that  the  optimal  canceling 
signal  in  the  single- channel  feedback  ANC  system  is  the  LLSE  of  -n{t)  based  on 

p'(t)  =  n{t)*gauit),  (3.35) 

where  Gaii{s)  is  the  aU-pass  part  of  G{s).  The  mean-square  estimation  error  will  be  small 
whenever  n{t)  can  be  estimated  accurately  based  on  For  a  relatively  narrow  band 

n(f),  Gaii{s)  can  be  approximated  as  a  delay.  Therefore  in  the  narrow  band  case,  the 
problem  reduces  to  predicting  n{t)  based  on  p'{t)  =  n{t  -  r),  where  r  is  the  delay.  More 
importantly,  the  mean-square  estimation  error  will  be  small  with  a  narrow  band  n(i),  since 
a  narrow  band  signal  is  highly  predictable. 

3.2.4  Frequency- Weighted  Optimal  Controller  for  ANC 

Assuming  time-invariant  system  matrices  and  stationary  acoustic  noise,  we  shall  drive  a 
feedback  controller  that  minimizes  a  cost  function  involving  the  frequency-weighted  power 
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of  the  residual  signal  at  the  microphone.  Referring  to  Fig.  3.2,  the  objective  in  this  section 
is  to  minimize  ^ 

=  ^Im^  jE  I  [el{t)  +  p  r\t))  +  e^(iV)| ,  (3.36) 

where  =  e{t)  *  s{t),  with  s{t)  being  a  pre-specified  weighting  filter.  The  relationship 
between  e(t)  and  can  be  expressed  in  state-space  form  as: 


Xyj{t)  =  Xu,{t  -  1)  +  Gy,e{t) 


(3.37) 

(3.38) 


The  weighting  filter  s{t)  might  represent  the  perceived  loudness  of  sound  at  different 
frequencies.  In  this  context,  the  filter  s(t)  is  typically  referred  to  as  an  A-weight  filter. 

The  state-space  representation  eqs.  (3.37)-(3.38)  of  the  weighting  filter  can  be  combined 
with  the  state-space  representation  eqs.  (3.9)-(3.10)  of  the  original  acoustic  system.  The 
state-space  representation  for  the  combined  system  is 


Xc{t)  =  Xc{t  -  1)  -1-  Gc  ^(<  -  1)  +  Tc 


m{t)  =  He  Xc{t)  -f-  u(t). 


w{t  —  1) 
u{t  —  1) 


(3.39) 


(3.40) 


where  Xc{t)  is  the  state  vector  of  the  combined  system 


Xc{t)  = 


is  the  corresponding  transition  matrix 


x^{t) 


Gc  is  the  column  vector 


GyjH 


(3.41) 


(3.42) 


(3.43) 
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Lc  is  the  matrix 


(3.44) 


and  He  is  the  row  vector 


Ha  =  [H  o]. 


(3.45) 


Next,  we  express  the  performance  criterion  J^o  in  terms  of  the  state  vector  Xc{t).  For 
this  purpose,  we  define 


J{N)  =  £  I  E  (e^ (t)  +  p  r\i))  +  el{N)  I . 


Recalling  that 


(3.46) 


(3.47) 


we  can  rewrite  J{N)  in  terms  of  the  state  vector  Xc{t)  as 


riv-i  'j 

J(iV)  =  F;|E  {^l{t)C'^C,Xe{t)+  pT\t))+xl{N)ClC.Xe{N)Y  (3-48) 

where  C*  =  [o  H-uj],  Note  that  J{N)  is  the  expected  value  of  a  quadratic  function  of  the 
state  vector  Xc{t)  and  the  control  signal  r(t). 

Since  J^o  =  limjv->oo  the  choice  of  {r(t)  :  t  >  0}  that  minimizes  J^o  is  a  linear 

feedback  of  the  optimal  state  estimates  of  the  combined  system,  provided  that  is 

stabilizable  and  He]  is  detectable.  The  exact  equations  for  the  controller  are  obtained 
following  steps  identical  to  the  ones  followed  to  get  the  infinite  horizon  controller  in  Section 
3.2.2. 


3.2.5  Implementation  and  Performance  of  the  Optimal  Feedback  Con¬ 
troller 

Two  sets  of  experiments  were  performed  to  evaluate  the  performance  of  the  optimal  single¬ 
microphone  feedback  ANC  system.  In  the  first  set,  computer  simulations  on  recorded 
aircraft  noise  were  used  to  compare  the  performance  of  the  optimal  feedback  system  to  the 
performance  of  the  feedforward  system  proposed  in  [10].  The  second  set  of  experiments 
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were  carried  out  to  measure  the  performance  of  the  hardware  implementation  of  the  op¬ 
timal  feedback  system  in  the  context  of  aircraft  noise.  The  performance  of  the  hardware 
implementation  of  the  optimal  feedback  system  was  also  compared  to  the  performance  of  a 
commercial  analog  feedback  system  built  by  the  BOSE  Corporation. 

Our  feedback  system  and  the  feedforward  system  in  [10]  were  evaluated  using  a  set  of 
noise  canceling  headphones.  These  headphones  were  equipped  with  a  reference  microphone, 
an  error  microphone,  and  a  canceling  loudspeaker.  The  reference  signal  for  the  feedforward 
controller  was  obtained  as  the  output  of  the  reference  microphone  that  was  attached  to  the 
outside  of  the  headphones,  and  the  error  microphone  was  placed  inside  the  headphones.  The 
canceling  loudspeaker  was  also  placed  inside  the  headphones,  about  1.2  centimeters  away 
from  the  error  microphone.  In  this  case,  the  transfer  function  G{z)  of  Figure  3.2  corresponds 
to  the  transfer  function  from  the  canceling  loudspeaker  to  the  error  microphone  inside  the 
headphones. 

A  separate  set  of  experiments  was  performed  initially  to  estimate  the  transfer  function 
G(z)  and  the  coefficients  of  the  disturbance  model.  It  was  found  that  an  ARMA  model  with 
three  zeros  and  four  poles  produces  a  very  good  approximation  to  the  transfer  function  from 
the  canceling  loudspeaker  to  the  error  microphone.  The  disturbance  model  was  an  all-pole 
system  function  of  order  four,  and  the  coefficients  of  this  model  {a'j.s  in  (3.1))  were  estimated 
using  the  algorithm  in  [40].  All  the  controllers  discussed  in  this  section  are  designed  based 
on  this  estimated  plant  model  and  this  estimated  disturbance  model. 

Simulations 

The  performance  of  the  feedforward  system  of  Burgess  [10]  was  compared  to  the  perfor¬ 
mance  of  the  infinite  horizon  feedback  controller  developed  in  Section  3.2.2.  In  computer 
simulations,  the  feedforward  system  of  Burgess  was  able  to  attenuate  the  noise  generated  by 
a  propeller  aircraft  by  about  17dB.  This  result  was  obtained  with  a  350-tap  finite  impulse 
response  controller  and  using  recordings  of  the  noise  made  at  the  reference  microphone  and 
at  the  error  microphone.  The  sampling  rate  was  6KHz.  Using  the  same  recorded  noise 
at  the  error  microphone,  our  infinite  horizon  feedback  controller  was  able  to  attenuate  the 
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noise  power  at  the  error  microphone  by  about  26dB.  Referring  to  Fig.  3.2,  this  attenuation 
was  calculated  according  to  the  following  equation: 


Un- Weighted  Attenuation  =  lOlogjo 


E{n\t)y 


(3.49) 


The  feedforward  controller  and  the  feedback  controller  were  designed  based  on  the  exact 
knowledge  of  the  plant  model  and  the  disturbance  model.  The  power  spectrum  of  the 
unwanted  noise  and  the  power  spectrum  of  the  residual  signal  obtained  with  the  infinite 
horizon  feedback  controller  are  depicted  in  Fig.  3.4(a). 

Another  set  of  simulations  was  performed  to  illustrate  the  performance  of  the  frequency- 
weighted  optimal  feedback  controller.  Again,  we  used  the  same  propeller  aircraft  noise  that 
was  used  in  the  simulations  described  in  the  previous  paragraph.  The  frequency  response 
of  the  weighting  filter,  s{t),  that  was  used  in  the  design  of  the  controller  is  depicted  in  Fig. 
3.5.  This  filter  models  the  variations  in  perceived  loudness  of  sound  at  different  frequencies. 
According  to  this  model,  a  pure  tone  of  unit  power  at  IKHz  sounds  70dB  louder  than  a  pure 
tone  of  unit  power  at  lOHz.  The  power  spectrum  of  the  unwanted  noise  at  the  microphone 
and  the  power  spectrum  of  the  residual  signal  obtained  with  the  frequency- weighted  optimal 
feedback  controller  are  depicted  in  Fig.  3.4(b).  The  un-weighted  attenuation  in  this  case  is 
about  20dB.  Referring  to  Fig.  3.2,  let  us  define  the  weighted  attenuation  as 


Weighted  Attenuation  =  10  log^o 


E{el{t)} 

E{nUt)y 


(3.50) 


where  e.w{t)  =  e{t)  *  s{t)  and  n^{t)  =  n{t)  *  s(t).  The  weighted  attenuation  in  this  case 
is  about  16dB.  For  comparison,  the  weighted  attenuation  corresponding  to  Fig.  3.4(a)  is 
about  lldB.  We  see  that  incorporation  of  the  weighting  filter  in  the  design  of  the  feedback 
controller  improves  the  weighted  attenuation  by  about  5dB. 


Hardware  Implementation 

The  infinite  horizon  feedback  controller,  designed  without  frequency  weighting,  was  imple¬ 
mented  in  hardware  using  a  single  AT&T  DSP32C  chip  with  6KHz  sampling  rate.  The 
resulting  noise  canceling  headphones  were  placed  in  an  enclosure,  and  propeller  aircraft 
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noise  was  injected  into  this  enclosure.  The  headphones  were  able  to  attenuate  the  noise 
power  at  the  error  microphone  by  25dB. 

The  spectrum  of  the  unwanted  noise  at  the  output  of  the  error  microphone  and  the 
spectrum  of  the  attenuated  noise  at  the  output  of  the  error  microphone  obtained  with  our 
feedback  controller  are  depicted  in  Fig.  3.6(a).  The  closed-loop  transfer  function  realized 
by  our  controller,  is  depicted  in  Fig.  3.6(b).  As  expected  from  the  discussion  in 

Section  3.2.3,  the  closed-loop  system  has  high  attenuation  over  the  90Hz-130Hz  band  where 
most  of  the  power  of  the  unwanted  noise  signal  is  concentrated. 

In  Figure  3.7,  the  performance  of  our  optimal  feedback  controller  is  compared  to  that 
of  the  analog  feedback  controller  developed  commercially  by  BOSE  Corporation.  The  spec¬ 
trum  of  the  original  noise  at  output  of  the  error  microphone  is  depicted  in  Fig.  3.7(a).  The 
spectrum  of  the  attenuated  noise  at  the  output  of  the  error  microphone  obtained  with  our 
controller  and  the  spectrum  of  the  attenuated  noise  at  the  output  of  the  error  microphone 
obtained  with  the  analog  controller  are  depicted  in  Fig.  3.7(b).  The  overall  attenuation  ob¬ 
tained  with  the  analog  feedback  system  is  17dB,  and  the  overall  attenuation  obtained  with 
the  optimal  feedback  system  is  25dB.  Looking  at  Fig.  3.7,  we  see  that  the  optimal  feedback 
controller  achieves  higher  attenuation  than  the  analog  controller  over  the  90Hz-130Hz  band 
where  most  of  the  energy  of  the  unwanted  noise  is  concentrated.  On  the  other  hand,  the 
optimal  feedback  controller  achieves  lower  attenuation  than  the  analog  controller  over  the 
350Hz-700Hz  band  where  the  original  noise  does  not  have  significant  power. 

The  main  drawback  of  the  above  formulation  is  that  the  disturbance  model  and  the 
plant  model  must  be  known  before  the  optimal  controller  can  be  computed. 

3.3  Multiple-Microphone  Optimal  Feedback  Controller  for 
ANC 

In  this  section,  the  optimal  control  framework  of  Section  3.2  is  extended  to  develop  a  feed¬ 
back  controller  for  attenuating  the  noise  at  the  locations  of  error  microphones.  The  error 
microphones  are  placed  where  noise  cancellation  is  desired,  and  the  canceling  loudspeaker 
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Power  Spectrum 


(a)  Power  spectrum  of  the  propeller  aircraft  noise  at  the  error  microphone  obtained  using  a  feedback 
controller  designed  without  a  weighting  filter. 


Power  Spectrum 


(b)  Power  spectrum  of  the  proper  aircraft  noise  at  the  error  microphone  obtained  with  a  frequency- 
weighted  feedback  controller. 

Figure  3.4:  Results  of  our  computer  simulations. 
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Weighting  Filter 


Figure  3.5:  Frequency  response  of  the  weighting  filter. 
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Power  Spectrum 


(a)  Power  spectrum  of  the  output  of  the  error  microphone  with  and  without  the  optimal  feedback 
controller. 


Closed  Loop  System  Function 


(b)  Closed-loop  transfer  function  T(e^“')  that  is  obtained  by  using  the  optimal  feedback  controller. 


Figure  3.6:  Performance  of  the  hardware  implementation  of  the  optimal  feedback  controller. 
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Power  Spectnitn  of  the  Original  NoLse 


(a)  Power  spectrum  of  the  propeller  aircraft  noise  at  the  error  microphone. 


Power  Spectrum  of  the  Residuals 


(b)  Comparison  of  the  analog  feedback  controller  to  the  optimal  feedback  controller. 


Figure  3.7:  Performance  of  the  hardware  implementation  of  the  optimal  feedback  controller. 
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is  placed  in  the  vicinity  of  these  microphones.  The  acoustic  objective  is  to  attenuate  the 
mean-square  sum  of  the  pressure  fluctuations  at  the  error  microphones.  The  noise  can¬ 
cellation  environment  is  depicted  in  Fi.g  3.8.  Let  us  point  out  that  our  formulation  can 
accommodate  multiple  canceling  loudspeaker;  however,  for  simplicity  we  will  only  consider 
the  single-loudspeaker  case  in  this  section. 


Figure  3.8:  Multiple-microphone  active  noise  cancellation  system. 


The  block  diagram  corresponding  to  the  system  in  Fig.  3.8  is  depicted  in  Fig.  3.9.  The 
plant  model  is  a  single-input /multiple-output  linear  system  with  input  r{t)  and  outputs 
{c,(t)}|_;^.  The  output  of  the  i-th  microphone  mi{t)  is  the  sum  of  the  unwanted  noise  at 
this  microphone  n,(t),  the  canceling  signal  at  this  microphone  Ci(t),  and  the  measurement 
noise  at  this  microphone  Vi{t).  The  unwanted  noise  signals  are  modeled  as  the 

outputs  of  a  multiple-input/multiple-output  linear  system  driven  by  the  white  Gaussian 
vector- valued  process  w{t),  where  w{t)  =  The  objective  is  to  choose  {r(t)  :  t  = 

0,. .  .,N  —  1},  based  on  the  observations  of  such  that 


J{N)  =  E 


t=o  \  \t=l 


I  +  P  ^^(0  )  + 


(3.51) 


is  minimized,  where  is  a  positive  weighting  factor  for  noise  cancellation  at  the  i-th  mi¬ 
crophone,  and  p  is  a  weighting  factor  for  the  overall  control  effort. 


Figure  3.9:  Block  diagram  of  the  multiple- microphone  active  noise  cancellation  system. 
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The  unwanted  noise  signals  at  the  microphones  are  modeled  as  the  I  outputs  of  the 
following  linear  system: 

Xd{t)  =  -  1)  +  Ldw(t  -  1)  (3.52) 

n(t)  =  HdXd{t  -  1),  (3.53) 

where  n(t)  is  the  (/  x  1)  vector  of  the  unwanted  noise  signals 

n{t)  =  [ni(t),  n2(t),  ...  ,n;(t)]^; 

Xd(t)  is  the  state  vector  of  this  system,  and  the  input  of  this  system  is  the  (/  x  1)  vector 
w(t)  with 

E{w(t)w(T)^}  =  S(t  —  r)  J. 

The  plant  model  is  the  following  single-input/multiple-output  linear  system  with  input 
r(t)  and  outputs  {cj(t)}^_i: 


— 

Xp(t  -  1)  +  Gp  r(t  -  1)  +  Lp  u(t  -  1) 

(3.54) 

c(t)  = 

Up  Xp{t  l), 

(3.55) 

where  c(t)  is  the  (/  x  1)  vector  of  the  canceling  signals  at  the  microphones 

c(t)  =  [ci(t),  C2(t),  ...  ,ci(t)f; 

Xp(t)  is  the  state  vector  of  the  plant  model;  and  u(t)  is  the  process  noise  representing  the 
unmodelled  dynamics  in  the  plant  model  with 

E{u(t)u(T)^}  =  6{t  —  t)I. 

Following  the  same  steps  as  the  ones  for  the  single-microphone  case,  we  now  derive  the 
state-space  model  for  the  overall  system  by  combining  (3.52-3.53)  with  (3.54-3.55): 

w(t  —  1) 

x(t)  =  ^x(t-l)  +  Gr(t-l)  +  L  (3.56) 

u(t  —  1) 

Tn(t)  =  ffx(t)  +  v(t),  (3.57) 
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where  x{t)  is  the  state  vector  of  the  overall  system 


x{t)  = 


Xdit) 

Xp[t) 


m(/)  is  the  (/  x  1)  vector  of  the  microphone  outputs 


m{t)  = 


mi{t) 


mi{t) 


$  is  the  transition  matrix  of  the  overall  system 


$  = 


0 

0 


matrix  G  is  defined  as 


G  = 


Gr. 


matrix  L  is  defined  as 


L  = 


Ld  0 


0  Lr. 


and  matrix  H  is  defined  as 


The  (I  X  1)  vector  v{t)  is 


H  =  [Hd  Hp]. 

Vi{t) 


v{t)  = 


Vl{t) 


where  Vi{t)  is  the  measurement  noise  at  the  i-th  microphone  and  is  assumed  to  be  a  white 
Gaussian  process  with  variance  a^.  Furthermore,  we  shall  assume  that  Vi{t)  and  Vj{t)  are 
uncorrelated  for  i  ^  j,  and  v{t)  is  uncorrelated  with  w{t)  and  u{t). 
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Next  we  shall  express  the  performance  index  J{N)  in  terms  of  the  state  vector  of  the 
overall  system.  For  concreteness,  let  e(<)  be  the  vector  of  the  residual  signals 


and  let  A  be  the  (/  x  /)  diagonal  matrix 


Recalling  that 

e{t)  =  H  x{t), 


we  can  rewrite  J{N)  as 


riv-i  >1 

J{N)  =  I  ^  {x{tfH^A  H  x{t)  +  p  +  x{nY H'^ A  Hx{N)  | . 


Since  the  performance  measure  J{N)  is  the  expected  value  of  a  quadratic  form  in  the 
state  and  the  control,  the  minimizing  choice  for  {r{t)  :  t  =  0, . . .,  iV  —  1}  is  again  a  linear 
feedback  of  the  optimal  state  estimates  [24].  The  optimal  state  estimate  can  be  computed 
using  a  Kalman  filter,  and  the  feedback  gain  can  be  computed  by  iterating  an  appropriate 
Riccati  difference  equation.  The  exact  equations  for  the  controller  are  obtained  following 
steps  identical  to  the  ones  followed  in  the  single- microphone  case. 


3.4  Comparison  of  the  Optimal  Feedback  Controller  and 
the  Optimal  Feedforward  Controller 

We  show  in  this  Section  that  the  optimal  feedback  controller  can  be  thought  of  as  a  gen¬ 
eralization  of  the  optimal  feedforward  controller.  To  this  end,  let  us  consider  the  acoustic 
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environment  depicted  in  Fig.  3.10.  We  assume  that  the  desired  acoustic  objective  is  to 
attenuate  the  pressure  at  Microphone  1  and  that  there  is  no  feedback  from  the  input  r(t) 
of  the  canceling  loudspeaker  r(i)  to  the  output  of  Microphone  2.  Furthermore,  we  assume 
that  the  noise  signals  at  the  two  microphones  are  jointly  stationary.  We  will  compare  the 
effect  of  using  a  two-microphone  feedback  controller  to  the  effect  of  using  a  single-channel 
feedforward  controller  for  attenuating  the  acoustic  pressure  at  Microphone  1.  All  signals  in 
this  Section  are  continuous-time  signals. 

For  compatibility  with  the  analysis  done  in  Chapter  1,  we  assume  in  this  Section  that 
the  output  of  each  microphone  is  simply  the  value  of  the  acoustic  pressure  at  the  location 
of  the  microphone,  i.e.  we  ignore  the  microphone  measurement  noise. 

Canceling  Loudspeaker 


n{t) 

2 


Figure  3.10:  Noise  cancellation  environment  with  one  canceling  loudspeaker  and  two  mi¬ 
crophones. 

In  Fig.  3.11,  the  two- microphone  feedback  controller  is  used  to  attenuate  the  noise 
at  Microphone  1.  The  block  diagram  for  this  feedback  system  is  depicted  in  Fig.  3.13, 
where  G(s)  is  the  transfer  function  from  the  input  r(t)  of  the  canceling  loudspeaker  to  the 
canceling  field  at  the  location  of  Microphone  1.  The  output  of  Microphone  1,  €i(t),  is  the 
sum  of  the  noise  at  this  microphone,  ni(t),  and  the  canceling  signal  at  this  microphone 
ci(t).  There  is  no  feedback  from  the  input  of  the  canceling  loudspeaker  to  the  output  of 
Microphone  2;  hence,  the  output  of  Microphone  2  is  simply  the  noise  at  this  microphone 
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Feedback  Controller 


Figure  3.11:  Two- microphone  feedback  controller. 


Canceling  Loudspeaker 


Figure  3.12:  Single- channel  feedforward  controller. 


Figure  3.13:  Block  diagram  of  the  two-microphone  feedback  ANC  system. 


e(l) 
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H 


Figure  3.14:  Restricted  class  of  two-input/single-output  controllers. 
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Q 


Figure  3.15:  Alternate  block  diagram  for  the  two- microphone  feedback  ANC  system. 


n,W 


njt) 

P^O) 

ff 

G 

H 

Figure  3.16:  Block  diagram  for  the  single- channel  feedforward  ANC  system. 

n2{t).  H(^  is  the  transfer  function  from  ei(t)  to  r{t),  and  is  the  transfer  function  from 
n2{t)  to  r{t). 

In  Fig.  3.13,  let  us  find  the  LTI  controller  that  minimizes  E{el{t)}.  Our  strategy 
for  finding  this  controller  is  similar  to  the  strategy  we  used  to  find  the  optimal  single¬ 
channel  feedback  controller  in  Chapter  2,  i.e.  we  first  find  the  minimizing  controller  within 
a  restricted  class  of  controllers,  and  then  show  that  no  other  controller  can  outperform  this 
controller.  We  shall  make  the  additional  assumption  that  G{s)  is  stable.  This  restricted 
clcLss  consists  of  all  two-input /single-output  LTI  systems  of  the  form  depicted  in  Fig.  3.14, 
where  G{s)  (the  same  as  (?(s)  in  Fig.  3.13)  is  fixed,  and  {Q11Q2}  can  be  varied.  For  a 
particular  choice  of  {Qi,Q2}  in  Fig.  3.14,  the  block  diagram  for  the  overall  ANC  system 
can  be  drawn  as  depicted  in  Fig.  3.15. 

Looking  at  Fig.  3.15,  we  see  that  the  minimizing  choice  for  the  canceling  signal  Ci(t) 
is  the  LLSE  of  — ni(f)  based  on  {pi(t),P2(7')  :  t  <  t}-  This  corresponds  to  choosing 
{Q15Q2}  so  that  the  resulting  two-input/single-output  system  Q  is  the  causal  Wiener  filter 
for  estimating  —ni{t)  based  on  {pi('r),P2('r)  '■  t  <t}.  An  argument  similar  to  the  one  given 
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in  Chapter  2  can  be  used  to  show  that  no  other  controller  can  outperform  this  controller 
and  that  the  optimal  canceling  signal  at  Microphone  1,  ci{t)  in  Fig.  3.13,  is  unique. 

In  Fig.  3.12,  a  single- channel  feedforward  controller  is  used  to  attenuate  the  noise  at 
Microphone  1.  The  block  diagram  corresponding  to  this  system  is  depicted  in  Fig.  3.16, 
where  G{s)  is  the  transfer  function  from  the  input  of  the  canceling  loudspeaker  r{t)  to  the 
value  of  the  canceling  field  at  the  location  of  Microphone  1.  The  output  of  Microphone 
1,  ei{t),  is  the  sum  of  the  canceling  signal  at  this  microphone,  ci(t),  and  the  noise  at  this 
microphone  ni{t).  The  reference  signal  for  the  feedforward  controller  is  n2(t),  the  output 
of  Microphone  2.  The  objective  is  to  find  the  LTI  controller  that  minimizes  £'{ef(t)}. 

According  to  the  analysis  presented  in  Chapter  2  for  the  optimal  single- channel  feed¬ 
forward  controller,  the  LTI  system  that  minimizes  E{el(t)}  in  Fig.  3.16  is  the  causal 
Wiener  filter  for  estimating  —ni{t)  based  on  P2it),  where  P2(0  =  ^2(0  *5(0-  The  resulting 
canceling  signal  in  this  case  is  the  LLSE  of  —ni{t)  based  on  {p2(t)  ■  t  <t}. 

Based  on  the  above  analysis,  we  can  now  compare  the  optimal  feedback  controller  to 
the  optimal  feedforward  controller.  With  the  optimal  two-microphone  feedback  controller, 
the  canceling  signal  at  Microphone  1  is  the  LLSE  of  —ni{t)  based  on  {pi{t),P2{t)  :  t  <  t}. 
With  the  optimal  single- channel  feedforward  controller,  the  canceling  signal  at  Microphone 
1  is  the  LLSE  of  — ni(t)  based  on  {p2{t)  <  t}-  Therefore,  with  the  optimal  two- 

microphone  feedback  controller  the  residual  signal  at  Microphone  1  is  the  estimation  error 
in  estimating  ni{t)  based  on  pi{t)  and  P2(0i  however,  with  the  optimal  single-channel  feed¬ 
forward  controller  the  residual  signal  at  Microphone  1  is  the  estimation  error  in  estimating 
ni(t)  based  on  P2(0  alone.  In  this  sense,  the  optimal  feedforward  controller  does  not  use 
aU  the  available  measurements. 

Let  us  note  in  passing  that  the  optimal  single- microphone  feedback  controller  in  this 
case,  using  Microphone  1  only,  will  result  in  a  canceling  signal  at  Microphone  1  that  is  the 
LLSE  of  -ni{t)  based  on  {pi(r) r  <t}  alone. 
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3.5  Summary 


An  optimal  feedback  controller  for  attenuating  the  acoustic  pressure  at  a  finite  number  of 
microphones  was  developed  in  this  chapter.  The  results  in  this  chapter  were  based  on  the 
following  assumptions; 

1.  the  unwanted  noise  signals  at  the  microphones  are  modeled  as  the  outputs  of  a  known 
linear  system  driven  by  a  white  process. 

2.  the  canceling  signals  at  the  microphones  are  modeled  as  the  outputs  of  another  known 
linear  system  driven  by  the  control  signal. 

3.  the  residual  signals  at  the  microphones  are  modeled  as  the  sum  of  the  unwanted  noise 
signals  and  the  canceling  signals. 

4.  the  observed  signals  are  modeled  as  the  sum  of  the  residual  signals  and  white  mea¬ 
surement  noise,  and  this  measurement  noise  is  uncorrelated  with  all  other  processes 
involved. 

Given  a  stochastic  model  for  the  unwanted  noise  at  the  microphones  and  a  linear  sys¬ 
tem  model  relating  the  canceling  loudspeaker  input  of  the  outputs  of  the  microphones,  our 
formulation  can  be  used  to  determine  the  minimum  achievable  value  of  the  statistical  per¬ 
formance  criterion  J  which  was  defined  in  eq.  (3.51).  We  also  showed  that  our  formulation 
can  be  modified  to  minimize  the  frequency  weighted  power  of  the  residual  signal.  More¬ 
over,  our  formulation  can  be  used  to  find  the  specific  feedback  controller  that  achieves  this 
minimum  value. 

The  feedback  controller  developed  in  this  chapter  has  a  few  advantages  compared  to 
the  existing  feedback  controllers.  Recall  that  the  analog  feedback  controller  developed  by 
Wheeler  [51]  does  not  take  advantage  of  the  statistics  of  the  noise.  However,  the  feedback 
controller  developed  in  this  chapter  takes  full  advantage  of  the  statistical  characteristics 
of  the  noise.  The  single-microphone  feedback  controller  developed  by  Graupe  [26]  can  not 
be  used  in  cases  for  which  the  plant  model  contains  one  or  more  delays.  However,  the 
feedback  controller  developed  in  this  chapter  can  be  used  in  these  cases.  Furthermore,  the 
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feedback  controller  developed  in  this  chapter  can  be  used  in  the  single-microphone  as  well 
as  multiple- microphone  configuration. 

A  comparison  of  the  optimal  feedback  controller  with  the  optimal  feedforward  controller 
was  presented  in  Section  3.4.  We  showed  that  both  controllers  generate  the  canceling  signal 
as  the  negative  of  the  LLSE  of  the  unwanted  noise  at  the  error  microphone.  However,  the 
feedback  controller  estimates  the  unwanted  noise  based  on  all  the  available  measurements, 
as  opposed  to  the  feedforward  controller  that  estimates  the  unwanted  noise  based  only  on 
the  reference  signal.  In  this  sense,  the  optimal  feedback  controller  can  be  thought  of  as  a 
generalization  of  the  optimal  feedforward  controller. 

The  main  drawback  of  the  design  procedure  presented  in  this  chapter  is  that  the  dis¬ 
turbance  model  and  the  plant  model  must  be  known  before  the  optimal  controller  can  be 
computed.  The  next  chapter  is  devoted  to  the  modeling  of  the  unwanted  noise  based  on 
the  measurements  of  this  noise  made  by  the  microphones. 
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Chapter  4 


Modeling  of  the  Unwanted  Noise 
at  a  Single  Microphone 


4.1  Introduction 

The  design  procedure  of  Chapter  3  is  only  useful  if  an  algorithm  for  identifying  the  distur¬ 
bance  model  and  the  plant  model  exists.  The  problem  of  identifying  the  disturbance  model 
in  the  single-microphone  case  is  considered  in  this  chapter.  We  assume  that  the  unwanted 
noise  at  the  location  of  the  microphone  is  a  non-stationary  autoregressive  process  and  de¬ 
velop  a  recursive/adaptive  procedure  for  estimating  the  parameters  of  this  process  based 
on  the  measurements  made  by  the  microphone. 

An  adaptive  feedback  ANC  system  can  be  obtained  by  combining  the  design  procedure 
of  Chapter  3  with  an  adaptive  parameter  estimation  algorithm.  Specifically,  an  adaptive 
algorithm  can  be  used  to  continuously  estimate  the  parameters  of  the  disturbance  model 
and  the  parameters  of  the  plant  model,  and  these  parameter  estimates  can  be  used  to 
continuously  adjust  the  controller  according  to  the  design  procedure  developed  in  Chapter 
3.  This  approach  is  commonly  referred  to  as  “indirect  adaptive  control”  [58].  Hence,  an 
optimal  control  algorithm  such  as  that  of  Chapter  3  is  an  essential  part  of  every  indirect 
adaptive  optimal  control  strategy  for  feedback  active  noise  cancellation. 

While  it  is  eventually  necessary  to  fully  develop  the  algorithm  by  adaptively  estimating 
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the  plant  model  and  the  disturbance  model,  we  focus  in  this  chapter  on  the  simpler  problem 
in  which  only  the  disturbance  model  is  estimated.  The  algorithms  presented  in  this  chapter 
can  be  applied  to  the  more  general  problem  of  identifying  non-stationary  autoregressive 
(AR)  processes  embedded  in  white  noise. 

4.2  Model  Specifications 

The  single- microphone  ANC  system  of  interest  and  the  corresponding  block  diagram  are 
depicted  in  Fig.  4.1.  System  G  relates  the  input  r{t)  of  the  canceling  loudspeaker  to  the 
value  of  the  canceling  field  c(t)  at  the  location  of  the  microphone.  The  microphone  output 
m{t)  is  the  sum  of  c{t)  and  the  value  of  the  noise  field  at  the  location  of  microphone  n{t) 
and  the  microphone  measurement  noise  v{t).  The  goal  in  this  chapter  is  to  develop  a 
recursive/ adaptive  algorithm  for  modeling  the  noise  n{t)  based  on  the  measurements  made 
by  the  microphone  m{t). 

The  overall  strategy  is  based  on  the  observation  that  if  G  is  known  exactly,  then  since 
r{t)  is  known  exactly,  an  estimate  of  the  unwanted  noise  at  the  location  of  the  microphone 
can  be  obtained  by  subtracting  out  the  component  of  the  microphone  output  due  to  the 
canceling  source.  Specifically,  if  we  define  z{t)  as 

z{t)  =  n{t)  +  v{t),  (4.1) 

it  is  clear  that  z{t)  can  be  computed  from  the  measurements  of  m,{t)  according  to 

z{t)  =  m{t)  -  G({r(t)}).  (4.2) 

Hence,  if  G  is  known  exactly,  the  problem  of  identifying  the  disturbance  model  based  on 
m{i)  reduces  to  the  problem  of  identifying  the  disturbance  model  based  on  z(t). 
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Loudspeaker 
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4.3  Single-Microphone  Recursive/ Adaptive  Identification  Al¬ 
gorithm 

In  the  reminder  of  this  chapter  we  focus  on  modeling  of  the  unwanted  noise  n{t)  based  on 
the  observations  of 

z{t)  =  n{t)  +  v{t),  (4.3) 

where  the  microphone  measurement  noise  v{i)  is  assumed  to  be  a  zero-mean  white  Gaussian 
process  with  average  power  cr^.  Furthermore,  it  is  assumed  that  n{t)  is  an  autoregressive 
process,  i.e. 

p 

n{t)  = -Y^akn{t  -  k)  +  w{t},  (4.4) 

fc=i 

where  w{t)  is  a  zero- mean  white  Gaussian  process  with  average  power  of  a^,  and  w{t)  is 
independent  of  v{t).  A  more  detailed  discussion  of  the  the  algorithms  developed  in  this 
chapter  can  be  found  in  [55]. 

Using  the  state-space  description  of  equations  (4.3)  and  (4.4),  we  will  develop  two  al¬ 
gorithms  for  estimating  the  parameters  of  the  disturbance  model  (oi, . . . ,  Op  and  cr^).  We 
shall  assume  that  the  average  power  of  the  measurement  noise,  al,  is  known.  Equations 
(4.3)  and  (4.4)  can  be  expressed  in  state-space  form  as  [34]: 

a:(t)  =  ^dx{t  -  1)  -f  Ldw{t)  (4.5) 

z(t)  =  Hdx(t)  +  v{t)  (4.6) 

where  x{t)  is  the  (p  1)  x  1  state  vector  defined  by: 

x{t)  =  [n(t),  n{t  -  1),  . . . ,  n{t-  p)]^;  (4.7) 
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the  optimal  state  estimate  based  on  data  up  to  time  “t”,  and  by 

P(t\t)  =  E{[x{t\t)  -  a;(0][®(t|t)  -  ®(t)]"^|2(l),...,2(t)}  (4.15) 

the  associated  error  covariance  matrix.  In  accordance  with  (4.13),  an  estimate  of  R  can  be 
computed  by  replacing  the  expectation  operator  in  eq.  (4.13)  with  the  weighted  averaging; 

X!  A‘"^a:(r)a:^(r)  (4.16) 

^T=l  ^  T  =  1 

where 

i(t|t)i^(t|f)  +  P(t|t)  (4-17) 

E{x{t)x^  {t)\z{l), .  (4.18) 

Then  using  the  standard  Kalman  filter  formulation  x{t\t)  and  P(tlt)  in  eq.  (4.17)  can 


be  computed  sequentially  in  time,  in  two  stages,  as  follows^: 

Propagation  Equations: 

i(t|t  —  1)  =  $da:(t  —  l|t  —  1)  (4-19) 

P(t|t-1)  =  ^dP{t  -  l\t  -  inl  +  alLdll  (4.20) 

Updating  Equations: 

x{t\t)  =  x{t\t  —  1)  +  k{t)[z{t)  —  H  dx{t\t  —  1)]  (4-21) 

P(t|t)  =  [I-kit)Hd]Pit\t-l),  (4.22) 

where  k(t)  is  the  Kalman  gain  given  by 

k{t)  =  — - - - P{t\t  -  l)Hd.  (4.23) 

H]Pit\t-l)Hd  +  (T^ 


We  would  like  to  generate  the  parameters  of  the  disturbance  model  in  such  a  way  that 
equation  (4.11)  is  satisfied;  however,  this  requires  knowing  matrix  R,  which  in  turn  requires 

^These  equations  can  be  simplified  by  exploiting  the  structure  of  Ld,  and  Hd,  ss  was  done  in  [56]; 
however,  this  is  not  essential  to  our  development  here. 


93 


knowing  the  exact  parameters  of  the  disturbance  model.  Instead,  we  propose  generating 
these  parameter  estimates  in  such  a  way  that  the  following  equation  is  satisfied; 


+  1) 


(4.24) 


+  1) 


where  R{t)  is  the  estimate  of  R,  obtained  by  performing  the  weighted  averaging; 


where 


2^r=l  T=1 


x{i)x  (t)  =  x(t\t)x  (t|t)  +  P(t|t), 


(4.25) 


(4.26) 


where  x{t\t)  and  P(t|t)  are  the  estimate  of  the  state  and  its  covariance,  computed  using 
the  Kalman  filtering  equations  (4.19)-(4.23),  where  instead  of  a.  and  we  use  the  most 
current  estimates  ot{t)  and  o'^(t),  i.e. 


Propagation  Equations: 


a;(t|t  —  1)  =  ^  d{t)x{t  —  l\t  —  1) 


Updating  Equations: 


a:(t|t)  =  x{t\t  —  1)  +  k{t)[z{t)  —  H  ^x{t\t  —  1)] 
P{t\t)  =  [I-k{t)H]]P{t\t-l), 


(4.27) 

(4.28) 


(4.29) 

(4.30) 


where  $d(t)  is  the  matrix  defined  in  eq.  (4.8)  computed  at  a  =  oc(t),  and  k{t)  is  the  vector 
defined  in  eq.  (4.23)  with  P{t\t  —  1)  obtained  from  eq.  (4.28). 

Our  goal  now  is  to  develop  a  sequential  procedure  for  generating  the  parameter  estimates 
in  accordance  with  equation  (4.24),  i.e.  the  parameter  estimates  at  time  t  +  1,  a{t  +  1)  and 
cr^{t  +  1),  must  satisfy  equation  (4.24).  To  obtain  this  sequential  procedure  for  updating 
the  parameter  estimates,  we  define 


^  1  Qii(i)  Qi2(0 

P  I  Q2i(0  ^22(0 


=  Q(0  =  ^®(T-)a:^(r) 


(4.31) 


1 


p 
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(4.32) 


=  x(t)x  (t)  +  AQ(/  —  1). 


Using  eq.  (4.31)  and  eq.  (4.25),  we  express  eq.  (4.24)  as 

1 


1 

<5ii(0 

^12(0 

EUl 

Q2i(0 

^22(0 

1) 


(4.33) 


d(t+l) 

which  leads  to  the  following  update  equations  for  the  noise  parameter  estimates: 


a(f  +  l)  =  _Q-'(t)Q2i(t)  (4.34) 

+  1)  =  — [(5ii(t)  +  Q^2(^)d(t  +  1)],  (4.35) 

where  we  have  used  the  fact  that 
the  non-gradient  algorithm. 

At  each  time  step  the  non-gradient  algorithm  first  estimates  the  current  state  using 
the  latest  parameter  estimates  in  eqs.  (4.27)  -  (4.30),  and  then  updates  the  parameter 
estimates  using  the  state  estimate  just  computed  and  its  covariance  in  eqs.  (4.34)-  (4.35). 
The  intermediate  quantity  Q(t)  is  computed  recursively  in  time  according  to  eq.  (4.32). 

The  factor  A  that  appears  in  the  cumulative  averaging  in  eq.  (4.31)  is  a  number  between 
0  and  1 .  To  maximize  statistical  stability  we  choose  A  =  1 .  Choosing  A  to  be  strictly  smaller 
than  1  corresponds  to  exponential  weighting  that  gives  more  weight  to  current  data  samples 
and  results  in  an  adaptive  algorithm  that  is  capable  of  tracking  non-stationary  changes  in 
the  structure  of  the  data. 

As  an  alternative  to  the  parameter  update  equations  (4.31)-  (4.35),  we  may  consider  a 
gradient-search  algorithm  for  solving  the  Yule- Walker  equation  (4.11).  In  this  case,  instead 
of  the  signal  correlation  matrix  R  we  use  its  estimate  given  by  eq.  (4.25),  and  we  proceed 
sequentially  through  the  data  using  a  stochastic  gradient  type  procedure  [43] .  Replacing  B. 
in  eq.  (4.11)  by  its  current  estimate  R(t)  =  yf^Q(t)  where  Q(t)  is  defined  in  eq.  (4.31), 
we  obtain: 


=  jzy-  We  shall  refer  to  eqs.  (4.34)-  (4.35)  as 


^2i(0  d-  ^22(0*^  ~  ° 


(4.36) 


cr,„  — 


1  -  A 
1  -  A< 


[Qii(0  +  Qi2(0<^]  —  0- 


(4.37) 
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From  eqs.  (4.36)  and  (4.37)  and  using  an  approach  similar  to  that  used  in  compound  decision 
problems  (e.g.,  see  [25,  38,  42,  45,  54]),  the  following  sequential  update  equations  for  the 
signal  parameters  are  suggested: 


+  1) 
+  1) 


«(0  -  7t[<32l(0  +  <?22Wa(0] 

‘522(0“(0] 


-  It 


(4.38) 

(4.39) 


where  Qij{t)  i,j  =1,2  are  computed  recursively  in  t  using  eq.  (4.32),  and  7<  is  the  step-sizes 
of  the  stochastic  gradient  algorithm.  We  shall  refer  to  eqs.  (4.38)-(4.39)  as  the  gradient 
algorithm. 

The  advantage  of  using  the  gradient  algorithm  specified  by  eq.  (4.38)  -  (4.39)  is  that 
it  does  not  require  matrix  inversion  in  contrast  to  the  non-gradient  algorithm  specified  by 
eqs.  (4.34)-(4.35),  and  therefore  the  gradient  algorithm  is  computationally  simpler  than  the 
non- gradient  algorithm. 


4.4  Algorithm  Performance 

The  results  of  applying  the  algorithms  developed  in  the  previous  section  to  synthetically 
generated  noise  and  to  recorded  aircraft  noise  are  presented  in  this  section. 

4.4.1  Synthetic  Noise 

Two  sets  of  experiments  were  performed  in  which  the  gradient  algorithm  of  the  previous 
section  was  applied  to  computer  generated  autoregressive  time  series.  The  first  set  of 
experiments  was  performed  to  show  how  close  the  estimated  spectrum  of  the  AR  process 
matches  the  actual  spectrum.  The  second  set  of  experiments  was  performed  to  calculate 
the  one-step-ahead  prediction  error  that  results  from  using  the  estimated  parameters.  This 
prediction  error  was  then  compared  to  the  one-step-ahead  prediction  error  that  results  from 
using  the  correct  parameters.  In  all  these  experiments  the  model  order  and  a-u  were  assumed 
to  be  known  exactly.  The  variance  of  the  AR  process  {E{n^{t)}  in  (4.3))  was  equal  to  1  in 
these  experiments,  and  the  forgetting  factor  A  was  set  equal  to  0.975. 
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The  first  set  of  experiments  presented  in  this  section  illustrates  how  accurately  the 
algorithms  developed  in  this  chapter  can  estimate  the  model  for  the  unwanted  noise  that  is 
needed  to  compute  the  optimal  controller  of  Chapter  3.  The  design  procedure  of  Chapter 
3  relies  solely  on  the  second-order  statistics  of  the  unwanted  noise  and  the  measurement 
noise.  Since  the  measurement  noise  is  assumed  to  be  white,  we  can  equivalently  say  that 
the  design  procedure  of  Chapter  3  relies  solely  on  the  power  spectrum  of  the  unwanted  noise 
and  the  variance  of  the  measurement  noise.  Therefore,  it  is  important  to  see  how  close  the 
estimated  power  spectrum  of  the  unwanted  noise  matches  the  actual  power  spectrum  of  the 
unwanted  noise. 

In  Fig.  4.2(a)  and  Fig.  4.3(a),  the  doted  curve  is  the  spectrum  of  the  actual  underlying 
AR  process,  and  the  solid  curve  is  the  average  estimated  spectrum  obtained  by  using  the 
gradient  algorithm.  The  solid  curve  is  computed  by  averaging  the  estimated  spectrum  over 
6000  consecutive  points.  The  actual  system  parameters  for  these  two  figures  are: 

Fig. I :  A{z)  =  53^-\+o  =  0.234  =  0.05 

1  1.53Z  +0.6Z 

Fig.2  :  A(z)  =  =  0-0501  =  0.05. 

In  other  words,  the  AR  process  was  generated  as  the  output  of  A(z)  with  a  white  noise 
input  of  variance  cr^. 

In  Fig.  4.2(b)  and  Fig.  4.3(b),  the  dotted  horizontal  line  is  the  average  one-step-ahead 
prediction  error  of  a  Kalman  predictor  that  uses  the  correct  system  parameters.  This  is 
the  optimum  predictor,  in  the  sense  that  no  other  predictor  can  achieve  a  lower  average 
prediction  error.  The  solid  curve  in  these  two  figures  is  the  average  prediction  error  obtained 
by  using  the  gradient  algorithm. 

4.4.2  Aircraft  Noise 

The  results  of  applying  both  the  gradient  and  the  non-gradient  algorithm  to  three  types  of 
aircraft  noise  are  presented  in  this  section.  These  two  algorithms  were  used  on  the  noise 
generated  by  a  propeller  aircraft,  a  helicopter,  and  a  jet  aircraft.  Recordings  of  these  three 
types  of  noise  were  made  from  a  microphone  placed  inside  a  set  of  headphones  approximately 
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two  meters  away  from  each  aircraft.  Note  that  these  recordings  correspond  to  z{t)  in  eq. 
(4.3). 

These  experiments  simulate  an  indirect  adaptive  feedback  controller  in  which  the  plant 
is  a  known  delay  of  M  samples,  i.e.  G{z)  =  z~^ .  In  Fig.  4.1,  with  G  —  z~^ ,  the  feedback 
controller  that  minimizes  jF{(n(t)  +  c(t))^}  is  the  M-step-ahead  predictor  that  generates  the 
LLSE  of  -n(t)  based  on  {m(r)  :  t  <  t  -  M},  i.e.  the  optimal  canceling  signal  c{t)  in  this 
case  is  the  LLSE  of  —n{t)  based  on  {m(r)  :  t  <  t  —  M}.  The  resulting  residual  signal  e(t) 
in  this  case  is  the  corresponding  M-step  ahead  prediction  error,  i.e. 

e(t)  =  n{t)  -  E{n{t)\m{T) :  t  <t  -  M}.  (4-41) 

The  performance  of  the  adaptive  feedback  controller  in  this  case  depends  on  how  well  n{t) 
can  be  predicted  based  on  {m(r)  :  t  <t  —  M},  or  equivalently  based  on  {z{t)  :t  <t  —  M}. 
In  short  the  residual  signal  in  this  case  is  simply  the  prediction  error  in  predicting  n{t) 
based  on  {z{t)  :  t  <  t  —  M}  and  based  on  the  estimated  parameters  of  the  noise  model. 

Several  simulations  were  performed  to  illustrate  how  well  n{t)  can  be  predicted  based 
on  measurements  of  {z(t)  :  t  <  t  —  M}  and  based  on  the  estimated  parameters  of  the 
model  of  the  noise  obtained  from  our  gradient  algorithm.  In  these  simulations,  the  order 
of  the  AR  model  for  the  unwanted  noise  was  assumed  to  be  five,  and  Cy  was  fixed  at  5% 
of  the  standard  deviation  of  z{t).  The  algorithms  were  tried  with  the  order  of  the  AR 
model  ranging  from  three  to  nine,  and  it  was  found  that  the  performance  of  the  algorithm 
improves  very  little  by  increasing  the  order  beyond  five.  The  variance  of  the  microphone 
measurement  noise,  can  be  estimated  in  the  following  way.  The  microphone  is  placed 
in  a  quiet  environment  with  the  ANC  turned  off  so  that  the  output  of  the  microphone  is 
approximately  v{t).  An  estimate  of  can  then  be  obtained  as 

(4-42) 

^  t=i 

Figure  4.4  is  a  plot  of  the  normalized  prediction  error  versus  the  prediction  time  for  the 
non-gradient  algorithm.  Similarly,  Figure  4.5  is  the  plot  of  the  normalized  prediction  error 
versus  the  prediction  time  for  the  gradient  algorithm.  In  these  two  figures,  the  normalized 
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prediction  error  is  calculated  as 


Normalized  Prediction  Error(dB)  =  lOlog^o 


E{z^it)r 


(4.43) 


where  E{z'^{t)}  is  the  average  power  of  the  original  recorded  noise,  and  E{e'^(t)}  is  the 
average  power  of  the  prediction  error. 

An  intuitive  explanation  for  why  the  normalized  prediction  error  of  the  propeller  noise 
is  less  than  the  normalized  prediction  error  of  the  helicopter  noise  and  the  normalized 
prediction  error  of  the  jet  noise.  The  propeller  noise  is  more  narrow-band  than  the  helicopter 
noise  and  the  jet  noise;  therefore,  the  propeller  noise  results  in  the  least  prediction  error. 

Furthermore,  it  is  evident  from  these  set  of  experiments  that  the  performance  of  the 
gradient  algorithm  is  slightly  better  than  the  non-gradient  algorithm.  Furthermore,  the 
gradient  algorithm  is  much  less  computationally  intensive  than  the  non-gradient  algorithm. 


4.5  Summary 

Two  recursive/ adaptive  algorithms  for  modeling  the  unwanted  noise  at  a  single  microphone 
were  presented  in  this  chapter.  We  argued  that  if  the  plant  model  is  exactly  known,  these 
algorithms  can  be  combined  with  the  design  procedure  of  Chapter  3  to  obtain  an  adaptive 
feedback  ANC  system.  Of  course,  it  is  eventually  necessary  to  fully  develop  these  algorithms 
so  that  the  disturbance  model  and  the  plant  model  are  both  adaptively  estimated. 

We  assumed  that  the  unwanted  noise  was  a  non-stationary  autoregressive  process,  and 
that  the  observed  signal  was  the  sum  of  this  autoregressive  process  plus  white  measurement 
noise.  We  then  developed  two  algorithms  for  estimating  the  parameters  of  this  autoregres¬ 
sive  process  based  on  the  observed  data.  The  performance  of  these  algorithms  were  studied 
by  applying  them  to  synthetically  generated  noise  and  recorded  aircraft  noise.  It  is  worth 
noting  that  these  algorithms  can  be  applied  to  the  more  general  problem  of  identifying 
non-stationary  autoregressive  processes  embedded  in  white  noise. 
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Ideal  and  Estimated  Spectrum 


Ideal  and  Realized  Prediction  Error 


(b) 


Figure  4.2:  Performance  of  the  gradient  algorithm  on  synthetic  noise 


Ideal  and  Estimated  Spectrum 
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Prediction  Time  in  Seconds 


Figure  4.4;  Performance  of  the  non-gradient  algorithm  on  recorded  aircraft  noise. 
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Prediction  Error 


Prediction  Time  in  |i  Seconds 


Figure  4.5:  Performance  of  the  gradient  algorithm  on  recorded  aircraft  noise. 
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Chapter  5 


Distributed  Optimal  Feedback 
Control  for  ANC 

5.1  Introduction 

Our  treatment  of  the  active  noise  cancellation  (ANC)  problem  up  to  this  point  has  been 
limited  to  pointwise  active  noise  cancellation.  In  Chapter  3,  we  developed  an  optimal 
feedback  controller  for  attenuating  the  acoustic  pressure  at  the  locations  of  a  finite  number  of 
error  microphones.  Specifically,  this  controller  minimizes  a  performance  criterion  involving 
the  average  power  of  the  residual  pressure  field  at  the  location  of  the  error  microphones 
and  the  average  power  of  the  control  signal.  Hence,  this  pointwise  criterion  depends  only 
on  the  interaction  between  the  primary  field  and  the  canceling  field  at  those  points  where 
the  microphones  are  located. 

The  distributed  active, noise  cancellation  problem  is  studied  in  this  chapter.  Specifically, 
we  develop  an  optimal  feedback  controller  for  minimizing  a  performance  criterion  involving 
the  total  acoustic  energy  in  an  enclosure  and  the  power  expended  by  the  control  signal.  To 
minimize  this  distributed  performance  criterion,  the  interaction  between  the  primary  field 
and  the  canceling  field  at  every  point  in  the  enclosure  needs  to  be  considered. 

Our  approach  is  to  formulate  the  distributed  ANC  problem  as  an  infinite  dimensional 
feedback  control  problem  and  use  existing  control  theory  to  solve  it.  The  deterministic  for- 
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mulation,  which  assumes  that  the  entire  state  is  observable,  and  the  stochastic  formulation, 
which  estimates  the  state  based  on  the  available  observations,  are  both  considered.  The 
main  assumption  is  that  the  boundary  surrounding  the  volume  of  interest  is  either  purely 
reflecting  or  purely  absorbing. 

Our  formulation  of  the  distributed  ANC  problem  in  this  chapter  is  remarkably  similar 
to  our  formulation  of  the  pointwise  ANC  problem  in  Chapter  3.  In  Chapter  3,  the  unwanted 
noise  signal  was  modeled  as  the  output  of  a  linear  system  driven  by  white  noise,  and  the 
canceling  signal  was  modeled  as  the  output  of  another  linear  system  driven  by  the  control 
signal.  In  this  chapter,  the  unwanted  noise  field  is  modeled  as  the  output  of  a  linear  system 
driven  by  white  noise,  and  the  canceling  field  is  modeled  as  the  output  of  yet  another  linear 
system  driven  by  the  control  signal.  In  Chapter  3,  the  Unear  systems  were  finite  dimensional 
and  were  characterized  by  differential  equations.  In  this  chapter,  the  linear  systems  are 
infinite  dimensional  and  are  characterized  by  partial  differential  equations.  The  optimal 
control  in  either  case  is  a  linear  feedback  of  the  state  estimates.  Furthermore,  in  both  cases, 
the  state  estimates  are  computed  using  a  Kalman  filter,  and  the  feedback  gain  is  obtained 
from  a  Riccati  equation. 

Section  5.2  is  essentially  background  material  and  terminology  that  is  needed  to  under¬ 
stand  our  formulation  of  the  distributed  active  noise  cancellation  problem  in  the  determin¬ 
istic  case.  The  deterministic  quadratic  regulator  problem  is  presented  in  Section  5.3.  The 
acoustic  system  that  we  focus  on  and  its  governing  equations  are  presented  in  Section  5.4. 
In  Section  5.5,  we  derive  an  abstract  differential  equation  representation  for  this  acoustic 
system  and  show  that  the  noise  cancellation  problem  for  this  system  is  equivalent  to  a 
deterministic  quadratic  regulator  problem. 

The  stochastic  formulation  of  the  distributed  active  noise  cancellation  problem  is  pre¬ 
sented  in  Section  5.6.  The  necessary  measure  theoretic  structure  is  introduced  in  Section 
5.6.1,  followed  by  the  stochastic  quadratic  regulator  problem  in  Section  5.6.2.  In  Sec¬ 
tion  5.6.3,  the  noise  cancellation  system  that  we  focus  on  is  modeled  as  a  linear  infinite¬ 
dimensional  system  driven  by  a  white  noise  process  and  a  control  process.  Furthermore, 
we  show  that  this  model  can  accommodate  an  observation  process  that  is  the  output  of  a 
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pressure  microphone,  and  a  control  process  that  is  the  input  to  an  ordinary  loudspeaker. 
Finally,  we  show  that  the  distributed  active  noise  cancellation  problem  for  the  acoustic 
system  that  we  focus  on  is  equivalent  to  a  stochastic  quadratic  regulator  problem. 

5.2  Preliminaries 

This  Section  contains  the  terminology  and  background  from  functional  analysis  that  is 
needed  to  understand  our  particular  formulation  of  infinite-dimensional  linear  systems. 
General  properties  of  semigroups  are  presented  in  Section  5.2.1,  and  abstract  diflFerential 
equations  are  presented  in  Section  5.2.2.  Proofs  for  almost  all  the  theorems  in  this  Section 
can  be  found  in  [5]. 

5.2.1  General  Properties  of  Semigroups 

Let  T{t),  t  >  0,  be  a  family  of  linear  bounded  transformations  mapping  a  Hilbert  space  H 
into  itself.  T{t)  is  said  to  be  a  “semigroup”  if 

(i)  T(0)  =  Identity 

(ii)  T(ti  +  t2)  =  T{t2)T{h)  =  Tih)T{t2). 

The  semigroup  T{t)  is  said  to  be  strongly  continuous  if  for  each  x  in  the  Hilbert  space  H 

(iii)  ||r(t)a:  —  xH  — >■  0  as  t  O’*". 

Using  (i)-(iii),  one  can  show  that  T{t)  is  bounded  on  bounded  domains;  that  is,  for  each 
i  <  00,  there  exists  an  M  such  that 

Sup  l|T(t)||  <  M  <  oo. 

0<t<L 

Moreover,  we  can  find  a  dense  subspace  of  H  on  which  T{t)  is  differentiable.  That  is,  we 
can  define  the  operator  A  (called  the  infinitesimal  generator  of  T{t))  as 

4;,  =  li„,  IXA^ 

A— 0+  A 

X>(A)  =  |x  €  .ff  I  T(A)^ — X  ^ 
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where  X>(A)  is  the  domain  of  A.  Note  that  operator  A  is  a  mapping  from  V{A)  into  H . 
From  the  above  definition  of  A  it  follows  that  T{t)x  is  differentiable  on  VA  for  all  t  >  0 
and 

^iT{t)x)  =  T{t)Ax  =  A(T(t)x)  \/x  €  V(A). 
at 

This  definition  of  A  also  implies  that  A  is  a  closed  linear  operator  and  its  domain  is  dense 
in  H. 

The  relationship  between  the  semigroup  T{t)  and  its  infinitesimal  generator  A  is  similar 
to  the  relationship  between  the  matrix  exponential  and  the  matrix  A  where 

^e^^x  =  Ae^^x. 
dt 

It  will  be  seen  later  that  this  analogy  can  be  carried  much  further.  For  example,  e'^*a;(0)  is  a 
solution  to  the  system  of  ordinary  differential  equations  ^x{t)  =  Ax{t).  Similarly,  T(t)x{0) 
is  a  solution  to  the  abstract  differential  equation  ^x{t)  —  Ax{t),  provided  that  A  is  the 
infinitesimal  generator  of  the  strongly  continuous  semigroup  T{t). 

We  now  state  two  theorems  regarding  the  necessary  and  sufficient  conditions  for  a  dissi¬ 
pative  linear  operator  to  be  the  infinitesimal  generator  of  the  strongly  continuous  semigroup 
T{t).  Recall  that  a  closed  linear  operator  A  mapping  a  Hilbert  space  H  into  itself  and  with 
dense  domain  'D{A)  in  H  is  said  to  be  dissipative  if 

[Aa:,x] -f  [x.  Ax]  <  0  Vx  €  P(A).  (5A) 

We  already  have  seen  that  the  infinitesimal  generator  of  a  strongly  continuous  semigroup 
must  necessarily  be  closed  with  dense  domain  in  Ff.  A  closed  linear  operator  that  maps  a 
Hilbert  space  into  itself  and  has  a  dense  domain  will  have  a  unique  adjoint  operator  whose 
domain  is  also  dense  [33]. 

Theorem  1  Let  A  be  a  closed  linear  operator  with  dense  domain  in  a  Hilbert  space 
H,  and  assume  that  A  maps  its  domain  into  H .  Suppose  further  that  both  A  and  A*,  the 
adjoint  of  A,  are  dissipative.  Then  A  generates  a  strongly  continuous  semigroup  T{t)  such 
that  ||r(t)||  <  1  for  aU  t  >  0. 
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Theorem  2  Suppose  A  is  a  closed  linear  operator  with  dense  domain  V{A)  in  a  Hilbert 
space  H .  Suppose  further  that  A  maps  its  domain  into  H,  and 

0  =  [Ax,x]A[x,Ax\  Vx  €  P(A) 

0  =  [A*x,x]a[x,A*x]^x  £V{A*), 

where  [., .]  is  the  inner  product  in  H.  Then,  A  generates  a  strongly  continuous  semigroup 
T{t)  such  that  ||T(t)x||  =  ||a;||  Vx  €  H. 

From  the  last  two  theorems,  we  see  that  the  semigroup  generated  by  a  dissipative  linear 
operator,  call  it  T{t),  is  a  contraction,  i.e. 

||r(t>||<||a:i|  Vxgf. 

A  contraction  semigroup  wiU  obviously  have  the  property  that  ||r(t2)a;ll  <  lir(ti)x||  for  all 
<  t2-  This  property  is  analogous  to  the  following  property  for  a  matrix  A:  if  [Ax,x]  + 
[a:.  Ax]  <  0  for  all  vectors  x,  then  ||e^‘2a:||  <  for  all  <  t2. 

5.2.2  Abstract  Differential  Equations:  Cauchy  Problem 

The  initial  value  (Cauchy)  problem  for  first  order  linear  abstract  differential  equations  is 
introduced  in  this  Section.  With  the  aid  of  semigroup  theory,  these  abstract  differential 
equations  are  put  in  a  framework  very  similar  to  that  of  ordinary  differential  equations.  A 
specific  example  is  given  to  illustrate  how  a  partial  differential  equation  can  be  viewed  as 
an  abstract  differential  equation. 

Let  A  be  the  infinitesimal  generator  of  a  strongly  continuous  semigroup  T{t)  in  a  Hilbert 
space  H .  The  initial  value  problem  for  the  abstract  differential  equation: 

dxit^ 

=  Ax{i)  x(0)  is  given  and  is  inX>(A)  (5-2) 

is  to  find  the  function  x{t)  with  range  in  Hilbert  space  H  that  satisfies  (5.2).  Then,  x{t)  = 
T{t)x{0)  gives  one  solution,  since  ^{T{t)x(0))  =  A(r(t)a;(0))  and  T(0)ar(0)  =  a;(0).  This  is 
the  only  solution  under  the  additional  assumptions  that: 

1.  x{t)  e  X>(A)  t  >  0 
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2.  a;(^)  is  absolutely  continuous  for  t  >  0 

3.  ||a:(t)  -  a;(0)l|  — 0  as  t  — »■  O"*'. 

From  the  above  discussion,  we  see  that  in  the  case  of  abstract  differential  equations  T{t) 
plays  a  role  analogous  to  the  state  transition  matrix  for  ordinary  differential  equations. 
Next  we  consider  the  nonhomogeneous  equation 

=  Ax{t)  +  u{t)  0  <  t  <  T  (5.3) 

with  a  given  initial  condition  ar(0).  By  analogy  with  the  finite  dimensional  case,  we  would 
expect  that 

x{t)  =  T{t)x{0)  +  f  T{t  —  s)u{s)ds 
Jo 

should  be  the  unique  solution  in  an  appropriate  sense.  Depending  on  the  smoothness  of 
u(t)  and  whether  x(0)  is  in  T>{A)  or  not,  the  sense  in  which  (5.3)  is  satisfied  will  vary.  The 
next  two  theorems  are  concerned  with  this  issue. 

Theorem  3  Suppose  a;(0)  is  in  T>{A),  and  the  function  u{t)  with  range  in  Hilbert 
space  H  is  strongly  continuously  differentiable  in  the  open  interval  (0,  T)  with  derivatives 
continuous  in  the  closed  interval  [0,7].  Then 

^x(t)  =  Ax{t)  +  u{t)  0  <  t  <  T 

has  a  unique  solution  satisfying 

1.  x(t)  €  V{A)  t  >  0 

2.  x{t)  is  absolutely  continuous  in  (0,T) 

3.  \\x{t)  —  x(0)||  0  as  t  — s-  O"*", 

and  this  solution  is  given  by 

x{t)  —  T{t)x{0)  +  /  T{t  —  s)u{s)ds. 

Jo 

If  we  wish  to  consider  the  case  where  neither  the  condition  that  x(0)  belongs  to  the 
domain  of  A,  nor  the  condition  of  smoothness  on  the  input  u(t)  hold,  we  have  to  change 
the  sense  in  which  equation  (5.3)  holds. 


109 


Theorem  4  Suppose  u(.)  €  £2((0,  T);  if),  then  there  is  one  and  only  one  function 
<  t  <T,  such  that  [x(t),y]  is  absolutely  continuous  for  each  y  in  T>(A*),  and 

^[x(t),  y]  =  [a;(t),  A*y]  +  [u(t),  y]  a.e.  0  <  t  <  T,  (5.4) 

and  such  that  for  a  given  x(0)  in  if, 

[a;(t),y]  =  [ar(0),y]  Vy  €  X>(A*). 

Moreover,  this  function  is  given  by 

x{t)  =  r(t)a;(0)  +  f  T{t  —  s)u{s)ds. 

Jo 

The  following  example  illustrates  how  a  partial  differential  equation  can  be  viewed  as 
an  abstract  differential  equation  in  an  appropriate  function  space. 

Example:  Application  to  Partial  Differential  Equations 

Consider  the  partial  differential  equation  for  one- dimensional  heat  transfer 

d 

—f{t,x)=-^f{t,x)  t>0]-oo<x<oo  ;/(0,a;)  is  given. 

Let  us  show  that  this  partial  differential  equation  can  be  represented  as  an  abstract  differ¬ 
ential  equation  of  the  type  discussed  in  the  previous  Section.  For  this,  we  need  to  show 
that  ^  is  an  infinitesimal  generator  of  a  strongly  continuous  semigroup.  Note  that  it  is 
not  necessary  to  explicitly  find  the  semigroup  generated  by  all  that  is  needed  to  is  to 
show  that  such  a  semigroup  exists. 

First  we  choose  the  underlying  Hilbert  space  to  be  if  =  jC2{—oo,  oo),  the  space  of  square 
integrable  functions  with  respect  to  the  Lebesgue  measure.  Next  we  take  A  =  ^  as  our 
linear  operator  with  the  following  dense  domain 

T>(A)  =  {f{x)  e  if  I/, /'are  absolutely  continuous  and  €  if}  . 

This  means  in  particular  that  /(— oo)  =  /(oo)  =  /'(oo)  =  /'(-oo)  =  0.  Domain  of  A  is 
dense  in  if,  since  the  C°°,  the  set  of  all  infinitely  differentiable  functions  in  X2(-oo,+oo), 
is  included  in  the  domain  of  A,  and  C°°  is  dense  in  H  as  shown  in  [1]. 
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Again  using  integration  by  parts,  we  see  that  A  is  dissipative  since 


dx  <  0. 


It  follows  that  operator  A  is  closed,  since  it  is  self  adjoint  and  A*  is  closed  by  the 
virtue  of  being  an  adjoint  operator  [33] .  A  is  a  closed  linear  operator  with  dense  domain 
in  H,  and  both  A  and  A*  are  dissipative;  therefore,  by  theorem  1,  A  generates  a  strongly 
continuous  contraction  semigroup  T{t).  The  original  partial  differential  equation  can  now 
be  represented  as  an  abstract  differential  equation: 


— /(t)  =  Af{t)  t  >  0;  /(O)  G  P(A)  is  given. 


where 


A 

P(A) 


dx'^ 

{/(x)  G  iT]/, /'are  absolutely  continuous  and  /', /"  £  H} . 


Note  that  f{t)  =  T{t)f{0),  where  T{t)  is  the  semigroup  generated  by  A,  is  a  solution  to 
(5.6)  since 

^(T(t)/(0))  =  A(nt)/(0)), 

and  r(/)/(0)  is  in  the  domain  of  A  for  all  values  of  t  with  r(0)/(0)  =  /(O).  In  summary,  we 
have  shown  that  the  operator  A  =  ^  is  the  infinitesimal  generator  of  a  strongly  continuous 
semigroup  T{t),  and  hence  been  able  to  represent  the  original  partial  differential  equation 
as  an  abstract  differential  equation  in  the  Hilbert  space  +oo)- 

Having  discussed  abstract  differential  equations,  we  will  present  the  control  problem  for 
systems  governed  by  these  equations  in  the  next  Section. 
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5.3  Deterministic  Linear  Quadratic  Regulator  Problem 


The  linear  quadratic  regulator  problem  for  systems  governed  by  abstract  differential  equa¬ 
tions  is  presented  in  this  Section.  We  will  see  that  under  very  mild  conditions,  a  unique 
feedback  solution  to  this  problem  exists. 

Let  A  denote  the  infinitesimal  generator  of  a  strongly  continuous  semigroup  of  operators 
T{t)  over  a  Hilbert  space  H.  The  dynamic  equation  governing  the  problem  is 

J 

—x{t)  =  Ax{t)  -i-  Bu{t)  0  <  t  <  T  <  oo  (5.7) 

at 

a;(0)  given, 

where  H  is  a  linear  bounded  transformation  mapping  a  Hilbert  space  into  H,  and  u(.) 
is  an  element  of  Wu  -  ^2((0,  T);  H^).  It  follows  that  Bu{.)  is  an  element  of  >C2((0,  T);  H). 
Hence,  a  unique  solution  to  (5.7),  in  the  sense  of  (5.4),  exists  and  is  given  by 

x{i)  =  r(t)a;(0)  -1-  /  T{t  —  s)Bu{s)ds. 

Jo 

The  optimal  control  problem  is  to  find  u{t)  in  Wu  so  that  the  cost  functional 

J  =  [  [Rx{t),x{t)]Hdt+  f  [u{t),u{t)]Hudt  (5.8) 

Jo  Jo 

is  minimized,  where  J?  is  a  linear  bounded  nonnegative  definite  operator  mapping  H  into 
H .  The  first  integral  in  this  functional  is  essentially  what  we  want  to  minimize,  and  the 
second  integral  is  the  amount  of  power  expended  by  the  control  input  u{t). 

It  is  well  known  [5]  that  the  above  problem  has  a  unique  solution  Uo{t)  which  can  be 
expressed  in  state  feedback  form  as 

Uoit)  =  -B*K{t)xoit),  (5.9) 

where  K{t),  for  each  t,  is  a  linear  bounded  transformation  from  H  into  H  and  is  absolutely 

continuous  in  t,  0  <  t  <  T.  Furthermore,  K{t)  satisfies  the  operator  Riccati  equation 

^[K{t)x,y]  =  -[Rx,y]-[K{t)x,Ay]-[K{t)Ax,y]  (5.10) 

+[K{t)BB*K{t)x,  y].  K{t)  =  0  Vx,  j/  e  V{A) 
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Hence,  if  we  can  formulate  the  distributed  active  noise  cancellation  problem  as  a  linear 
quadratic  regulator  problem,  we  will  be  able  to  use  the  above  result  to  find  the  unique 
solution  to  the  distributed  active  noise  cancellation  problem.  Furthermore  if  such  a  solution 
exists,  it  will  be  expressible  in  feedback  form. 

The  acoustic  system  that  we  focus  on  in  this  chapter  and  its  governing  equations  are 
presented  in  the  next  Section.  Our  ultimate  goal  is  to  represent  these  equations  as  abstract 
differential  equations  of  the  type  discussed  in  Section  5.2.2. 


5.4  Dynamics  of  Two  Deterministic  Acoustic  Systems 


The  governing  equations  of  the  acoustic  systems  that  we  focus  on  in  this  chapter  are  pre¬ 
sented  in  this  Section.  The  second  order  wave  equation  for  pressure  is  broken  up  into  two 
first  order  partial  differential  equations,  so  that  it  can  be  later  represented  as  an  abstract 
differential  equation  in  Section  5.5. 

Sound  essentially  involves  a  weak  motion  of  a  non-viscous  fluid  from  an  initial  state  of 
rest  [16].  Particle  velocity  u(t,x)  and  pressure  p{t,x)  are  the  two  physical  quantities  whose 
values  completely  characterize  any  acoustic  system,  where  the  underline  denotes  a  vector¬ 
valued  quantity;  t  is  the  continuous- time  variable;  and  x  is  the  spatial  variable.  These  two 
quantities  are  related  according  to  the  following  partial  differential  equations: 


V  p{t,x)  = 


dv(t,x)  , 

-Po . .  +  lit,  x) 

1  dp{t,x) 


V.v{t,x)  =  - 

C-^po  ot 


(5.11) 

(5.12) 


where  f^t,x)  is  an  external  force  field;  po  is  the  ambient  density;  c  is  the  speed  of  sound; 
V.  u  is  the  divergence  of  v;  and  VP  is  the  gradient  of  P.  The  first  equation  is  obtained  by 
applying  Newton’s  law  to  a  small  control  volume,  and  the  second  equation  is  obtained  by 
applying  the  conservation  of  mass  principle  to  this  volume.  The  familiar  wave  equation  for 
p{t,x)  can  be  derived  from  the  above  equations  by  differentiating  eq.  (5.12)  with  respect 
to  time,  taking  the  divergence  of  eq.  (5.11),  and  eliminating  V.(^"^’--^-)  to  get: 


1  d'^p{t,x) 
^  dt^ 


V'^p(t,x)  -  V.  (/(t,x)) 


(5.13) 
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,  where  is  the  three  dimensional  Laplacian.  Equations  (5.11)  and  (5.12)  are  two 
first  order  partial  differential  equations  equivalent  to  the  second  order  partial  differential 
equation  (5.13).  For  us,  it  is  more  convenient  to  work  with  the  following  arrangement  of 
(5.11)-(5.12): 


+  Bu{t), 


(5.14) 


r 


Figure  5.1:  Enclosed  region  in  with  boundary  F. 

Next  let  us  specialize  the  above  equations  to  an  enclosed  region  in  B?.  We  shall  do 
this  specialization  for  two  types  of  boundary  conditions.  The  first  kind  of  boundary  to 
be  considered  is  a  perfectly  reflecting  one,  i.e.  the  normal  component  of  velocity  at  the 
boundary  surface  is  assumed  to  be  zero.  The  second  kind  of  boundary  to  be  considered 
is  an  absorbing  one,  i.e.  the  ratio  of  the  pressure  at  the  boundary  surface  to  the  normal 
component  of  particle  velocity  at  the  boundary  surface  is  assumed  to  be  a  known  positive 


114 


real  number.  This  ratio  is  called  the  acoustic  impedance  of  the  the  boundary;  hence,  in 
the  second  case  we  are  assuming  that  the  impedance  of  the  boundary  is  purely  resistive. 
The  most  general  type  of  boundary  condition  is  obtained  when  this  ratio,  in  the  frequency 
domain,  is  a  complex  valued  function  where  w  is  the  frequency  variable.  The  most 

general  case  is  not  considered  in  this  chapter. 

5.4.1  Purely  Reflecting  Boundaries 

In  Fig.  5.1,  let  Q  be  a  bounded  open  domain  in  with  smooth  boundary  F.  Furthermore, 
assume  that  eq.  (5.14)  is  satisfied  for  aU  the  interior  points  of  with 

Vn{t,x)  =  0  X  e  r,  (5.15) 

where  is  the  normal  component  of  velocity  at  the  boundary.  Finally,  assume  that  the 
following  initial  conditions  are  given: 

■P(0,2.)  =  -Fb(£)  and  u(0,s)  =  Vo{x). 

The  total  acoustic  energy  of  the  system  in  Fig.  5.1  at  time  t,  denoted  by  E{t),  is 

=  i  +  i  f 

where  the  first  integral  is  the  total  acoustic  potential  energy  and  the  second  integral  is  the 
total  kinetic  energy. 

With  the  boundary  condition  eq.  (5.15)  and  with  no  input  u{t),  it  can  be  easily  shown 
[16]  that  E{t)  is  constant  cis  a  function  of  t.  This  is  consistent  with  our  assumption  that 
the  boundary  surfaces  are  lossless  and  that  there  is  no  input  to  the  system;  hence,  there  is 
no  mechanism  to  increase  or  decrease  the  acoustic  energy  in  this  case. 

5.4.2  Purely  Absorbing  Boundaries 

Again  in  Fig.  5.1,  let  fl  be  a  bounded  domain  in  with  smooth  boundary  F.  Furthermore, 
assume  that  eq.  (5.14)  is  satisfied  for  all  the  interior  points  of  fl  with  the  following  boundary 
condition 

pit,x)  = +Z{x)  Vn{t,x)  xeT,  (5.16) 
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where  Vn  is  the  normal  component  of  the  velocity  field  at  the  boundary,  and  Z(x)  is  a 
positive  real  number.  Again,  we  assume  that  the  following  initial  conditions  are  given: 

^(0,1.)  =  Pq{x)  and  u(0,^)  =  Uo(3:)- 

The  purpose  of  active  noise  cancellation  is  to  choose  the  signal  u{t)  in  such  a  way  that 
the  resulting  force  field  /(t,x)  does  negative  work  on  the  system  and  thereby  reduces  the 
acoustic  energy  of  the  system.  Naturally,  one  would  like  to  decrease  E{t)  without  expending 
too  much  power;  hence,  a  reasonable  performance  criterion  for  an  active  noise  cancellation 
system  is  the  sum  of  the  acoustic  energy  E{t)  and  the  power  expended  by  the  control 
input  u{t).  We  would  like  to  formulate  the  active  noise  cancellation  problem  as  an  abstract 
quadratic  regulator  problem  with  the  cost  function  equal  to  this  sum.  To  this  end,  we  first 
need  to  derive  an  abstract  differential  equation  representation  for  the  acoustic  system  of 
Figure  5.1. 


5.5  Abstract  Representation  of  Two  Deterministic  Acous¬ 
tic  Systems 


The  goal  of  this  Section  is  to  derive  an  abstract  differential  equation  representation  for  the 
acoustic  systems  of  Fig.  5.1  subject  to  the  boundary  condition  specified  by  eq.  (5.15)  or 
specified  by  eq.  (5.16).  This  derivation  is  achieved  by  showing  that  the  differential  operator 


A  = 


0  -c  V. 
— c  V  0 


(5.17) 


is  an  infinitesimal  generator  of  an  absolutely  continuous  semigroup  in  an  appropriate  Hilbert 
space.  Once  this  abstract  representation  is  established,  it  will  be  easy  to  show  that  the 
noise  cancellation  problem  in  Fig.  5.1  is  equivalent  to  a  quadratic  regulator  problem.  The 
existence  and  uniqueness  of  the  solution  to  the  noise  cancellation  problem  wiU  then  follow 
from  the  existence  and  uniqueness  theorems  for  the  solutions  of  the  abstract  quadratic 
regulator  problem. 
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The  first  step  toward  establishing  the  abstract  formulation  is  to  choose  the  underlying 
Hilbert  space  and  the  domain  of  A  in  this  space.  For  this,  let  us  define 

=  |/  e  £2(J^)1^7  exists  and  for  i  =  1, 2,  sj  , 

where  £2(f^)  is  the  class  of  square  integrable  functions  with  respect  to  the  Lebesgue  measure 
that  are  defined  on  fi,  and  we  further  define 


M  =  H^{a)  X  H\n)  X  H\Q.)  X  H\n). 


(5.18) 


The  underlying  Hilbert  space  H  is 


H  =  C2{n)  X  £2(0)  X  £2(0)  X  £2(11), 


(5.19) 


with  the  following  inner  product 


i=l 

5.5.1  Purely  Reflecting  Boundaries 

In  the  case  of  purely  reflecting  boundaries,  the  linear  operator  A  is  deflned  as: 


— c  V. (v) 
-c  V(p) 


P(A)  = 


P  ,,,  dv{x) 


eViA) 


0;  Vx  G  r  V  , 


(5.20) 


(5.21) 


where  is  the  normal  component  of  the  velocity  field  at  the  boundary.  Using  the 

divergence  theorem  in  R^,  we  can  show  that  A*  equals  the  negative  of  A: 


9i  fi 

,A 


51  -cV.(/,) 

£,  -cV(/i) 


V  €  ViA) 


=  -c  J^gi{x)V .  {f^{x))dx  -  c  {fi)]fi3dx  (5.22) 


=  ^  (9i)]mdx-c  J^-^{gJi)ds  + 

c  f  /iV.  (g  )dx 
JQ 

c,V.  (5^) 

—  5 

<^^91)  I2 
LL  1  I  ^  Mh 


(5.23) 


(5.24) 


where  in  transition  from  eq.  (5.22)  to  eq.  (5.23)  we  have  used  the  divergence  theorem  in  B? 

Q  j 

and  the  fact  that  on  the  boundary  =  0  .  This  is  essentially  a  three  dimensional  version 
of  the  integration  by  parts  that  was  done  to  get  eq.  (5.5)  in  Section  5.2.2.  The  domain 
of  operator  A  defined  in  eq.  (5.21)  is  dense  in  H,  since  C^(fi),  the  set  of  all  infinitely 
differentiable  function  with  compact  support  in  fi,  is  dense  in  £2(fi)  as  shown  in  [1]. 

We  see  that  the  surface  integral  in  equation  (5.24)  must  equal  zero;  hence,  A*  is 


V{A*)  = 


c  V.  {g  )  gi 

=  -A 

c  V(5i)  g 


dgJx) 

e  M\  =  0  Vx  e  r 

on 


e  x>(A*) 


Closedness  of  the  operator  A  follows  immediately  since  A*  =  —  A  and  all  adjoint  opera¬ 
tors  are  closed  [33].  Next,  using  the  divergence  theorem,  we  show  that  A  is  dissipative: 


/i  /i 


i.2  i.2 


V.  (Z^)  /i 


=  -c  f  V.  (/  )/idx-c  f  [V(/i),/  jjqsdx  (5.25) 

JQ 

=  c^[/2,V(/i)]^j3dx  -  c^[V(/i),/2]R3dx  (5.26) 


=  0, 


(5.27) 


where  in  transition  from  eq.  (5.25)  to  eq.  (5.26),  we  have  used  the  divergence  theorem  and 

Q  y 

the  fact  that  on  the  boundary  =  0. 
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Combining  eq.  (5.27)  with  A“  =  —  A,  we  get 


(5.28) 

(5.29) 


0  =  [AfJ]  +  [f,Af]  'ifeViA) 

0  =  +  =  0  V5GP(A*). 

We  have  therefore  shown  that  >1  is  a  closed,  linear  and  dissipative  operator  with  dense 
domain  in  H .  Moreover,  A  satisfies  (5.28)  and  (5.29).  Applying  Theorem  2,  it  follows  that 
A  generates  a  strongly  continuous  semigroup  T{t)  on  H  such  that 

lir(0/llH=li/il//  v/G^. 

From  the  last  expression  we  see  that  T{t)  is  norm  preserving.  This  can  be  interpreted  as 
saying  that  the  acoustic  energy  of  the  system  in  Fig.  5.1  with  the  boundary  condition 
specified  by  eq.  (5.15)  is  constant,  since  the  norm  in  H  is  proportional  to  the  acoustic 
energy  in  the  enclosure. 


5.5.2  Purely  Absorbing  Boundaries 

In  the  case  of  purely  absorbing  boundaries,  the  linear  operator  A  is  defined  as: 


P 

-c  V.  (u) 

P 

A 

= 

V 

P. 

-c  V(p) 

V 

V(A)  =  { 


P 


eV(A) 


G  M\  p{x)  =  Z{x)  Vx  G  r 


(5.30) 


(5.31) 


Using  the  divergence  theorem  in  72^,  we  find  A*  and  its  domain: 


9i 

/i 

5i 

-cv.(/2) 

h 

,A 

= 

9 

V 

92 

H 

9.2  _ 

-c  V(/i) 

H 

to 

1 _ 

=  -c  J^gi{x)V.  {l^{x))dx  -  c  (/i)]fl3 

=  (9i)]R^dx  -  c  J^^{l^gi)ds  + 


GP(A) 
dx  (5.32) 
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■-lhV.(g^)dx-cJ^^(g^A)d, 


(5.33) 


c  V.  (£,)  A 

—  9 

‘^^(51)  L2 

LL  1  I Mh 

a  +  i  I;(22A)<'^}  .  (5-34) 

where  in  transition  from  eq.  (5.32)  to  eq.  (5.33)  we  have  used  the  divergence  theorem  in 
R^.  The  domain  of  operator  A  defined  in  eq.  (5.31)  is  dense  in  H,  since  C'^(fi)  is  included 
in  the  domain  of  A,  and  is  dense  in  .ff  as  shown  in  [1]. 

From  eq.  (5.34),  we  see  that  the  expression  in  the  brackets  must  equal  zero.  By  equating 
this  expression  to  zero  and  recalling  eq.  (5.16),  we  get  the  following  expression  for  A*: 


<^^■(9.2) 

c^{9i) 


=  -A 


€  V{A*) 


P(A*)  = 


%2(^) 


G  M|  g^ix)  =  ~Z(x)  Vx  G  r  ^  . 


Using  the  above  derivation  for  A*,  we  see  that  operator  A  is  the  adjoint  of  the  linear 


operator  B  defined  as 


C  V.  (£) 

cV(/) 


eV{B) 


V{B)  = 


f  dg(x) 

G  M\  fix)  =  -Z(x)  Vx  G  r 


(5.35) 


(5.36) 


Operator  A  must  then  be  closed,  since  the  adjoint  of  B  is  closed.  Recall  that  all  adjoint 
operators  are  closed  [33]. 

Next  using  the  divergence  theorem,  we  show  that  A  is  dissipative: 


A  A 


-cV.(/„)  A 


Z2  i-2 


-cV(A) 
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=  — c 

[  V.  if^)hdx  -  c  1 

IQ  JQ 

[V{f,)J^]R3dx  (5.37) 

=  c/ 

Jq 

[/2>V(/i)]fi3c!a:  -Cj 

—c 

ljyih),f:2Wdx 

(5.38) 

=  —c 

w 

I  o 

II  VI 

where  in  transition  from  eq.  (5.37)  to  eq.  (5.38)  we  have  used  the  divergence  theorem. 

To  get  the  last  expression,  we  have  used  the  assumption  that  on  the  boundary  fi{x)  = 

Zix)i^f2{x). 

The  operator  A*  is  also  dissipative,  and  this  can  be  shown  following  steps  very  similar 
to  the  one  used  to  show  that  A  is  dissipative. 

We  have  shown  that  yi  is  a  closed,  linear  and  dissipative  operator  with  dense  domain 
in  H.  Moreover,  A*  has  also  been  shown  to  be  dissipative.  Applying  Theorem  1,  it  follows 
that  A  generates  a  strongly  continuous  semigroup  T{t)  on  H  such  that 

||T(t)/i|H<||/||H  V/€F. 

The  last  expression  can  be  interpreted  as  saying  that  the  mechanical  energy  of  the  system  in 
Fig.  5.1  with  the  boundary  condition  specified  by  eq.  (5.16)  is  decreasing.  This  is  consistent 
with  the  fact  that  the  boundaries  in  this  case  are  assumed  to  be  absorbing. 

We  have  shown  that  the  governing  equation  for  the  acoustic  system  of  Fig.  5.1,  with 
purely  reflecting  boundaries  or  purely  absorbing  boundaries,  can  be  represented  as  an  ab¬ 
stract  differential  equation  of  the  form  ^x{t)  =  Ax(t)  -f  Bu{t),  where  A  is  the  infinitesimal 
generator  of  a  strongly  continuous  semigroup. 

We  now  write  down  the  abstract  differential  equation  representation  for  the  system 
in  Fig.  5.1,  subject  to  the  boundary  condition  specified  by  eq.  (5.15)  or  by  eq.  (5.16). 
Assuming  that  n(.)  is  an  element  of  >V„  =  £2((0,T);  Hu)  and  recalling  that  5  is  a  bounded 
linear  operator  from  the  Hilbert  space  Hu  into  H,  we  arrive  at  the  following  abstract 
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representation  of  the  acoustic  system  in  Fig.  5.1: 

=  Ax{t)  +  Bu{t)  (5.39) 

a;(0)  is  given  in  if, 

where  H  is  the  underlying  Hilbert  space.  Moreover,  the  operator  A  is  defined  by  eq.  (5.17), 
and  its  domain  is  defined  by  eq.  (5.21)  in  the  case  of  purely  reflecting  boundaries  and  by 
eq.  (5.31)  in  the  case  of  purely  absorbing  boundaries. 

The  distributed  active  noise  cancellation  problem  can  now  be  viewed  as  a  linear  quadratic 
regulator  problem  applied  to  the  system  governed  by  eq.  (5.39).  Specifically,  the  objective 
is  to  choose  the  input  u{t)  such  that  the  the  cost  functional 

X*  T 

J  =  f  E{t)dt+  f  [u{t),u{t)]H^dt  (5.40) 

Jo  Jo 

=  [02;(t),2;(t)]Hdt  +  ^  [u{t),u{t)]ffjt  (5.41) 

is  minimized,  where  E{t)  is  the  mechanical  energy  of  the  system  at  time  t,  and  [u{t),  u(t)] 
is  the  power  expended  by  the  control  input  at  time  t.  From  eq.  (5.9),  we  see  that  a  unique 
minimizing  control  Uo{t)  in  Wu  exists  and  can  be  expressed  in  the  state  feedback  form  as 

Uo{t)  =  -B*K{t)xo{t), 

where  the  linear  bounded  operator  I{{i)  satisfies  the  operator  Riccati  equation 

^[K{t)x,y]  =  -[:^x,y]-[K{t)x,Ay]-[K{t)Ax,y] 

+[K{t)BB*K{t)x,y].  K(t)  =  0  Vx,y  e 'D(A)  (5.42) 

In  summary,  by  formulating  the  distributed  active  noise  cancellation  problem  as  a  lin¬ 
ear  quadratic  regulator  problem,  we  have  been  able  to  show  that  the  distributed  active 
noise  cancellation  problem  has  a  unique  solution.  Furthermore,  we  have  given  an  explicit 
expression  for  this  solution  as  a  linear  feedback  of  the  state  of  the  system. 

The  main  shortcoming  of  the  above  solution  is  that  the  control  signal  is  generated  based 
on  the  exact  knowledge  of  the  entire  state  of  the  system.  In  a  practical  system,  the  state 
is  not  known  exactly  and  the  control  signal  must  be  generated  based  on  the  observed  data 
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alone  (typically  the  observed  data  is  the  output  of  one  or  more  microphones  that  are  placed 
in  the  volume  of  interest).  It  turns  out  that  the  optimal  control  in  this  case  is  obtained 
by  feeding  back  the  estimate  of  the  state.  The  resulting  estimation  and  control  problem 
for  infinite-dimensional  linear  systems  is  considered  in  the  next  Section.  We  will  show  that 
the  distributed  active  noise  cancellation  problem  in  which  the  control  signal  is  based  on 
the  observed  data  alone  can  be  viewed  as  a  stochastic  quadratic  regulator  problem.  The 
noise  cancellation  problem  can  then  be  solved  using  the  stochastic  optimal  control  theory 
for  distributed  parameter  systems. 

5.6  Stochastic  Optimal  Control  for  ANC 

In  this  Section,  the  stochastic  analog  of  the  linear  quadratic  control  problem  of  Section  5.3 
is  presented.  For  this  purpose  a  measure  theoretic  structure  is  introduced  on  the  top  of  the 
topological  structure.  The  measure  theory  is  nonclassical  in  the  sense  that  the  measures 
are  only  finitely  additive  on  the  field  of  cylinder  sets.  The  notion  of  a  weak  random  variable 
suffices  for  the  stochastic  extension  of  the  control  problem  of  Section  5.3  with  a  crucial 
notion  being  that  of  “white  noise”.  The  development  of  the  stochastic  control  problem 
presented  here  is  entirely  due  to  Balakrishnan  [5]. 

5.6.1  Preliminaries 

The  measure  theoretic  background  that  is  needed  to  understand  our  formulation  of  the 
stochastic  quadratic  regulator  problem  is  presented  in  this  Section.  Measures  on  cylinder 
sets,  weak  random  variable,  “white  noise”  processes,  and  linear  infinite-dimensional  systems 
with  white  noise  input  are  discussed. 

Measure  on  Cylinder  Sets 

Cylinder  sets  and  cylinder  set  measures  are  introduced  in  this  Section.  We  start  with  H  a 
separable  real  Hilbert  space.  Let  us  take  any  finite  dimensional  subspace  Hm  in  H.  By  a 
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cylinder  set  we  mean  any  set  of  the  form 


where  5  is  a  Borel  subset  of  ‘Hm,  and  is  the  orthogonal  complement  of  Hm-  The  Borel 
set  B  is  then  called  the  “base”  of  the  cylinder,  and  Tim  is  called  the  “base”  space.  It  can 
be  shown  that  the  class  of  cylinder  sets  C  forms  a  field  of  sets.  Furthermore,  Ti  itself  is  a 
cylinder  set,  and  the  Borel  sets  of  Ti  are  the  smallest  cr-algebra  that  contains  C. 

Next  we  consider  measures  on  the  field  of  cylinder  sets,  these  measures  will  be  referred 
to  as  “cylinder  set  measures”  and  denoted  by  /r.  Let  Z  be  a  cylinder  set  with  base  B  in 
Tim-  Then  by  definition 

H{Z)  =  Vm{B), 

where  Vm{-)  is  a  countably  additive  probability  measure  on  the  a-algebra  of  Borel  subsets 
of  Tim-  Furthermore  Vm  must  satisfy  the  following  condition. 

Compatibility  Condition:  In  order  that  p,  be  well  defined,  it  is  necessary  that  if 

Z  =  B  Tim  ~  {B  +  Tip)  +  {Tim  +  'Bp)^, 

where  Tip  is  orthogonal  to  Tim,  then  Vm{B)  =  Vm+p{B  +  Tip)  where  Vm-\-pi-)  is  a  countably 
additive  probability  measure  on  the  Borel  subsets  Tim  +  Tip. 

It  is  important  to  realize  that  a  cylinder  set  measure  is  not  necessarily  countably  additive 
on  the  field  of  cylinder  sets.  Consequently,  a  cylinder  set  measure  can  not  in  general  be 
extended  to  the  Borel  sets  of  Tt. 

The  following  example  of  a  cylinder  set  measure  wiU  be  of  prime  importance  in  our 
development.  Let  R  be  any  self  adjoint  nonnegative  definite  operator  mapping  Ti  into  Ti  - 
We  define  the  measure  on  the  Borel  sets  of  any  finite  dimensional  space  Tim  in  the  following 
way.  Let  ci, . . . ,  be  an  orthonormal  basis  in  Tim-  The  Borel  sets  in  Tim  are  in  one-to-one 
correspondence  with  the  Borel  sets  in  via  the  transformation 


where  x  6  Ti-m  a-nd 


€  -R'".  The  measure  on  Borel  sets  of  R”^  is  defined  to  be  a 


[x,em]  J 

Gaussian  measure  with  a  moment  matrix  given  by 

Tij  —  [^e^,  e^]  i,  j  —  1, . . . ,  TTi. 

The  measure  of  any  Borel  set  in  Hm  is  now  defined  to  be  the  Gauss  measure  of  the  cor¬ 
responding  Borel  set  in  R™-.  This  measure  is  independent  of  the  particular  basis  used.  It 
can  be  verified  that  the  compatibility  condition  is  also  satisfied.  The  resulting  cylinder 
set  measure  is  usually  denoted  by  fi  and  is  referred  to  as  the  measure  induced  by  R.  In 
the  special  case  that  R  is  the  identity  operator,  the  resulting  cylinder  set  measure  on  Ti  is 
referred  to  as  the  standard  Gauss  measure  on  H. 

Weak  Random  Variables 

The  usual  definition  of  a  random  variable  requires  a  probability  triple  where  fi  is 

an  abstract  space,  5  is  a  <7-field  of  subsets  of  fi,  and  p  is  a  countably  additive  probability 
measure  on  B.  A  random  variable  is  then  merely  any  function  (usually  with  range  in  iZ") 
that  is  measurable  with  respect  to  B.  Typically,  there  is  a  field  of  subsets  of  fl,  call  it  F, 
which  generates  the  a-field  B.  In  this  case,  specifying  p  as  a  countably  additive  set  function 
on  F  uniquely  determines  p  on  B.  In  other  words,  p  can  be  uniquely  extended  from  F  to 
B,  provided  that  p  is  countably  additive  on  F. 

A  cylinder  set  measure  p  can  not  necessarily  be  extended  to  the  Borel  subsets  of  Ti, 
since  a  cylinder  set  measure  is  not  necessarily  countably  additive  on  the  field  of  cylinder 
sets.  This  implies  that  if  5  is  a  Borel  set  in  Ti,  the  measure  of  B  is  not  generally  defined 
based  on  the  specification  of  p  as  a  cylinder  set  measure.  This  in  turn  implies  that  we  can 
not  adopt  the  usual  definition  of  a  random  variable  as  any  function  /(.)  defined  on  71  that 
is  measurable  with  respect  to  the  Borel  sets  of  H.  The  probability  of  the  event  that  /(.)  is 
an  element  of  a  Borel  set  in  the  range  space  is  simply  not  defined  in  general. 


Hence,  the  notion  of  a  “weak  random  variable”  is  introduced.  Let  [H,  C,  pi]  denote  the 
probability  triple  for  which  is  a  real  separable  Hilbert  space,  C  is  the  field  of  cylinder 
sets  in  and  is  a  cylinder  set  measure  there  on.  We  shall  denote  points  in  by  u;  in 
this  context.  Let  /(.)  be  a  function  from  71  into  another  Hilbert  space  Tir-  /(.)  is  called  a 
weak  random  variable  if  for  each  <!>  in  Tir ,  is  defined  and  countably  additive  on  inverse 
images  of  one- dimensional  Borel  sets  of  the  form 

=  {u;|w  €  n,  [f{w),(l>]  €  E}, 

where  jE  is  a  one-dimensio'nal  Borel  set,  and  ^  is  fixed  element  of  Hr-  This  definition 
essentially  ensures  that  inverse  images  of  cylinder  sets  in  Tir  are  cylinder  sets  in  7i. 

The  canonical  example  of  a  weak  random  variable  is  f{w)  =  Lw  where  L  is  any  bounded 
linear  transformation  from  Ti  into  Hr-  Let  fii,  denote  the  cylinder  set  measure  induced  by 
L  on  Tir’,  that  is,  for  any  cylinder  set  C  in  Tir,  we  define 

=  fj,{{w\w  e  H  and  Lw  G  C})  =  /i(L"’^(C)), 

note  that  L~^{C)  is  a  cylinder  set  in  H.  If  fi  is  induced  by  the  self  adjoint  nonnegative 
operator  R  in  H,  then  m  is  the  cylinder  set  measure  induced  by  LRL*  in  Ht- 

Second  Order  Characterization 

Second  order  characterization  of  weak  random  variables  are  reviewed  in  this  Section.  No¬ 
tions  of  mean,  covariance,  and  cross  covariance  for  weak  random  variables  are  presented. 
Let  ^  be  a  Hilbert  space  valued  weak  random  variable.  We  shall  say  it  has  a  finite  first 
moment  if 

1.  £{1[^,  (;i»]|}  <  oo  for  every  (j)  in  H,  and 

2.  is  continuous  in  <f>, 

where  E{.)  denotes  the  expectation.  Then  there  exists  an  element  m  in  Ti  such  that 
jE{[^,  <f>]}  =  [m,  4>].  We  shall  use  the  notation  E{^)  =  m  and  call  m  the  mean  of  Similarly, 
we  shall  say  that  ^  has  a  finite  second  moment  if 
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1.  J5{[^,  <  oo  for  every  (j>  in  7i,  and 

2.  E{[^,<f)]‘^}  is  continuous  in  4>- 

Let  I  =  ^  -  m,  where  m  =  then  ^  also  has  a  finite  second  moment,  with  first  moment 

equal  to  zero.  For  any  two  elements  x,y  in  H,  y]}  is  a  continuous  bilinear  form 

over  H.  Hence,  there  exists  a  bounded  linear  nonnegative  self  adjoint  operator  R  :  Ti  — >•  H, 
called  the  covariance  of  such  that 

=  [R^,y]- 

In  particular,  if  ^  is  a  Gaussian  weak  random  variable  so  that 

E{exp{i[^,<j)]))  =  exp{-^-^^^)exp(i[m,(f)]), 

it  follows  that  ^  has  a  finite  second  moment,  with  covariance  R  and  mean  m. 

The  concept  of  cross  covariance  for  weak  random  variables  is  introduced  next.  Let  p 
be  two  Hilbert  space  valued  weak  random  variables  with  finite  second  moment  and  zero 
first  moments.  Then  Q{x,y)  =  £'([^,  a;][»?,  y])  defines  a  continuous  bilinear  map,  and  hence 
there  exists  a  bounded  linear  operator  S  such  that  Q{x,y)  =  [5x,y].  We  shall  write 

If  R^^  denotes  the  covariance  of  and  i?,,,,  that  of  rj,  then  it  is  consistent  to  use  the  following 
notation,  analogous  to  the  finite  dimensional  case, 


%  =  E{^n,  (5.43) 

Rin  =  (5.44) 

Rr,,  =  E{rjp*).  (5.45) 

In  the  case  that  R^^  is  nuclear,  we  can  define  jE(||^|P)  as 

OO 

£(l|{||'')  =  (5-«) 

k—1 

oo 

=  (5.47) 

fc=i 

=  rr.(%),  (5.48) 
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where  {4>k}  is  an  orthonormal  basis  for  H.  Note  that  the  nuclearity  of  guarantees  that 
Tr.{R^^)  exists  and  is  finite. 

White  Noise 

The  crucial  notion  of  white  noise  is  introduced  in  this  Section.  Let  H  be  a  real  separable 

Hilbert  space,  and  let  W  =  £2((0,T);?f)  with  0  <  T  <  +oo.  Note  that  >V  is  also  a 

separable  Hilbert  space  with  the  following  inner  product 

fT 

[u,u]w=  /  [u{t),v{t)]-Hdt.  u,v£yV 
Jo 

Let  denote  the  standard  Gauss  measure  on  >V,  and  let  w  denote  any  element  of  W.  That 
is  to  say,  fj.  is  the  measure  induced  by  the  identity  operator  on  W.  Then,  the  “function 
space”  process  so  obtained  is  called  “white  noise”.  Each  w(t)  in  W  is  called  a  white  noise 
sample  function.  For  any  u{t)  in  W,  [tt7,u]w  defines  a  zero  mean  Gaussian  random  variable 
with  variance  [u,  u]w  H  is  natural  that  any  physical  observation  such  as  [tr,  u]w  be  a 
random  variable  and  not  merely  a  weak  random  variable. 

There  is  a  close  analogy  between  the  Karhunen-Loeve  expansion  of  a  white  noise  process 
in  the  finite-dimensional  case  and  the  projections  of  an  infinite-dimensional  white  noise 
process  along  an  orthonormal  basis  in  W.  Recall  that  if  n{t)  is  a  unit  variance,  scalar, 
Gaussian  white  noise  process,  and  orthonormal  basis  for  £2((—oo,  -f-oo);  R), 

then 

Cfc(n)  =  ln(t),t/;k(t)]c2 

are  mutually  independent  zero  mean,  unit  variance  Gaussian  random  variables. 

Similarly,  if  is  an  orthonormal  basis  in  W  and  w{t)  is  a  standard  white  noise 

process  defined  on  W,  then 

77fc(w)  =  [w{t),<pk{t)]w 

are  mutually  independent  zero  mean,  unit  variance  Gaussian  random  variables. 

It  is  important  to  note  that  the  definition  of  the  white  noise  process  restricts  the  sample 
space  to  be  the  class  of  square  summable  functions,  since  W  —  £2((0,T);'?f).  This  restric¬ 
tion  becomes  important  when  we  consider  systems  governed  by  partial  differential  equations 
with  white  noise  input. 
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Abstract  Differential  Equations  with  White  Noise  Input 

By  invoking  the  theory  of  semigroups  we  obtain  a  natural  extension  of  linear  systems  with 
finite-dimensional  state  space  to  the  infinite-dimensional  case  embracing  in  particular,  sys¬ 
tems  described  by  partial  differential  equations  with  white  noise  input.  Let  A  denote  the 
infinitesimal  generator  of  a  strongly  continuous  semigroup  T(t)  over  the  real  separable 
Hilbert  space  Let  W  =  £2((0,T);?f)  with  T  finite.  Let  Tin,  where  subscript  n  stands 
for  noise,  be  another  real  separable  Hilbert  space,  and  Wn  =  £2((0,r);  Let  w{t)  be 
a  standard  white  noise  process  in  W„,  and  let  B  denote  a  bounded  linear  transformation 
mapping  into  V..  Consider  the  abstract  differential  equation 

=  Axii)  d-  Bw(t).  0  <  t  <T;  x(0)given 
at 

For  each  w(t)  in  W„,  we  have  seen  (Theorem  4,  Section  5.2.2)  that  this  equation  has  a 
unique  solution  such  that  [x(t),y]  is  absolutely  continuous  for  each  y  in  T>(A*).  Moreover, 
the  solution  is  given  by 

x(t,  tu)  =  T(t)x(0) -h  f  T{t  —  s)Bw(s)ds, 

Jo 

where  we  have  written  the  solution  as  x{t,w),  instead  of  x{t),  to  emphasize  the  dependence 
of  the  solution  on  w{t).  Having  discussed  linear  infinite-dimensional  systems  with  white 
noise  input,  we  wiU  present  the  stochastic  control  problem  for  these  systems  next.  Our  goal 
is  to  ultimately  formulate  the  active  noise  cancellation  problem  as  such  a  control  problem. 

5.6.2  Stochastic  Quadratic  Regulator  Problem 

Stochastic  analog  of  the  quadratic  regulator  problem  of  Section  5.3  is  presented  in  this 
Section.  The  main  difference  is  that  the  control  has  to  be  based  on  the  observed  data  only; 
hence,  the  control  is  forced  to  be  “feedback”  or  “closed  loop”.  Consider  the  state  evolution 
equation 

d 

—x{t,-w)  =  Ax{t^w)-\-Bu{t,w)-{-Fw{i)  0  <  t  <  T;  x(0)  =  0,  (5.49) 

at 

where  A  is  the  infinitesimal  generator  of  a  strongly  continuous  semigroup  T{t)  over  the 
Hilbert  space  ?f,  w{t)  is  a  standard  white  noise  process  in  W„  =  £2((0,  T);  ?f„),  Hn  is  a 
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separable  real  Hilbert  space,  and  is  a  linear  bounded  transformation  mapping  Tin  into 
H.  The  control  process  is  u{t,w)  where  u{.,w)  is  in  =  £2((0,r);7fu),  Tiu  is  a  separable 
Hilbert  space,  and  we  assume  that  5  is  a  Hilbert-Schmidt  operator  mapping  Tin  into  Ti. 
The  observation  process  is  defined  by 

y{t,w)  =  Cx{t^w) Gw{t)  0  <  f  <  T,  (5.50) 

where  C  is  assumed  to  be  a  linear  and  bounded  mapping  from  Ti  into  Tig-,  where  Tig  is  a 
real  separable  Hilbert  space.  Finally  G  maps  Tin  into  Tig,  with 

GG*  =  Identity  (5.51) 

and 

FG*  =  0.  (5.52) 

Note  that  eq.  (5.51)  implies  that  Gw{t)  is  a  white  noise  process,  and  eq.  (5.52)  implies 
that  the  observation  noise,  Gw{t),  is  independent  of  the  process  noise  Fw(i): 

The  control  is  required  to  be  such  that 

E{u{t,w)u{t,w)*) 

is  nuclear,  and  further  be  determined  for  each  “t”  in  terms  of  the  observations  up  to  time 
“t”  by  a  linear  transformation;  specifically, 

u{t,w)  =  f  W{t,s)y{s,w)ds  0<t<T, 

Jo 

where  for  each  s  <  t,  IF(t,  s)  is  a  linear  bounded  mapping  from  Tig  into  Tin,  furthermore, 
IF(t,5)  is  strongly  continuous  in  0  <  5  <  t  <  T. 

The  optimal  control  problem  is  to  find  the  control  u{t,  w)  that  minimizes  the  cost  func¬ 
tional 

q{u)  =  f  Tr.E{{Qx(t,w)){{Qx(t,w))*)dt  +  f  Tr.E{u{t,w)u{t,w)*)dt,  (5.53) 
Jq  Jo 

where  Q  is  Hilbert-Schmidt.  We  are  only  interested  in  the  controls  for  which 

/  Tr.E{{u{t,w)u{t,'wy))dt  <  oo. 

Jo 
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Utilizing  the  “separation  principle”,  we  can  separate  the  filtering  problem  and  the  control 
problem.  More  precisely,  the  optimal  control  Uo(t,  w)  is  unique  and  can  be  expressed  in  the 
form 

Uo{t,w)  = —B*  K{t)x{i,w),  (5.54) 

where  the  operator  K{t)  satisfies  the  following  operator  Riccati  equation 

^[K{t)x,y]  -  -[Q''Qx,y]-[K{t)x,Ay]-[K{i)Ax,y]  (5.55) 

A[K{i)BB*K{t)x,y]  A'(7^)  =  0,  Vx,t/6P(A), 


and 


x(t,  w)  =  E{x{t,  w)|y(s,  w)  0  <  s  <t}. 
Furthermore,  x{t,  w)  is  calculated  using  a  Kalman  filter,  i.e. 

^x(t,  w)  =  Ax{t,  w)  +  Pf{t)C*{y{t,  tn)  -  Cx{t,  w))  +  Bu(t,  w) 


x(0,w)  =  0, 


where 

Ff(t)  =  E{(x(t,  w)  -  f  (t,  w)){x{t,  w)  -  x{i,  w))*} 

and  Pf(t)  satisfies 

^[Py(t)x,  y]  =  [Pfit)A*x,  y]  +  [Pf{t)x,  A*y]  +  [FF^x,  y]  (5.56) 

-[Pj{t)C*CPj{t)x,  y]  P(0)  =  A  Vx,ye  V(A*). 


Next  we  show  the  close  relationship  between  the  solution  to  the  optimal  control  problem 
in  the  deterministic  case  and  stochastic  case.  Recall  from  Section  5.3  that  the  deterministic 
optimal  control  problem  is  to  find  u(t)  that  minimizes  the  cost  functional 

fT  rT 

9d(w)  =  /  [Rx(t),x(t)]dt+  /  [w(f),u(t)](it,  (5.57) 

Jo  Jo 

where  iZ  is  a  bounded  nonnegative  definite  operator  used  as  a  weighting  factor.  Setting 
R  =  Q*Q,  we  can  rewrite  qd{u)  as 


Next  using  the  notation  introduced  in  eq.  (5.46),  we  can  rewrite  eq.  (5.53)  as 

q{u)=  j  E{\\Qx{t,w)\\'^]dt+  I  i;{||u(t,  iu)|p}dt,  (5.59) 

JQ  JQ 

since  both  Qx{t,  w)  and  u(t,  w)  have  nuclear  covariance  operators.  Comparing  eq.  (5.58) 
to  eq.  (5.59),  we  see  that  the  only  difference  between  the  two  is  that  the  integrands  in 
the  stochastic  cost  functional  are  the  expected  values  of  the  integrands  in  the  deterministic 
cost  functional. 

According  to  eqs.  (5.9)  and  (5.54),  the  optimal  control  in  the  deterministic  case  is 
given  by 

Uo  =  -B*K{i)x(t),  (5.60) 

and  the  optimal  control  in  the  stochastic  case  is  given  by 

Uo  = —B*  K{t)x{t,w).  (5.61) 

Comparing  eq.  (5.10)  to  eq.  (5.55),  we  see  that  K{t)  in  eq.  (5.60)  is  identical  to  K{t)  in 

eq.  (5.61),  provided  that  R  =  Q*Q.  In  other  words,  the  optimal  control  in  the  deterministic 
case  is  obtained  by  feeding  back  the  state,  x(<),  through  a  linear  operator,  and  the  optimal 
control  in  the  stochastic  case  is  obtained  by  feeding  back  the  state  estimate,  x{t,  w),  through 
the  same  linear  operator. 

The  results  of  this  Section  will  be  very  important  to  us  later,  since  these  results  imply 
that  if  we  can  formulate  the  noise  cancellation  problem  as  a  stochastic  control  problem,  we 
will  be  able  to  show  that  the  noise  cancellation  problem  has  a  unique  solution.  Furthermore 
if  such  a  solution  exists,  it  will  be  expressible  in  feedback  form  according  to  eq.  (5.61). 

The  details  of  the  acoustic  system  that  we  focus  on  and  its  governing  equations  are 
presented  in  the  next  Section.  Our  goal  is  to  ultimately  represent  these  equations  in  the 
format  of  eqs.  (5.49)-(5.50)  and  use  the  results  of  this  Section  to  find  the  closed-loop 
feedback  controller  that  minimizes  the  total  acoustic  energy  of  this  system. 

5.6.3  Abstract  Representation  of  a  Stochastic  Acoustic  System 

The  acoustic  system  that  we  focus  on  is  depicted  in  Fig.  5.2.  In  this  figure,  O  is  an  open 
domain  in  with  smooth  boundary  T,  and  for  concreteness,  we  assume  that  F  is  perfectly 
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Figure  5.2:  Bounded  region  in  with  loudspeaker  and  microphone. 

reflecting.  The  control  signal  u{t)  drives  the  loudspeaker  that  generates  the  canceling  field, 
and  the  output  of  the  microphone  y{t)  is  the  observed  data.  Loosely  speaking,  the  goal 
of  the  noise  cancellation  system  is  to  minimize  the  expected  value  of  the  acoustic  energy 
in  f),  without  expending  too  much  power.  The  main  restriction  is  that  the  control  signal 
at  time  “f’  must  be  generated  based  on  the  observed  data  up  to  time  “t”,  i.e.  based  on 
{y(r) :  r  <  t}. 

First,  we  develop  an  abstract  model  for  the  acoustic  system  of  Fig.  5.2  as  an  infinite¬ 
dimensional  linear  system  driven  by  a  white  noise  process  and  a  control  process.  Our 
approach  is  to  model  the  acoustic  field  in  Q,  as  the  sum  of  two  separate  fields:  the  unwanted 
noise  field  and  the  canceling  field.  The  unwanted  noise  field  is  modeled  as  the  state  of  an 
infinite-dimensional  linear  system  driven  by  white  noise,  and  the  canceling  field  is  modeled 
as  the  state  of  another  infinite-dimensional  linear  system  driven  by  the  control  signal.  These 
two  fields  are  then  added  together  to  obtain  the  residual  acoustic  field  that  is  present  when 
the  noise  canceling  system  is  on. 

Having  represented  the  acoustic  system  of  Fig.  5.2  as  an  infinite-dimensional  linear 
system,  we  use  the  results  of  the  previous  Section  to  find  the  closed-loop  feedback  controller 
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that  minimizes  the  total  acoustic  energy  in  12. 

Unwanted  Noise  Field 

The  unwanted  noise  field  in  is  modeled  as  the  state,  Xd(t),  of  an  infinite-dimensional 
linear  system  driven  by  white  noise,  i.e. 

d 

=  AdXdit)  +  Fdw(t),  (5.62) 

where  Ad  is  an  infinitesimal  generator  of  a  strongly  continuous  semigroup  Td{t)  over  the 
Hilbert  space  H,  with  H  defined  in  eq.  (5.19);  w(t)  is  a  standard  white  noise  process  in 
Wn  =  ^2((0,T);fr);  and  Fd  is  a  linear  bounded  transformation  mapping  H  into  H.  The 
state  vector  Xd{t)  is  composed  of  the  unwanted  pressure  field  and  the  unwanted  particle 
velocity  field  in  the  enclosure,  i.e. 

Pd(t,x) 

Xd{t)  =  ,  (5.63) 

C3u(t,S.) 

where  Pd{t.,x)  is  the  unwanted  pressure  field,  and  is  the  unwanted  velocity  field. 

Equation  (5.62)  determines  the  second  order  statistics  of  the  unwanted  field,  i.e.  it 
determines  E(xd{t,w))  and  E{xd(t,w)x2is,w)). 

Canceling  Field 

The  canceling  field  is  modeled  as  the  state,  Xp{t),  of  another  infinite-dimensional  linear 
system  driven  by  the  control  signal  u{t).  The  canceling  field  is  generated  by  driving  the 
canceling  loudspeaker  with  the  control  input  u{t);  hence,  the  exact  nature  of  this  infinite¬ 
dimensional  system  depends  on  the  dynamics  of  the  loudspeaker  as  well  as  the  dynamics  of 
the  acoustic  enclosure. 

The  simplest  model  for  a  loudspeaker  is  a  collection  of  point  sources  of  pressure  [6]. 
Denoting  the  input  to  the  loudspeaker  by  u{t),  we  model  the  loudspeaker  as  “m”  point 
sources,  located  at  {j?i, . . .,  77m}  respectively.  Each  one  of  these  point  sources  has  the 
strength  r(t)6{r]  —  r]i),  where  r{t)  is  related  to  the  control  signal  u{t)  according  to  the 


134 


following  finite-dimensional  state-space  equations: 

d 

=  A/a:/(t)  -1-  Biu{t) 
at 

r(t)  =  Cixi(t), 


(5.64) 

(5.65) 


where  xi(t)  is  the  (n  x  1)  state  vector;  At  is  the  (n  x  n)  state-transition  matrix;  5;  is  an 
(n  X  1)  vector;  and  C/  is  a  (1  x  n)  row  vector.  Equations  (5.64)-(5.65)  model  the  dynamics 
of  the  loudspeaker. 

In  this  case,  the  canceling  pressure  filed  p{t,x)  satisfies  the  following  wave  equation 

-  r)i),  (5.66) 

1=1 

where  r{t)  is  related  to  the  control  signal  u{t)  according  to  eqs.  (5.64)-(5.65),  and  ^(.)  is  a 
sufficiently  smooth  approximation  to  the  delta  function.  We  claim  that  equation  (5.66)  can 
be  expressed  in  standard  state-space  form  as: 


(5.67) 


where  the  linear  operator  Ap  is 


P 

— cV.  (v) 

P 

Ap 

= 

V 

V 

-cV(p) 

V. 

V{Ap)  = 


P 


€  ViAp) 


(5.68) 


(5.69) 


and  Bp  is  a  bounded  linear  operator  from  into  H  defined  as 


0 

t'^Thi{v)z 
t  hi{v)z 


y  zeR\ 


(5.70) 
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where  is  the  ratio  of  the  speed  of  sound  to  the  ambient  density,  and 


K'n  ~  Vi)  with  hi{r})  G 


(5.71) 


This  claim  can  be  verified  by  substituting  eq.  (5.70)  in  eq.  (5.13)  and  observing  that  the 
resulting  .pressure  field  satisfies  eq.  (5.66). 

The  existence  of  a  function  hi{Tj)  that  satisfies  eq.  (5.71)  can  be  argued  in  the  following 
way.  Let  the  three-dimensional  Fourier  transform  of  /i,(T?)  be  denoted  by  FT{hi{Tj)}  and 
observe  that  eq.  (5.71)  implies 


mhiiv)} 


j{uJi  +W2-\-  W3)  ’ 


where  w  =  [wi  W2  is  the  frequency  variable.  Therefore,  hi(7j)  =  FT~ ^ ^ } 
is  the  desired  function  except  that  h{(rf)  is  in  £2(^1)  and  not  necessarily  in  However, 

ki(Tj)  can  be  arbitrary  closely  approximated  by  functions  in  since  is  dense  in 

£2(^1)  [39]. 

By  combining  the  finite-dimensional  state-space  equations  (5.64)-  (5.65)  with  the  infinite¬ 
dimensional  state-space  equation  (5.67),  we  will  obtain  a  single  infinite-dimensional  state- 
space  representation  for  the  system  relating  the  control  input  u{t)  to  the  canceling  field  Xp{t). 
Note  that  the  combined  system  is  obtained  by  cascading  the  finite- dimensional  system  with 
the  infinite-dimensional  system.  Specifically,  the  combined  system  can  be  expressed  in 
state-space  form  as: 

^Xp;(0  =  ApiXpi{t)  -i-  Bpiu{t),  (5.72) 

where  Xpi{t)  is  the  state  vector  of  the  combined  system 


(5.73) 
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the  linear  operator  Api  is  defined  as 


Api 


XI 


Ap  BCi 
0  Ai 


*1 

• 

Xp 

Xi 

ViApi)  = 


XI 


\xp  €  T>(Ap)  and  xi  € 


; 


(5.74) 


(5.75) 


and  Bpi  is  a  bounded  linear  operator  from  into  (H  x  R^)  defined  as 


Bpiy  = 


yyeR\ 


(5.76) 


[Biy  ^ 

It  remains  to  be  shown  that  operator  Api  is  an  infinitesimal  generator  of  a  strongly  contin¬ 
uous  semigroup  in  an  appropriate  Hilbert  space.  The  appropriate  Hilbert  space  in  this  case 
is 


H'  =  H  X  R^, 


(5.77) 


where  H  is  defined  in  equation  (5.19),  and  the  inner  product  on  H'  is  defined  as 

=  [xp,yp]H  +  [xi,yi]Rn.  (5.78) 


“ 

‘ 

*■ 

Xp 

Vp 

Xi 

yi 

Let  Ti{t)  denote  the  strongly  continuous  semigroup  generated  by  A/,  and  let  Tp{t)  denote 
the  strongly  continuous  semigroup  generated  by  Ap.  We  claim  that  operator  Api  is  the 
infinitesimal  generator  of  the  strongly  continuous  semigroup  Tpi{i)  defined  as 


Tpiii) 

_ 

Tp{t)xp  +  /q  Tp{t  -  s)BCiTi(s)xids 

V 

Xp 

Xl 

Ti{t)xi 

XI 

(5.79) 

Note  that  for  a  fixed  t  and  a  fixed  xi,  the  integral  in  (5.79)  is  a  Riemann  integral,  since 
||Tp(t  —  s)BC/Ti(s)l|  is  bounded  and  Tp{t  —  s)BCiTi{s)xi  is  a  continuous  function  of  s. 
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First,  we  show  that  Tpi{t)  is  a  strongly  continuous  semigroup  on  Hilbert  space  H',  and 
then  we  show  that  operator  Apt  is  the  infinitesimal  generator  of  this  semigroup. 

Clearly,  rp/(0)  is  the  identity  operator,  and  it  is  easy  to  show  that  limt_o+  l|rp/(0x-x||  = 
0  for  all  X  in  H'.  To  show  that  Tpi{ti  +  ^2)  =  Tpi{t2)Tpi{t\)  =  Tpi{ti)Tpi{t2).,  we  observe  that 


for  any 


in  H'  we  have 


Tpi{t2)  I  Tpi{t\) 


=Tpl{t2) 


Tp{t\)xp  +  Tp{h  -  s)BCiTi{s)xids 

Ti{ti)xi 


Tp{ti  +  t2)xp  +  Tp(t2)  Tp{ti  —  s)BCiTi{s)xids^  + 

Tp{t2  -  s)BCiTi{s)Ti{h)xids 


Tl{ti  +  t2)xi 

Tp{ti  +  t2)xp  +  Jq^  Tp{t\  +  ^2  ~  s)BCiTi[s)xids-{- 
Tp{t2  -  s)BCiTi(s)Tiih)xids 


Tl{t\  +  t2)xi 


Tp(ti  +  t2Xp  +  Tp{ti  A  t2  —  s)BCiTi{s)xids 
Ti{t\  +  t2)xi 


(5.80) 


(5.81) 


(5.82) 


(5.83) 


=?pi(^i  T  ^2) 


(5.84) 


where  in  transition  from  eq.  (5.81)  to  eq.  (5.82),  we  have  used  the  fact  that 


Tp{t2)  N  rp(s)xds  j  =  Tp{t2)Tp{s)xds. 


This  fact  can  be  proved  by  noting  that  the  Tp{t2)  is  a  hnear  and  continuous  operator  and 
the  integral  is  a  Riemann  integral. 
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Similarly,  we  can  show  that  T-pi{t\)T.pi{t-z)  =  +  t2).  Hence,  Tpi{t)  is  a  strongly 

continuous  semigroup  in  the  Hilbert  space  H'.  We  now  show  that  operator  Api  is  the 


infinitesimal  generator  of  the  semigroup  Tpi{t).  To  this  end,  take  any 


Xp 


in  H'  and 


xt 


observe  that 


rw(A) 


lim 

A-*o+ 


" 

Xp 

Xp 

XI 

XI 

lim 

A-»0+ 


+  i  jA  _  s)BCiTi{s)xids 

Ti(A)xt-xi 

A 


Axp  +  BCiXi 
Aixi 


A  BCi 
0  Ai 


(5.85) 


(5.86) 


XI 


(5.87) 


Combined  State-Space  Representation  for  the  Unwanted  Noise  Field  and  the 
Canceling  Field 


A  single  state  representation  for  the  acoustic  field  in  is  derived  in  this  Section  by  combining 
the  state-space  representation  for  the  unwanted  noise  field  and  the  state-space  representa¬ 
tion  for  the  canceling  field. 

We  propose  the  following  state-space  representation  for  the  acoustic  field  in  fl; 

=  AcXcit)  -t-  Bc{t)u{t)  -1-  Fcw{t),  (5.88) 


where  a:c(i)  is  the  state  vector 


Xc{t) 


Xd{i) 


(5.89) 
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We  claim  that  Tdt)  is  a  strongly  continuous  semigroup  on  He-,  and  that  operator  Ac  is  the 
infinitesimal  generator  of  this  semigroup. 

Using  the  fact  that  Td{t)  and  Tpi{t)  are  both  strongly  continuous  semigroups,  it  is  easy 
to  show  that  that  Tc{t)  is  a  strongly  continuous  semigroup  on  He-  To  show  that  Ac  is  the 
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In  a  practical  active  noise  cancellation  system,  the  observed  signal  y{t)  is  most  likely  the 
output  of  a  microphone.  A  simple  model  relating  the  output  of  a  pressure  microphone  to 
the  acoustic  field  surrounding  it  is 

m{t)=  f  p{t,T])g{r])dT]+  v{t),  (5.98) 

JQ 

where  m{t)  is  the  output  of  the  microphone  at  time  t;  p{t,  rj)  is  the  pressure  field  in  Q  at 
time  t;  g{r])  is  a  weighting  function  in  and  v(t)  is  white  measurement  noise,  uncor¬ 

related  with  all  other  processes  involved.  The  region  of  support  of  g{T])  is  the  surface  of  the 
diaphragm  of  the  pressure  microphone.  In  this  model,  the  output  of  the  microphone  with¬ 
out  the  measurement  noise)  is  proportional  to  the  weighted  integral  of  the  pressure  exerted 
at  the  diaphragm  of  the  microphone.  This  model  implicitly  assumes  that  the  presence  of 
the  microphone  does  not  alter  the  acoustic  field  in  fi,  and  that  the  frequency  response  of 
the  microphone  is  constant  over  those  frequencies  for  which  Jq  p(t,  T))g{T})drj  has  significant 
energy. 
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Recalling  that  the  first  component  of  Xd{t)  is  the  unwanted  pressure  field  at  time  t,  and 
that  the  first  component  of  the  Xpi{t)  is  the  canceling  pressure  field,  we  see  that  y{t)  can 
be  expressed  in  terms  of  the  state- vector  Xc{t)  as 


y{t)  =  CcXcit)  +  Gcw{t) 


(5.99) 


=  a 


-f  Gcw{t), 


(5.100) 


®p/(0 


where  Cc  is  a  linear  bounded  operator  form  H'  into  R^,  and  Gc  is  a  bounded  linear  operator 
form  H  into  with 

GcG*  =  Identity,  (5.101) 


Fcg:  =  0. 


(5.102) 


Note  that  eq.  (5.101)  ensures  that  the  measurement  noise  Gcw{t)  is  white,  and  eq.  (5.102) 
ensures  that  the  measurement  noise  GcW{t)  is  uncorrelated  with  the  process  noise  Fcw{t). 
The  explicit  operator  Cc  that  results  from  the  microphone  model  in  eq.  (5.98)  is 


Cc  =  [(a;d  +  4)’5'(^)]^'2  ^ 


(5.103) 


where  Xd  £  H  and  Xpi  G  H';  and  x^  is  the  first  component  of  Xd  and  x^i  is  the  first 
component  of  Xpi.  So  defined,  operator  Cc  is  clearly  linear  and  bounded,  since  the  inner 
product  operation  is  linear  and  bounded. 

Optimal  Feedback  Controller 

The  governing  equations  of  the  acoustic  system  in  Fig.  5.2  can  be  expressed  in  state-space 
form  as 


— Xc(0  =  AcXc{t)  -I-  Bcu{t)  -f  FcW{t) 

y{t)  =  CcXcit)  +  Gcw{t). 
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(5.104) 

(5.105) 


The  noise  cancellation  problem  is  to  generate  the  control  signal  u{t),  based  on  the  observa¬ 
tions  of  y{i)i  so  that  the  sum  of  the  acoustic  energy  in  Q,  and  the  power  expended  by  the 
control  signal  is  minimized.  Specifically,  the  cost  functional  to  be  minimized  is 

J  =  £  E{\\Q{xd{t)  +  (5.106) 

where  Q  is  bounded  linear  operator  mapping  H  into  H.  According  to  eq.  (5.41),  -f  Xp 

is  the  total  acoustic  energy  in  the  enclosure  at  time  t.  Hence,  ||Q(a;d(t) -h  Xp(t))||^  is  the 
weighted  acoustic  energy  in  the  enclosure  at  time  t. 

The  cost  function  J  in  eq.  (5.106)  can  be  expressed  in  terms  of  the  state  vector  Xc{t)  as 

J  =  E{\\Q,xS)\\Hjsdt  +  E{\\u{t)\\^]dt,  (5.107) 

where  Qc  is  a  linear  bounded  operator  from  He  into  defined  as 

^  d 

Qc  Xp  =  Qi^d  +  Xp)  Xp  ^  He,  (5.108) 

XI  XI 

with  Xd  £  H,  Xp  £  H,  and  xi  £  i2". 

From  eq.  (5.54)  we  see  that  for  this  problem,  a  unique  optimal  control  Uo(t)  exists  and 
can  be  expressed  in  feedback  form  as 

Uoit)  =  -B*K{t)xcit), 

where  the  linear  operator  K{t)  satisfies 

^[K(t)x,y]  =  -[Q*e.QeX,y]-[K(t)x,Acy]-[K{t)AcX,y]  (5.109) 

+[K{t)BeB:K{t)x,y]  K{T)  =  0,  Vx,j/eP(A,), 

and 

Xc{i,w)  =  E{xc{t,w)\y{s,w)  0<s<t). 

Furthermore,  Xc{t,  w)  is  calculated  using  a  Kalman  filter,  i.e. 

w)  =  AcXc{t,  w)  +  Pf{t)C*{y{t,  w)  -  CcXcit,  w))  +  Bcu{t,  w)  xdO,  w)  =  0, 


143 


where 

P}{t)  =  w)  -  Xc{t,  w))(xc(<,  w)  -  ic(i,  w)y), 

and  Pfit)  satisfies 

|[P/(t)x,  y]  =  [Pfit)A:x,  y]  +  [Py  (t)x,  A^y]  +  [P.P^x,  y]  (5.110) 

-[Pfit)C:CcPf{t)x,y]  Py(0)  =  A  x,ye  ViAl). 

An  identical  derivation  can  be  carried  out  for  the  case  in  which  the  boundary  F  of  the 
enclosure  is  purely  absorbing,  i.e.  eq.  (5.16)  is  satisfied  at  the  boundary  in  Fig.  5.2. 

5.7  Summary 

The  distributed  active  noise  cancellation  problem  was  studied  in  this  chapter.  We  focused 
on  minimizing  the  total  acoustic  energy  in  an  enclosure.  To  minimize  this  acoustic  energy, 
the  interaction  between  the  primary  field  and  the  canceling  field  at  every  point  in  the 
enclosure  needed  to  be  considered.  The  acoustic  wave  equation  was  therefore  utilized  to 
model  this  interaction.  The  resulting  governing  equations  for  the  acoustic  system  were  put 
in  the  frame  work  of  infinite-dimensional  stochastic  linear  systems.  Using  this  frame  work, 
we  then  derived  an  optimal  feedback  controller  for  attenuating  the  total  acoustic  energy  in 
the  enclosure. 

The  main  contribution  of  this  chapter  is  showing  concretely  that  under  certain  assump¬ 
tions  the  distributed  ANC  problem  can  be  put  in  the  frame  work  of  distributed  linear 
control.  Based  on  physical  models  for  the  canceling  loudspeaker  and  the  acoustic  field  in 
the  enclosure,  a  set  of  coupled  ordinary  and  partial  differential  equations  relating  the  con¬ 
trol  signal  to  the  canceling  field  were  derived.  We  then  showed  that  these  equations  can 
be  represented  as  abstract  differential  equations  in  an  appropriate  Hilbert  space.  This  rep¬ 
resentation  was  derived  by  proving  that  certain  operators  were  infinitesimal  generators  of 
strongly  continuous  semigroups. 

The  primary  field  was  modeled  as  the  state  of  an  infinite- dimensional  linear  system 
driven  by  a  white  process.  The  residual  field  in  the  enclosure  was  modeled  as  the  sum  of 
this  primary  field  and  the  canceling  field.  Based  on  a  physical  model  for  the  measurement 
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microphone,  we  showed  that  the  microphone  output  could  be  expressed  as  a  linear  function 
of  the  residual  field. 

We  showed  that  the  distributed  active  noise  cancellation  problem  in  an  enclosure  with 
perfectly  reflecting  or  perfectly  absorbing  boundaries  can  be  put  in  the  above  frame  work, 
provided  that  the  canceling  loudspeaker  is  modeled  as  a  collection  of  point  sources  of  pres¬ 
sure,  and  the  microphones  are  modeled  by  eq.  (5.98).  In  our  formulation,  the  control 
signal  is  pointwise  and  is  used  to  drive  an  ordinary  loudspeaker,  and  the  observed  signal  is 
pointwise  and  is  obtained  as  the  output  of  an  ordinary  pressure  microphone.  Specifically, 
our  formulation  was  used  to  find  the  feedback  controller  that  minimizes  the  distributed 
performance  criterion  J  that  was  defined  in  eq.  (5.106).  More  precisely,  this  performance 
criterion  involves  the  sum  of  the  acoustic  energy  in  the  enclosure  and  the  power  expended 
by  the  control  signal. 

The  formulation  in  this  chapter  is  the  first  mathematically  rigorous  formulation  of  the 
distributed  active  noise  cancellation  problem  known  to  us.  However,  this  formulation  needs 
to  be  extended  to  the  case  in  which  the  boundary  conditions  are  more  general. 


Chapter  6 


Conclusions  and  Future  Directions 


This  thesis  dealt  with  the  design  of  optimal  feedback  controllers  for  active  noise  cancellation. 
We  applied  optimal  feedback  control  theory  and  considered  pointwise  noise  cancellation  as 
well  as  distributed  noise  cancellation.  By  formulating  the  ANC  problem  as  an  optimal  feed¬ 
back  control  problem,  we  developed  a  single  approach  for  designing  controllers  for  pointwise 
and  distributed  active  noise  cancellation. 

The  key  strategy  is  to  model  the  residual  signal/field  as  the  sum  of  the  outputs  of  two 
linear  systems.  The  unwanted  noise  signal/field  is  modeled  as  the  output  of  a  linear  system 
driven  by  a  white  process.  Similarly,  the  canceling  signal  is  modeled  as  the  output  of  another 
linear  system  driven  by  the  control  signal.  Finally,  the  residual  signal/field  is  modeled  as 
the  sum  of  the  outputs  of  these  two  linear  systems.  These  two  models  are  combined  to 
represent  the  overall  acoustic  system  as  a  single  linear  system  driven  by  a  control  process 
and  a  white  noise  process.  We  showed  that  the  optimal  control  signal  is  a  linear  feedback 
of  the  estimated  states  of  this  overall  system.  These  state  estimates  are  computed  using  a 
Kalman  filter,  and  the  feedback  gain  is  obtained  from  a  Riccati  equation.  Note  that  in  the 
pointwise  case  and  the  distributed  case,  the  control  signal  is  not  distributed  and  is  used  to 
drive  an  ordinary  loudspeaker.  Moreover,  in  both  these  cases,  the  control  signal  is  generated 
based  on  the  outputs  of  ordinary  microphones. 

While  we  focused  specifically  on  the  acoustic  noise  cancellation  problem,  the  results 
developed  in  this  thesis  can  be  applied  to  other  active  cancellation  problems.  Vibration 
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control  is  an  example  of  a  non-acoustic  problem  to  which  our  results  can  be  applied. 


6.1  Summary  and  Contributions 

We  started  by  developing  an  optimal  feedback  controller  for  attenuating  the  acoustic  pres¬ 
sure  at  the  locations  of  a  finite  number  of  microphones.  The  error  microphones  are  placed 
at  those  locations  where  noise  attenuation  is  desired,  and  the  canceling  loudspeaker  is 
placed  in  the  vicinity  of  these  microphones.  We  model  the  residual  signals  at  the  micro¬ 
phones  as  the  sum  of  the  outputs  of  two  finite-dimensional  linear  systems.  The  unwanted 
noise  signals  at  the  microphones  are  modeled  as  the  outputs  of  a  multiple-input/multiple- 
output  finite-dimensional  linear  system  driven  by  a  white  process.  This  finite-dimensional 
system  provides  a  compact  representation  for  the  second-order  statistics  of  the  unwanted 
signals.  Similarly,  the  canceling  signals  at  the  microphones  are  modeled  as  the  outputs  of  a 
single-input /multiple-output  finite-dimensional  linear  system  driven  by  the  control  signal. 
This  linear  system  corresponds  to  the  matrix  of  transfer  functions  from  the  input  to  the 
loudspeaker  to  the  canceling  signals  at  microphones.  Finally,  the  residual  signals  at  the 
microphones  are  modeled  as  the  sum  of  the  outputs  of  these  two  finite-dimensional  linear 
systems.  These  two  linear  system  models  are  then  combined  to  represent  the  overall  acous¬ 
tic  system  as  a  single  finite-dimensional  linear  system  driven  by  a  control  process  and  a 
white  noise  process. 

We  showed  that  a  certain  class  of  performance  criterion  involving  the  mean-square  sum 
of  the  residual  signals  and  the  mean-square  sum  of  the  control  signal  could  be  expressed  as  a 
quadratic  form  in  the  state  of  the  overall  system  and  in  the  control  signal.  The  control  signal 
that  minimizes  this  performance  criterion  is  found  to  be  a  linear  feedback  of  the  estimated 
states  of  the  overall  system.  These  state  estimates  can  be  computed  using  a  Kalman  filter, 
and  the  control  gain  can  be  computed  by  iterating  a  Riccati  matrix  equation.  We  also 
showed  that  the  same  formulation  can  be  used  to  minimize  the  frequency  weighted  power 
of  the  residual  signal.  Such  a  frequency  dependent  weighting  might  be  important  because 
of  the  difference  in  the  sensitivity  of  the  ear  to  sounds  at  different  frequencies. 

Our  pointwise  feedback  controller  has  a  few  important  advantages  compared  to  existing 
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feedback  controllers  for  ANC.  Recall  that  analog  feedback  controllers  (e.g.,  [51])  do  not  take 
advantage  of  the  statistics  of  the  noise.  However,  our  feedback  controller  takes  full  advantage 
of  the  statistical  characteristics  of  the  noise.  The  discrete-time  single-microphone  feedback 
controller  proposed  by  Graupe  [26]  can  not  be  used  in  the  case  that  the  plant  model  contains 
one  or  more  delays.  However,  the  feedback  controller  developed  in  this  thesis  can  be  used 
in  this  case.  Furthermore,  the  feedback  controller  developed  in  this  thesis  can  be  used  in 
the  single-microphone  as  well  as  multiple-microphone  configuration. 

Our  pointwise  feedback  controller  achieves  the  same  acoustic  objective  as  the  existing 
multiple- channel  feedforward  controllers  for  ANC,  without  the  problems  associated  with 
acquiring  an  appropriate  reference  signal  that  is  commonly  experienced  by  feedforward 
controllers. 

A  few  experiments  were  presented  to  illustrate  the  performance  of  the  optimal  point- 
wise  feedback  controller  in  the  context  of  aircraft  noise.  The  optimal  feedback  controller 
outperformed  the  feedforward  controller  proposed  by  Burgess  by  about  lOdB.  Similarly, 
the  hardware  implementation  of  the  optimal  feedback  controller  outperformed  an  analog 
feedback  controller  by  about  9dB. 

In  the  context  of  distributed  active  noise  cancellation,  we  focused  on  developing  an 
optimal  feedback  controller  for  attenuating  the  total  acoustic  energy  in  an  enclosure.  The 
distributed  ANC  system  utilizes  one  microphone  to  monitor  the  residual  field  and  one 
loudspeaker  to  generate  the  canceling  field.  The  control  signal  driving  the  loudspeaker 
must  be  generated  based  solely  on  the  microphone  output.  To  minimize  the  total  acoustic 
energy  in  the  enclosure,  the  interaction  between  the  canceling  field  and  the  unwanted  noise 
field  at  every  point  in  the  enclosure  needs  to  be  considered. 

The  residual  field  in  the  enclosure  is  modeled  as  the  sum  of  the  outputs  of  two  infinite¬ 
dimensional  linear  systems.  The  unwanted  noise  field  is  modeled  as  the  output  of  an  infinite¬ 
dimensional  linear  system  driven  by  a  white  process.  This  infinite- dimensional  system 
provides  a  compact  representation  for  the  second-order  statistics  of  the  unwanted  noise 
field.  Similarly,  the  canceling  field  in  the  enclosure  is  modeled  as  the  output  of  another 
infinite-dimensional  linear  system  driven  by  the  input  to  the  loudspeaker.  We  showed  that 
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this  model  can  accommodate  the  finite-dimensional  dynamics  of  the  loudspeaker  as  well  as 
the  infinite-dimensional  acoustical  dynamics  of  the  enclosure,  provided  that  the  loudspeaker 
is  modeled  as  a  collection  of  point  sources  of  pressure  and  the  boundary  surrounding  the 
enclosure  is  perfectly  absorbing  or  perfectly  reflecting.  This  was  shown  by  proving  that  a 
certain  operate  was  an  infinitesimal  generator  of  a  strongly  continuous  semigroup. 

We  then  illustrated  that  the  linear  system  representation  of  the  unwanted  noise  field 
and  the  linear  system  representation  of  the  canceling  field  can  be  combined  to  represent 
the  overall  acoustic  system  as  a  single  infinite-dimensional  linear  system  driven  by  a  control 
process  and  a  white  noise  process.  Again,  this  was  illustrated  by  proving  that  a  certain 
operator  was  an  infinitesimal  generator  of  a  strongly  continuous  semigroup. 

Next,  we  showed  that  the  microphone  output  at  time  “t”  can  be  expressed  as  a  linear 
functional  of  the  state  of  the  overall  system  at  time  “t”.  This  was  shown  based  on  a  physical 
model  for  a  pressure  microphone  suggested  in  [6]. 

Finally,  we  considered  a  class  of  performance  criterion  involving  the  sum  of  the  total 
acoustic  energy  in  the  enclosure  and  the  power  expended  by  the  control  signal.  We  demon¬ 
strated  that  this  performance  criterion  can  be  expressed  as  the  expected  value  of  a  quadratic 
form  in  the  state  of  the  overall  system  and  in  the  control  signal.  The  control  input  that 
minimizes  this  performance  criterion  is  a  linear  feedback  of  the  estimated  states  of  the  over¬ 
all  system.  These  state  estimates  can  be  computed  using  a  Kalman  filter,  and  the  feedback 
gain  can  be  computed  form  a  Riccati  operator  equation. 

Our  results  for  the  distributed  ANC  problem  are  important  because  they  are  the  first 
mathematically  rigorous  formulation  of  the  distributed  ANC  problem.  Furthermore,  it  has 
been  demonstrated  in  numerical  studies  that  attenuation  of  noise  at  a  finite  number  of 
spatial  points  can  result  in  amplification  of  the  noise  at  other  points  [15]. 

Comparing  our  formulation  of  the  pointwise  ANC  problem  to  our  formulation  of  the 
distributed  ANC  problem,  we  see  the  fundamental  importance  of  being  able  to  represent 
the  acoustic  system  as  a  linear  system  driven  by  a  control  process  and  a  white  noise  pro¬ 
cess.  Similarly,  to  use  this  formulation,  it  is  essential  that  the  performance  criterion  be 
an  expected  value  of  a  quadratic  form  in  the  state  of  this  linear  system  and  in  the  control 
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process.  In  the  pointwise  case,  it  was  straight  forward  to  represent  the  acoustic  system  in 
this  form;  however,  in  the  distributed  case,  rather  sophisticated  mathematical  machinery 
was  needed  to  represent  the  overall  acoustic  system  in  this  form. 

Several  recursive/adaptive  algorithms  for  modeling  the  unwanted  noise  at  a  single  micro¬ 
phone  were  developed  in  Chapter  4.  We  assumed  that  the  unwanted  noise  at  the  microphone 
is  the  output  of  an  all-pole  transfer  function  driven  by  a  white  process  and  develop  a  re¬ 
cursive/adaptive  algorithm  for  estimating  the  parameters  of  this  transfer  function  based 
on  the  measurements  made  by  the  microphone.  An  adaptive  feedback  ANC  system  can 
be  obtained  by  combining  the  pointwise  optimal  control  law  with  these  recursive/adaptive 
algorithms  for  modeling  the  unwanted  noise.  Note  that  the  resulting  adaptive  controller 
only  adapts  to  changes  in  the  noise  statistics  and  not  to  changes  in  the  plant.  Although  we 
focused  specifically  on  developing  estimation  algorithms  for  modeling  the  unwanted  noise 
in  the  ANC  problem,  the  algorithms  presented  in  Chapter  4  can  be  applied  to  the  more 
general  problem  of  identifying  non-stationary  autoregressive  (AR)  processes  embedded  in 
white  noise. 


6.2  Future  Directions 

We  foresee  several  immediate  branches  down  which  the  work  in  this  thesis  may  continue. 

Developing  adaptive  optimal  feedback  controllers  for  pointwise  ANC  is  a  potentially 
promising  area  for  future  work.  Our  derivation  of  the  optimal  pointwise  feedback  con¬ 
troller  was  based  on  the  assumption  that  the  noise  model  and  the  plan  model  were  known. 
However,  in  some  applications  both  the  plant  characteristics  and  the  unwanted  noise  char¬ 
acteristics  might  be  unknown  or  changing.  By  properly  combining  an  adaptive  parameter 
estimation  algorithm  with  the  control  law  developed  in  Chapter  3,  an  adaptive  optimal 
feedback  controller  for  ANC  can  be  obtained.  Adaptive  optimal  control  is  a  mature  field, 
and  many  of  the  results  developed  in  this  field  should  be  directly  applicable  here. 

Our  results  for  the  distributed  ANC  problem  were  derived  for  enclosures  with  perfectly 
reflecting  or  perfectly  absorbing  boundaries.  Extending  these  results  so  that  they  apply 
to  enclosures  with  more  general  boundary  conditions  would  be  very  useful.  In  the  more 
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general  case,  the  normal  component  of  particle  velocity  at  the  boundary  and  the  pressure 
at  the  boundary  are  related  through  an  ordinary  differential  equation.  It  needs  to  be  shown 
that  operator 


0  -c  V. 
-c  V  0 


(6.1) 


,  with  domain  consistent  with  the  above  boundary  conditions,  is  an  infinitesimal  generator 
of  a  strongly  continuous  semigroup  in  an  appropriate  Hilbert  space. 

A  computationally  efficient  method  for  implementing  the  optimal  control  law  in  the 
distributed  case  is  lacking.  More  generally,  we  feel  that  a  careful  study  of  the  computational 
aspects  of  the  control  law  in  the  distributed  case  needs  to  be  carried  out. 
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